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Abstract 

The measurements of the hadronic final state in deep inelastic scattering at HERA 
are reviewed and discussed in the context of QCD. Covered are the general event 
properties in terms of energy flows, charged particle production, and charm and 
strangeness production. Quark fragmentation properties are studied in the Breit 
frame. Event shape measurements allow "power corrections" to be applied and 
the strong coupling a s to be extracted. Other a s measurements are based on dijet 
rates. Jet rates as well as charm production have been used to determine the gluon 
density in the proton. Indications have been found in the hadronic final state for 
unconventional, non-DGLAP evolution at small x, which could be explained with 
BFKL evolution. Signatures for QCD instanton effects are discussed and first search 
results are presented. 
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Chapter 1 
Introduction 

1.1 Overview 

Motivation 

HERAQ is the world's first electronQ-proton collider. The large centre of mass (CM) 
energy of 300 GeV allows to explore new regimes: new particles with masses up 
to 300 GeV can be produced, the structure of the proton can be studied with a 
resolving power varying over 5 orders of magnitude down to dimensions of 10~ 18 m, 
and partons with very small fractional proton momenta (Bjorken x down to 10~ 6 ) 
become experimentally accessible. The hadronic final state which emerges when the 
proton breaks up has an invariant mass up to 300 GeV. It provides a laboratory 
to study quantum chromodynamics (QCD) under varying experimental conditions, 
which can be controlled by measuring the scattered electron. The hadronic final 
state carries information on the structure of the proton. It is made use of to study 
the dynamics of the proton's constituents, complementary to structure function 
measurements. 

In contrast to "clean" e + e~ interactions, the initial state in ep collisions contains 
already a strongly interacting particle. That makes the physics more complicated, 
but leads also to interesting effects that cannot be studied in e + e~ collisions. At 
HERA, well established fields are being studied, exploiting the tunable kinematic 
conditions. These comprise jet physics and comparisons with perturbative QCD to 
extract the strong coupling a s and the density of gluons in the proton, or for example 
the measurement of fragmentation functions. Other topics have only blossomed with 
the advent of HERA. To name a few, a class of events with large rapidity gaps and 
small x physics allow longstanding problems in QCD to be addressed, which are 
connected with scattering cross sections at high energies and confinement. There 
are also exotic effects like instanton induced reactions, which, if discovered at HERA, 



1 "Hadron-Elektron-Ring-Anlage" 

2 HERA can operate with either electrons or positrons. In the following, the generic name 
electron is used for electrons as well as for positrons. 



would alter significantly our view of particle physics. 

Purpose 

This work provides a review of the hadronic final state measurements at HERA in 
deep inelastic scattering (DIS). The emphasis is on experimental results, because in 
many cases at HERA, experiment is driving theory. Many measurements are being 
performed without any theoretical prediction, and often they have not yet found an 
unambiguous theoretical explanation. Nevertheless, the results are discussed in the 
context of the theory, where possible. 

The review of the experimental situation is complete up to the fall of 1997. It 
can therefore be consulted for quick access to the HERA data. Basic concepts are 
explained in a partly pedagogical fashion to serve physicists from outside the HERA 
community and newcomers to HERA physics. 

Contents 

In chapter |1| the HERA machine and the HI and ZEUS experiments are introduced, 
and the kinematic variables are discussed. In chapter |2] the theoretical framework 
of deep inelastic scattering is set up with evolution equations and a special section 
on the interest in small x physics. The hadronic final state should not be discussed 
without knowledge of the inclusive ep cross section and the proton structure func- 
tions (chapter |3|). In chapter [| we start with simple models for hadron production. 
They are subsequently being refined in chapter |] and serve as the basis for the 
discussion of the data. 

Measurements of basic event properties are presented in chapter [5]: energy flows, 
charged particle spectra, charm and strangeness contents, Bose-Einstein correla- 
tions. The fragmentation of the scattered quark is studied in chapter ^ and compared 
with quark fragmentation in e + e _ annihilation and QCD calculations. Measure- 
ments of event shape variables allow a new view on hadronziation properties with 
"power corrections", offering a potentially powerful tool for measurements of the 
strong coupling a s . Jet production has been compared to perturbative QCD predic- 
tions to measure a s and the gluon density in the proton (chapter 0). By measuring 
jet rates, regions of phase space have been identified where the measured jet rates 
are not well understood yet, and where the underlying physics is possibly departing 
from the conventional picture of deep inelastic scattering. 

Also the energy flow measurements at small x have not yet found an unambiguous 
theoretical interpretation. In chapter [5] on low x physics dedicated searches for "foot- 
prints" of new QCD effects (BFKL) are discussed: energy flows, high p T particles 
and "forward jets". Chapter || deals with the possibility to discover QCD instantons 
at HERA, which would have far reaching consequences for our understanding of field 
theories and for cosmology. 

Some related topics and neighbouring fields could only be touched upon. The 
interested reader is referred to other reviews on structure functions fll], ||, rapidity 



gaps [Q, photoproduction [[| and on hadron production at fixed target experiments 

§• 

Throughout, unless stated otherwise, the data shown have been corrected for de- 
tector effects and QED radiation, and the errors comprise statistical and systematic 
errors added in quadrature. 

1.2 The HERA Machine 

In the HERA machine, electrons and and protons are accelerated and stored in two 
separate rings. The circumference of the machine is 6.3 km. The magnets of the 
proton ring are superconducting, the magnets of the electron ring are conventional. 
The final beam energies are E e = 27.5 GeV for electrons and E p = 820 GeV for 



protons with a collision centre of mass energy of ^fs = JAE e E p = 300 GeV. Early 
data were collected with E e = 26.7 GeV. 

The beams are collided head-on in two interaction regions occupied by the exper- 
iments HI and ZEUS. There are 220 bunch positions in the beam, of which typically 
190 are filled with a few 10 10 particles per bunch. The time between bunch crossings 
is 96 ns. The longitudinal bunch length is about 60 cm, leading to an approximately 
Gaussian distribution of interaction points along the beam line with width 10 cm. 
The transverse beam size is 300 /mi horizontally by 70 /zm vertically. 

In 1997 the average peak luminosity was 8.4 ■ 10 30 cm~ 2 s~ 1 with average beam 
currents at the beginning of a fill of 77 mA for protons and 36 mA for electrons. The 
total integrated luminosity for the 1997 run was 35 pb _1 . Most analyses in this 
review are based on data from the runs 1992-1994, corresponding to an integrated 
luminosity of (9(1 — 2) pb _1 . 

1.3 The HI and ZEUS Detectors 

The HI [|, and ZEUS detectors || serve to detect the scattered electron in ep col- 
lisions and to measure the emerging hadrons. The individual detector components 
are mounted concentrically around the beam line. Due to the asymmetric beam 
energies, the hadronic system is boosted into the proton direction (+z). Therefore 
the detectors are also asymmetric with respect to the interaction point, with en- 
hanced instrumentation for hadrons in the +z (forward) direction. The acceptances 
and resolutions of the main detector components for the analyses presented here are 

an 



given in table |1.1| . Fig. [TTT] shows a drawing of the ZEUS detector, and fig. fL2 
event display from HI. 

Closest to the beam line are wire chambers for measuring charged particle tra- 
jectories. The particles' momenta are determined from their track curvature in a 
longitudinal magnetic field provided by a superconducting coil. 

Electromagnetic and hadronic showers are measured in calorimeters surround- 
ing the inner tracking devices. HI emphasizes electron detection with a finely seg- 




' : 



Figure 1.1: The ZEUS detector. Protons come towards the observer. Shown arc central, for- 
ward, backward and vertex tracking detectors (CTD, FDET, RTD, VXD), the uranium calorime- 
ter (BCAL, RCAL, FCAL), the muon system (RMUON, BMUON, FMUON) and the backing 
calorimeter (BAC). The dimension of the whole detector is roughly 10 x 10 x 18 m 3 . 



merited lead (inner electromagnetic part) and steel (outer hadronic part) liquid argon 
calorimeter (LAr) with good energy resolution for electrons, supplemented by a ded- 
icated electromagnetic backward calorimeter. From 1992-1994 an electromagnetic 
lead/scintillator sandwich calorimeter was installed in the backward region (BEMC), 
and from 1995 onwards a lead/scintillating fibre calorimeter (SPACAL). A copper 
calorimeter with silicon readout (PLUG) covers part of the forward beam hole. 
Compensation for the different response to hadronic and electromagnetic showers 
in the LAr is done offline by a software weighting technique. ZEUS achieves better 
hadronic energy resolution with a self-compensating uranium/scintillator calorime- 
ter (U), and compromises on electromagnetic energy resolution. 

The calorimeters are surrounded by chambers and absorber plates for measuring 
shower leakage and for muon detection. Further specialized detectors are installed 
very close to the beam line to detect particles that are scattered under small angles. 
Silicon detectors that have already been installed as vertex detectors or are being 
planned have not yet been used in physics analyses. 



jet 



remnant 




BEMC 



Figure 1.2: DIS event with the scattered electron e' and two well separated jets detected in the 
HI detector. The proton remnant leaves mostly undetected in the +z direction. Shown are the 
central and forward tracking chambers (CT,FT), the electromagnetic and hadronic sections of the 
liquid argon calorimeter (EMC, HAC), and the backward electromagnetic calorimeter (BEMC). 
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acceptance 
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electromagnetic 


4° - 153° (LAr) 
155.5° - 174.5° (BEMC) 
151° - 177.5° (SPACAL) 


o.ii/Ve 

0.1/VE 
0.075/VE 


2.2° - 176.5° (U) 


0.18/VE 


hadronic 


0.7° - 3.3° (PLUG) 

4° - 153° (LAr) 

153° - 178° (SPACAL) 


« 1.5/ VI 

0.5/VE 
m 0.3/VE 


2.2° - 176.5° (U) 


0.35/VE 



Table 1.1: Acceptances and resolutions of the main detector components from HI and ZEUS. In 
the resolution formulae, energy E and transverse momentum pr are to be taken in GeV. Additional 
constant contributions to the resolutions of (9(1 — 3%) are not shown. 



1.4 Kinematics 

Definition of kinematic variables 



The kinematics of the basic ep scattering process in Fig. |1.3| can be characterized 
by any set of two Lorentz-invariants out of Q 2 ,x,y and W, which are built from 
the 4-momentum transfer q = k — k' mediated by the virtual boson and from the 4- 



momentum P of the incoming proton. The ep invariant mass squared is s = (e+P) 2 . 
These kinematic variables are then: 



Q 2 ■= -i 



'IX 



which gives the transverse resolving power of the probe with wavelength A = 1/Q 
(we set H= c = 1); 

QI (1-2) 



x :- 



2Pq 



the Bjorken scaling variable (0 < x < 1), which can be interpreted as the momentum 
fraction of the proton which is carried by the struck quark (in a frame where the 
proton is fast, and assuming the quark-parton model to be a good approximation); 



Pq 



Q 1 



Pk x(s — m 2 ) 



1.3) 



the transferred energy fraction from the electron to the proton in the proton rest 
frame (0 < y < 1); and 



x p 

the invariant mass squared of the outgoing hadronic system H . The invariant 

Pq 



v \- 



m„ 



;i-4) 



1.5) 



is rarely used at HERA. In the proton rest frame it gives the energy transfer from 
the lepton to the proton. 




Figure 1.3: Basic diagram for DIS in 0(a°) (quark parton model - QPM). 

The large CM energy gives access to kinematic regions both at very small x 
and at large Q 2 (Fig. [T|). The HERA data cover roughly Q 2 = 0.2 - 10 4 GeV 2 , 
x = lfr 5 - KT 1 and W = 40 - 300 GeV . 



o 
o 




Figure 1.4: Coverage in the kinematic plane (x, Q 2 ) of various DIS experiments. The kinematic 
boundary for HERA is given by the line y = 1. There exist also HERA data in the region 
10 4 GeV 2 < Q 2 < 10 5 GeV 2 . 

Experimental reconstruction of the event kinematics 

The kinematics can be determined either from the electron alone, or from the mea- 
sured hadronic system alone, or from a combination of both, permitting important 
systematic cross checks. The hadronic measurement relies mostly upon the calorime- 
ters. At large y the precision of the hadron method can be improved by momentum 
measurements in the trackers. 

The electron method The kinematic variables are calculated from the energy E' e 
and angle 9 e of the scattered electron (measured with respect to the proton 
direction): 



E' 



E P 



■ 2 "e 

sin — 



2 1-J/e 



P 2 Te 



• (1.6) 



The hadron method The kinematics is measured entirely with the hadronic sys- 



tern: 

Eh — PzH ~2 PxH+PyH ,-, -s 

Here En,p X H,PyH and p z n denote the 4- vector components of the hadronic 
system H, which are calculated as the 4-momentum sum over all final state 
hadrons h. Jacquet and Blondel || have shown that the contribution from 
hadrons lost in the beam pipe is insignificant. 



'Zl)l 



The Sigma method ||10|| Here the denominator of y^ is replaced with J2i(Ei~p z 
where % runs over all final state particles, including the scattered electron. 
This expression equals 2E e due to energy momentum conservation. In case 
the incident electron had radiated off photons which escape detection, the sum 
yields the true electron energy which goes into the ep interaction. This method 
relies on both electron and hadron measurements. With 

E = J £(E h -p th ), (1.8) 

h 

where the sum runs over all hadronic final state particles, the kinematic vari- 
ables can be written as 

E ^ , E' 2 sin 2 6 e , , 

T, + E' e (l -cos6» e ) 1-Z/s 

The denominator E + E' e (l — cos8 e ) is twice the energy of the "true" incident 
electron, after QED radiation from the incoming electron beam. 

The double angle method We define the angle 7 by 

Pth ~ (Eh ~ PzH? ,, in x 

cos7=^ — -. (1.10) 

Pth + ( E H -PzhY 

In the simple quark parton model 7 would be the angle of the scattered (mass- 
less) quark. The kinematic variables can be calculated from 7 regardless of its 
interpretation: 

Qi,=wi . ; in ^ 1+co 7> (lid 

sin 7 + sin v e — sm(^ e + 7) 

and 

__ (E e \ sin 7 + sin 6> e + sin(9 e +7) 

\ E p J sin 7 + sin 9 e — sm(9 e +7) 

For most of the phase space the electron method is superior. At small y the 
hadron method has a better resolution than the electron method. The "mixed 
method" uses Q 2 reconstructed from the electron method and y reconstructed with 
the hadron method. Also the double angle method and the sigma method use infor- 
mation from both the electron and the hadronic system, thus interpolating between 
the pure electron and hadron methods. The sigma method has the advantage that 
it corrects for initial state radiation. 

10 



Chapter 2 
Theoretical framework 

2.1 Deep Inelastic Scattering 

The fundamental measurement in DIS concerns the cross section for ep — > e'H 
as a function of the kinematic variables (any pair of two independent ones). The 
quark parton model (QPM) offers a physical picture: the scattering takes place via 
a virtual photon which is radiated off the scattering electron, and which couples to 
a pointlike constituent inside the proton, that is a quark or antiquark. The cross 
section is then proportional to the quark density inside the proton. 

The differential cross section ep — > e'H can be expressed in terms of twoQ inde- 
pendent structure functions Fi(x,Q 2 ) and F 2 (x,Q 2 ): 

{l-y)F 2 + y 2 xF 1 } 

i-y + ^rfs]^ (2-i) 

(l-V + *)-F 2 -£.F L . 

a is the electromagnetic coupling constant. Here we have expressed the cross section 
also in terms of the longitudinal structure function Fl(x, Q 2 ) and the ratio R, defined 
as 

F L :=F 2 -2xF 1 ^j_ = F 2 - 2xF, = a^ 

F 2 -F L 2xF 1 a T K J 

R can be interpreted as the ratio of the cross sections ay and <jl for the absorption 
of transversely and longitudinally polarized virtual photons on protons, with a 
0\l + &t ■ The structure function F 2 can be expressed^ in terms of ot and ol 



d 2 ff 


Ana 2 


cLrdQ 2 


xQ 4 




4-rra 2 




xQ 4 




Ana 2 



7 V 
tot 



Q (1 — x) Q 7 *„ Q , » , . 

4n 2 a Q 2 + 4m 2 p x 2 ' ^ " 4^ ( ^ + ° T } ' (2 ' 3) 



-'^For the sake of simplicity, Z exchange (a 1% correction for Q 2 rs 1000 GeV 2 ) has been 
neglected, and the structure function F 3 thus been omitted. 

2 We use the Hand convention |ll| for the definition of the virtual photon flux. 
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where the small x approximation has been applied. Similarly, 

Fl Q 2 (i-x) ^ ^_^. OL _ (24) 

In the "DIS" scheme F 2 can be written in terms of the quark and antiquark 
densities, qi and q i} and their couplings to the photon, i.e. their charges e qi : 

F 2 (x,Q 2 )=xJ2el [%0r,Q 2 ) + g 4 (a;,Q 2 )] , (2.5) 

i 

where the sum runs over all quark flavours^. In other schemes (for example the 



"MS"' scheme) the relation between F 2 and the parton densities (eq. |2.5|) holds 
only in leading order perturbation theory. The longitudinal structure function Fl 



vanishes in zero'th order a s , and will be discussed in section pl3 . 

In the simple quark parton model the proton consists just of 3 valence quarks. 
Their distribution functions in fractional proton momentum x, xq(x), would peak 
at x ~ 1/3 and tend towards zero for x — > 0, 1. In a static model of the proton, 
they would not depend on Q 2 . It follows that F 2 should not depend on Q 2 , just on 
x (Bjorken scaling). 

When QCD is "turned on" the quarks may radiate (and absorb) gluons, which 
in turn may split into quark - antiquark pairs or gluon pairs. More and more 
of these fluctuations can be resolved with increasingly shorter wavelength of the 
photonic probe, A = l/Q. With Q 2 increasing, we have a depletion of quarks at 
large x, and a corresponding accumulation at lower x. In addition, "sea quarks" 
from g —>■ qq splittings populate small x. In fact, at small x it is the gluon content 
with distribution function g(x, Q 2 ) that governs the proton and gives rise to the DIS 
cross section via the creation of qq pairs. 

2.2 Evolution Equations 

It has not yet been possible to calculate the structure of the hadrons from first 
principles, involving the building blocks of hadronic matter, the quarks and gluons, 
and their mutual interactions as given by QCD. Therefore also the lepton-nucleon 



3 Equation 2.5 represents the "leading twist" (called twist 2) contribution to the structure 
function F 2 , when expanded in powers of Q 2 B , 



F 2 (x,Q 2 ) = yC n (x,Q 2 )/(Q 2 ) n . (2.6) 



The coefficients C n (x,Q 2 ) are varying logarithmically with Q 2 . The "higher twist" terms (n > 0, 
called twist 4, 6 etc.) arise from interactions of the struck parton with the remnant and are 
suppressed by (l/Q 2 )™. We shall not pursue higher twist effects any further, but note that they 
may not be negligible at small x for Q 2 as large as a few GeV 2 M, O] . 



12 



scattering cross section cannot be calculated from first principles. Due to the factor- 
ization theorem of QCD one can however split the problem. The cross section can be 
calculated by folding initial parton distribution functions fi/ p , giving the density of 
partons % in the proton p, with a perturbatively calculable lepton-parton scattering 
cross section. Symbolically 

'<* = J2 [fi/p ® a ™] ( 2 - 7 ) 



(J, 



for ep scattering (see fig. |2.1| a). The initial parton distributions cannot be calculated. 
They have to be determined experimentally. They are however universal in the sense 
that once they have been measured in one reaction, they can be used for calculations 
of other processes. Rather than employing the rigorous operator product expansion 
(OPE) technique for the evolution of the structure functions (see e.g. [01), we shall 



use in this section the more intuitive picture of Feynman diagram summation. 





a) 

Figure 2.1: a) Deep inelastic ep scattering. The ep cross section is factorized into elec- 
tron - parton cross sections <r e i and parton densities fi/ p with the factorization scale /if'- 
a ep = J2ilfi/p(PF) ® a ei(np)]- b) The lowest order diagram (Born graph) contributing to cr e ,; 
in a). 

The expansion parameter for the perturbation series is the strong coupling a s . 
The coupling is scale dependent according to the "renormalization group equation" 
(see [Q for a concise summary). The two- loop expression (next-to-leading order = 
NLO) for the "running" coupling a s as a function of the renormalization scale [xr is 



where 



a s (LiR 2 ) 



0o 



An 



/3 ln(/x R 2 /A 2 ) 



11 "3»/ 



1- 



20, ln[ln(^ 2 /A 2 )] 
Pi ln(/V/A 2 ) 



Pi = 51 



19 



-Tli 



(2- 



(2.9) 



The renormalization scale fiR is set by the length scale (1///r) over which the inter- 
action takes place, given for example by the virtuality of the probing photon Q 2 in 
DIS, or by the px of a parton. By eq. |2l| the QCD scale parameter A is introduced. 
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Figure 2.2: Diagrams contributing to the DIS cross section a. 



ep • 



In contrast to a s , its definition depends on the number of active flavours rif and on 
the renormalization scheme (for example the "minimal subtraction scheme" MS). 
The strong coupling decreases with increasing scale - at short distances partons 
become asymptotically free. a s grows beyond all bounds for small scales (/zr — > A) 
or large distances, when perturbation theory breaks down and confinement sets in. 
The size of a hadron w 1 f m provides an estimate when that happens. Therefore 
A = 0(l/fm) ~ 0.2 GeV. The current world average for a s at the scale set by the 



Z mass is a s (m 2 z ) = 0.118 ± 0.003 or A- 



(5) 



0.209 



+0.039 



'MS —0.033 

The calculation of the cross section <t p? - is a formidable task 



GeV for 5 flavours IT 



It turns out that 
there is no fast convergence of the perturbation series (fig. 2^). Many diagrams 



contribute and have to be summed up. One encounters two types of divergencies. 
Divergencies due to the radiation of soft quanta with small momenta k — > are 
exactly cancelled by virtual corrections to graphs where that radiation is absent 
( "no emission" ) . Divergencies due to collinear radiation (so called collinear or mass 
singularities for kx — > 0) can be absorbed (factorized off) into the "bare" parton 
distribution functions. Thereby they are redefined and depend now on the (mass) 
factorization scale fip, and so does the electron-parton scattering cross section with 
the singularities removed: 



a ep = Ylifi/P&F) ® ^«(/*f)] 



(2.10) 



(see fig. |2.1| a). The choice of the factorization scale is arbitrary. The physical cross 
section a ep is of course independent of /ip. Often one chooses /i F = Q 2 , because 
then a e i reduces to the Born graph (here ei — > ei by 7 exchange, fig. |2.1| b). In the 
parton picture fi/ p {Q 2 ) is then interpreted as the parton density in the proton as 
seen by a photon with virtuality (resolving power) Q 2 . 

Since <J e i(fi 2 F ) can in principle be calculated perturbatively for any scale, one can 
also calculate the change of the redefined parton distribution function with a change 
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of scale. These are the evolution equations. Once a parton distribution function (or 
structure function) is known at one scale, it can be calculated for any other scale. 
For the derivation of the evolution equations, one has to perform the perturbative 
calculation of a e i, taking into account all contributing graphs (fig. |2.2| ). In order 
to carry out the calculation in practice, one applies certain approximations, thus 
restricting the phase space for radiation. Such approximations are then valid in 
regions of x and Q 2 where the selected contributions are the dominant ones. In the 
following, evolution equations will be discussed which differ in their approximations, 
and therefore in their x, Q 2 regions of validity. Always, in order to allow perturbation 
theory to be valid, a s {Q 2 ) <C 1 is required. 

In a "physical" gauge, in which only the physical transverse gluon polarization 
states need to be taken into account, the individual contributions of the perturbation 
series can be represented by so-called ladder diagrams (see fig. J2l3| ). We work in a 
frame that moves parallel to the proton, and where the proton is fast. The transverse 
momenta of the emitted quanta are denoted with pxi- Similarly, the transverse 
momenta carried by the quanta that constitute the side rails of the ladder are kxi- 
The longitudinal components are given in fractions of the proton energy Ei/E p and 
are labelled with £j for the emitted quanta and with Xi for the internal quanta. 
Energy- momentum conservation requires Xi = Xi+\ + £j, and therefore X{ > Xi+\. 




\n 



Figure 2.3: The notations for a ladder diagram with n rungs. 
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2.3 The DGLAP Equations 

In the approximation leading to the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli- 



Parisi) equations [15| all ladder diagrams are summed up, in which the transverse 
momenta along the side rails of the ladder are "strongly ordered", Ql <C •••^tj "C 
kq-i+i "C ... <C Q 2 . This condition implies strong ordering also for the emitted 
quanta, p T i < Pn+i- 

Where can we expect such an approximation to be valid? A good pedagogical 



discussion can be found in [Tj| . We shall sketch the main arguments. The evaluation 
of a ladder diagram with n rungs requires integrations over the internal momenta 
exchanged between rungs of the form 



a s ■ 



k 2 ' 



dxi 



J j 



(2.11) 



where the dots represent functions which depend on the actual nature of the emitted 
quanta and their dynamics. With strong kx ordering, the nested integration over all 
n rungs in the ladder can be carried out. The result is an expression oc (a s In %z) n . 
We see that the kx integration yields large logarithms when the kx 's are strongly 
ordered. They compensate the smallness of a s . Clearly, since a s decreases only log- 
arithmically with Q 2 and is compensated by a logarithmically growing term in Q 2 , 
in a perturbative expansion all graphs with rungs up to n = oo need to be summed 
up. (Often the expression "re-summation" is used, because one re-arranges the per- 
turbation series such that the largest terms come first). This is called a leading log 
approximation (here in 1ii(Q 2 /Qq)), since each power n in a s is accompanied by the 
same (maximal) power of \n(Q 2 /Ql). Subleading terms would be oc a"(ln §2-) n_1 . 

We expect this to be a good approximation when Q 2 is large, but x is not too 
small in order not to produce also large logarithms, 



a s (Q 2 )ln-«a s (Q 2 )ln^<l. 
x C^o 



(2.12) 



In this approximation, the evolution equations for the quark density qi for flavour 
% and the gluon density g are 



dqj{x,Q 2 ) 
dlnQ 2 

dg(x,Q 2 ) 
dlnQ 2 



a s r l dz 
2n Jx z 
a s f 1 dz 
2% Jx z 



q t (z,Q 2 )P, 



m 



x 



g(z,Q 2 )P qg 



22 qi (z,Q 2 )P gq (^j + g(z,Q 2 )P gg (- z 



(2.13) 



These are the famous DGLAP equations |15| , describing the scaling violations of the 
structure functions. They involve the calculable Altarelli-Parisi splitting functions 



Pij((). f^-fy(C) gives the probability per unit of In g^ for parton branchings q 



qg, 
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g — > gg and g — > gg, where the daughter parton i carries a fraction 1 — £ of the 
mother's (j) momentum. The splitting functions in LO are given by 

Pgq\() = 3 c = 32^ — ^ f2 14) 

P qg (0 = KC 2 + (i - C) 2 ) = P qg (i - C) 

^w(C) = 6(i^ + ¥ + C(l " 0) = ^ 9 (i - 0- 

The singularities for soft emissions ( — > 1 are cancelled by virtual corrections to 
the "no-emission" graphs. This is physical, because arbitrary soft emissions can- 
not be distinguished from the "no-emission" caseQ. The coupled integro-differential 
equations for the quark and gluon densities (eq. |2.13|) can be solved, allowing to 
calculate them for any value of Q 2 and x > Xq, once they are known at a particular 
value Ql for x > x . 

A special case for which the DGLAP equations can be solved analytically (see for 
example fi~3f ) occurs when in addition to the above conditions also strong ordering 
in x is required, x <C ---Xi+i « ii « ...Xq- The large logarithmic terms arising 
from the integration are then of the form oc (a s (Q 2 ) In g? In ^) n , which need to be 

resummed. This is the double leading log approximation (DLL). It is expected to 
hold when the DLL terms dominate over the others, 



^{■f U^h,*..!*!. (2,6) 

a s (Q 2 )\n± J Qt x 



X 

Y2 



This is the case for large Q and small x. At small x the parton content of the 
proton is expected to be dominated by gluons, because Pij(( = x/z) is largest when 
gluons are being produced (i = g). When quarks are neglected, and P gg — | is 



approximated, the DGLAP equations can be solved to yield the DLL solution fll~T 
xg(x, Q 2 ) ~ xg(x, Ql) exp 



144, 
In 



\ 25 



ln(Q 2 /A 2 ) 



\n(Q 2 /A 2 ) 



ln(l/x), (2.17) 



provided the gluon density is not too singular at small x (needs to be quantified). 
At small x a fast rise of the gluon density with decreasing x is predicted. That 
is, xg increases faster than (ln-) A , but slower than (-) A for any powers A, A' > 0. 
Apart from these shape restrictions the actual rate of the growth is not predicted, 
it depends on the "evolution length" from Qq to Q 2 . 

4 Technically the singularities can be regularized by the following prescription. Replace 1/(1 — C) 
with 1/(1 — C) + , where the "+ - prescription" defines how integrals are to be carried out: 

V «o : . r d< mzm (2 , 15) 

o I 1 — U+ Jo L s. 

for any function f(z). A term 2<5(1 — Q has to be added to P qq , and a term (11/2 — nt/3)8(l — C) 
to P gg . 
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2.4 The BFKL Equation 

When x is small, but Q 2 not large enough to reach the DLL regime, the DGLAP 
approximations cease to be valid. For the limit 1/x large and Q 2 finite and fixed the 
BFKL (Balitsky-Fadin-Kuraev-Lipatov) |18| equation has been derived. It takes into 
account diagrams in which the Xi are strongly ordered, xq <C ...x-i <C Xi+i <C ... <i. 
No ordering on kxi is imposed. Large logarithms oc (a s ln-) n are thus generated 
that need to be resummed, leading to the leading log approximation in In-. The 
region of validity is 

Q 2 ,„,,, l 



a s (Q 2 )\n^- <a s (Q 2 )ln-<l. 



Ql 
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The BFKL equation is expressed in terms of the "unintegrated" gluon density 



f(x, kp), which is related to the usual gluon density by 



xg(x, Q 2 



Q 2 ^fc 2 

Krp 



(2.19) 



The BFKL equation is an evolution equation in x. It is formulated for gluons which 
dominate at small x, 



df{x, fcfj 

<91n(l/x) 



3a s 2 

rvn-i 

TT JO 



™dk' 2 



7,/2 
fvrp 



j[x,k T ) j{x,k T ) 



\h' 2 



rCnn 



+ 



j(x, k T ) 

t: rCrp ~\ rCnn 



K®f, (2.20) 



where K is the BFKL kernel. f(x, k\) can be calculated for any (small) x, once it 
is known at some x$ for all k\. We note in passing that if one requires strong kx 
ordering (kj> 3> k^) when solving eq. 2.20 , the DLL result is retained 



For fixed a s the equation can be solved analytically. The result is (in the saddle 
point approximation) 



J{X, Krp) 



OC 



(OC 



[x/Xq) 



-\ 



27rA"ln(x /x) 

X N _A 



: CXP 



T-2 



-ln 2 (fc£/fcr) 
2A"ln(a;o/a;) 



x 



with A 



n c a 



- -4 In 2 fh 0.5 (for n c = 3 colours and a s = 0.19), and A" 



»r<-V.s 



:2.2r 



28C(3) 



T-2 



(the Riemann ( function gives £(3) = 1.202). k T specifies the starting point for the 
evolution. Therefore the gluon density is expected to rise like a power of 1/x for 
decreasing x, xg(x,Q 2 ) oc x~ x , faster than the DLL result eq. |2.17| . However, the 

PI- 



0.5 



running of a s and higher order corrections decrease the value of A 

Another characteristic prediction of the BFKL equation is u kp diffusion", in 
contrast to kp ordering for DGLAP (see fig. |2.4| ). The kp distribution function 

"evo- 



f(x, kp) / Jkp is Gaussian in In kp with a width that increases with the BFKL 



lution length" J\u.(xq/x). An individual evolution path will follow a kind of random 
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Figure 2.4: Possible evolution paths for the parton cascade a) for DGLAP evolution with strong 
kx ordering and b) for BFKL evolution without kf ordering. Shown are the Xj and kxi in the 
ladder for fixed start (xo, k T ) and end (x, Q 2 ) points. The infrared region k\ i < Q\ is to the left 
of the dotted line. The dashed "cigars" represent the ±ler contours for the {xi,kri) taken from 
a large ensemble of evolution paths. One of the BFKL paths shown "diffuses" into the infrared 
region k\ % < Q\. 



walk in k\. An ensemble of evolution paths exhibits a diffusion pattern according 
to the Gaussian in lnfc^. 

hr diffusion poses a difficulty for the application of the BFKL equation, because 
kx may diffuse into the infrared region (kx < Qo) where perturbation theory cannot 
be applied. One therefore usually introduces a lower cut-off kxo for the h? integra- 
tion, and studies the dependence of the result on that cut-off. Due to kr diffusion 
BFKL looses much of its predictive power when applied to the structure function F 2 . 
The inclusive structure function F 2 is probably not a good place to identify BFKL 
effects unambiguously. It is however possible to study special final state configura- 
tions where diffusion into the infrared region can be minimized by fixing both start 
and end point of the evolution far enough above the infrared region. The search for 
signs of BFKL evolution in the hadronic final state is presented in chapter || 

The CCFM equation |20| developed in recent years unifies the BFKL and DGLAP 
approaches [^TJ and takes into account coherence effects by angular ordering. The 
CCFM approach leads to a reduction of the exponent A and a reduction of the 
kx diffusion |32|, ^]. The Linked Dipole Chain model [21 provides an implementa- 
tion of the CCFM equation which is suited for final state predictions. 
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2.5 The Interest in Small x 



Orthogonal evolution equations 



In fig. |2l| the regions of validity of the different evolution equations are sketched. 
They are not predicted precisely by the theory, but have to be explored experimen- 
tally. 

With DGLAP evolution a parton density p(x,Qq) known for x G [xq, 1] can be 
evolved to any value of Q 2 for x G [xo, 1]. The behaviour for x < xq cannot be 
predicted with DGLAP. Similarly, with BFKL evolution a parton density p(xo, k?) 
known for < k\ < oo can be evolved to any value of x. The new feature, orthogonal 
to the DGLAP evolution, is that the low x behaviour is predicted by the theory. In 
principle DGLAP and BFKL evolution together could be used for a Munchhausen 
trick (bootstrapping), to predict the structure of the proton for all x as long as 
Q 2 > Ql, above a cut-off to avoid the non-perturbative region. The problem with 
kx diffusion into the non-perturbative regime however poses a severe obstacle for 
that goal to be reached. 
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Figure 2.5: Schematic map of the kinematic plane. Indicated are the actions and regions of 
validity of the DGLAP, DLL and BFKL evolutions. The non-perturbative region (Q 2 < Q 2 ,) and 
the saturation limit are given by the shaded areas. 

Another motivation is connected with the cross section for hadron-hadron scat- 
tering at high energies. Before we can make the point, we need to make an excursion 
to Regge theory (see for example [25] as a review related to HERA physics; |26 and 
2?J as a modern language textbook introduction; |28| for an in depth discussion of 



Regge theory). 
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Regge theory 

Consider the elastic scattering of hadrons a and b, ab — > ab (fig. |2.6| a). Their 4- 
momenta are denoted with a, b for the initial and a', b' for the final state. The cross 
section can be expressed as a function of the Mandelstam variables s = (a + b) 2 and 
t = {b — b') 2 . It is the squared sum over the scattering amplitudes due to the quanta 
X (conventionally mesons) that can be exchanged, 



£ = 5S? £**<•.*> 



A' 



(2.22) 



In Regge theory, where the mesons are connected via a so-called Regge trajectory 
(explained below), the sum yields 

EK(M) = Mffi(M)«^)e(«W)s aW . (2.23) 

x 

(3(t) is an unknown real function. The complex phase is given by 

£(«(*)) = e~ ina{t)/2 (2.24) 

for exchanged particles with C parity +1, and receives an extra factor i for negative 
C parity. The Regge trajectory a(t) gives the relationship between the mass m and 
the spin J of the exchanged mesons, J = a{m 2 ) (fig. |2^b). 

Empirically, Regge trajectories can be parametrized as straight lines with 

a{t) = a + a't. (2.25) 

a is called the intercept (with the ordinate), and a' the slope of the trajectory. As an 
example, a Regge trajectory for mesons is shown in fig. |2.6[ b. Due to the confinement 
problem, meson trajectories (hadron masses) could not yet be calculated from first 
principles in QCD. 

The total cross section 

Starting from the elastic cross section (in principle one has to sum over all Regge 
trajectories whose resonances can be exchanged in the reaction) 

^oc(/3(t)) 2 -^)- 2 , (2.26) 

we can use the optical theorem, relating the total cross section to the forward scat- 
tering amplitude 

a tot = -lmM cl {s,t = 0), (2.27) 

s 

to predict the behaviour of the total hadron-hadron scattering cross section 

a tot ex s aa -\ (2.28) 
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a) 



b' 



b) 



m 2 = t (GeV 2 ) 



Figure 2.6: a) The elastic scattering of hadrons a and b, mediated by the exchange of a quantum 
X. The 4-vectors before and after the scattering are a, b and a' , b' . b) Regge trajectories. The 
trajectories for the p, u, f% and the 02 resonances almost coincide; they are represented here with 
a solid line. The indicated resonances are the p(770), the w(782), the /2(1270), the a2(1320), the 
P3(1690) and the W3(1670). The Pomeron trajectory (dashed line) is shown together with the 
I{J PC ) = 0(2++) state X(1900) observed by WA91 g§. 

The total cross sections for pp, pp, jp and 77 reactions are plotted as a function 
of the CM energy y/s in fig. |2.7| . Their behaviour is surprisingly similar (and also for 
other hadron-hadron scattering cross sections like np, Kp etc. that are not shown 
in fig. ^J\ UTS]). They fall at small CM energy \/s < 10 GeV, and rise towards large 
energy. All these cross section can be parametrized with the universal ansatz |3(J 



a tot = As* 1 ^- 1 + Bs aip ( 0) -\ 



(2.29) 



A and B are process dependent constants, whereas ozm(0) — 0.5322 ± 0.0059 and 
«ip(0) = 1.0790 ± 0.0011 ||14[ are universal, process independent constants. 



The Pomeron 

The fall off at small energies is readily interpreted as due to meson exchange, whose 
Regge trajectories have an intercept am (0) < 1 (compare fig. p76| b; here only trajec- 
tories that dominate high energy scattering are shown; other meson trajectories have 
smaller intercepts and therefore do not contribute much to high energy scattering). 
Correspondingly, the rise at high energy is attributed to an exchange described by 
a Regge trajectory with intercept «/p(0) > 1. There exists however no established 
set of particles with such a Regge trajectory. Nevertheless, the Regge ansatz gives 
a successful parametrization of the scattering process. Therefore a hypothetical ob- 
ject, the Pomeron IP is postulated, which has the quantum numbers of the vacuum 
(electrically and colour neutral, isospin and C parity +1), and whose exchange 
is described by the Pomeron trajectory ajpit) pa 1.08 + 0.25 • t, see fig [2.6|b. The 
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Figure 2.7: The total cross section for pp (or pp), -fp and 77 interactions as a function of 
the total CM energy ^/s = Eqm PM- ZEUS and HI measure a lp in ep scattering with Q 2 i=s 



(photoproduction). The curves represent the DL parametrizations with aip 
=1.112 (dashed) and =1.088 (dotted). 



= 1.0808 (solid), 



object is suspected to be of gluonic nature, perhaps a glue ball. A possible glue 
ball candidate with J PC = 2 ++ in fact would fall on the Pomeron trajectory, see 
fig. |2.ti|b. Because ajp(0) > 1, physical states with J = 0, 1 belonging to the Pomeron 
trajectory cannot exist. An up-to-date textbook on QCD and the Pomeron is p7j . 
Deep inelastic scattering at small x can be viewed as virtual photon - proton 
scattering at high energy A /s 7 » p = W ~ JQ 2 /x. We had connected the structure 
function F 2 with the total cross section for 7*p scattering, 

„ Q 2 



4:TT 2 a 



a 



TV 
tot • 



(2.30) 



If the total cross section behaviour a tot oc s found for hadron-hadron scatter- 
ing and real photon-hadron scattering continues to hold for virtual photon-hadron 
scattering, one expects F2 to rise as F^ oc (1/x) 008 with decreasing x. 

We note that the power growth of the total cross section cr to t oc s x will eventually 



violate the Froissart bound [32] for hadron-hadron scattering 



Otot 



<^-fln^ 



mi 



•so 



X31] 
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where Sq is an unknown constant. The power growth of Otot must be dampened 
by some mechanism at large energies. 

The small x behaviour of F 2 

At small x, F 2 will be determined by the dominant gluon content of the proton, 
because the quarks to which the photon couples are pair created by the gluons. The 
BFKL prediction for small x was 

xg(x) oc I — ) =>- a7ot P oc F 2 oc xg(x) oc I — ) oc s x . (2.32) 

The (LO) BFKL expectation for the small x behaviour of F 2 is F 2 oc (1/s) 05 . 
This power growth is faster than the growth expected from eq. [2.17| , the DLL 
approximation. 

This situation is very interesting. The experience from total cross sections at 
high energies would suggest that F 2 should rise oc (1/x) 008 at small x, a behaviour 
long known and parametrized with the soft Pomeron, but whose origin is not under- 
stood from QCD. On the other hand, in the BFKL approximation, QCD does make 
a prediction for small x, which is different from the past experience: F 2 should rise 
much faster, oc (l/x) a5 . The DLL expectation is in between. The slow rise is often 
said to be due to the "soft Pomeron" or the "non-perturbative" Pomeron, or the 
"Donnachie-Landshoff" Pomeron. The fast rise would be attributed to the "hard", 
or "perturbative" , or "Lipatov" , or "BFKL" Pomeron, if one still wants to use the 
Regge language (see fig. |2.8|) . Under which conditions will we see which behaviour? 
How about the transition region? Do HERA data extend into a kinematic regime 
where the steep rise of F 2 predicted by BFKL can be seen? And if a steep rise is to 
be seen, is it really to be attributed to BFKL dynamics? Some of the answers will 
be given by the HERA data presented in the next chapter |3] on structure function 
measurements. It will turn out however that F 2 is too inclusive a quantity to re- 
solve the question of BFKL dynamics. The search for specific signatures of BFKL 
evolution in the hadronic final state is presented in chapter M. 



IP 




Figure 2.8: The exchange of a (BFKL) Pomeron in the Regge language is equivalent to a sum 
over graphs with gluon ladders between the interacting particles in perturbative QCD. 

Ultimately the hope is that the BFKL equation offers a way to approach the 
confinement problem of QCD. It allows to make predictions for F 2 at small x, thus 
for the structure of hadrons at small x, something which otherwise has to be assumed 

24 



as non-perturbative input for the parton densities. It also makes a prediction for the 
total j*p cross section at high energies, where previously we had just a prediction 
based upon a parametrization of "soft" phyics, the "soft" Pomeron. That the two 
predictions do not coincide makes it all the more interesting. 

2.6 Hadronic Final States 

So far we have developed the theory for the total inclusive cross section, loosely 
speaking a sum over everything that can happen inside the proton. When the 
proton is being probed by the virtual photon, one out of all the possible virtual 
fluctuations in the proton is selected by the measurement process. The remnants 
of the fluctuations materialize in the hadronic final state and become observable 
when the proton wave function is projected onto a specific state. Also the scattered 
quark and radiation thereoff contribute to the hadronic final state. Matters are 
even more complicated: the distinction between initial state fluctuations and final 
state radiation may be practical and justifiable in many cases, but is quantum 
mechanically not rigorous. 

QCD predictions for the hadronic final state are in general much more difficult 
and less rigorous than for the inclusive cross section. On the other hand the hadronic 
final state can provide much more detailed information on the QCD processes in ep 
scattering than just the total cross section. 

Different observables and approximations 

It depends very much on the final state observable which technique and approxi- 
mation is appropriate for a theoretical description. For example, for a global event 
property like the amount of energy on average emitted in a certain solid angle, 
it may be a good approximation to interpret the parton evolution equations in a 
probabilistic wayQ. The evolution equations define a cascade of parton emissions (a 
parton shower) with specified emission probabilities. One has to sum up all emis- 
sion cross sections, weighted with the energy. Techncially the calculation could be 
done analytically, numerically, or with a "Monte Carlo" program. Of course such a 
calculation will only be valid where the approximations are valid that went into the 
evolution equations. 

For less global quantities, one may find that one has selected a piece of the 
phase space which is unimportant for the total inclusive cross section, and which 
might have been rightfully neglected in the evolution equations. For example, for 
the process 7*g — > qq with high transverse momenta of the quarks, the DGLAP 
approximation with strong pt ordering will not be sufficient to describe the cross 



5 Strictly speaking the evolution equations have been derived for the inclusive cross section. 
Cancellations between different diagrams are necessary to ensure that the total cross section stays 
finite. It is therefore not a priori clear to what extent the evolution equations can be used for 
hadronic final state predictions. 
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section accurately. In that case one will instead use the exact, fixed order QCD 
matrix element for the process 7*g — > qq to cover the full phase space, and fold it 
with the probability to find a gluon in the proton. In practice the matrix element 
will be calculated only to a few orders in perturbation theory, in leading (LO) or 
next-to-leading order (NLO)[]. In turn, this approximation will be insufficient for our 
first example, where many parton emissions contribute. This case is hardly covered 
by the fixed order matrix element (compare section |5.1|) ; in NLO there are at most 3 
partons in the final state. The actual implementations of the different perturbative 
QCD approximations are presented in section [4.4| . 

A special role play "infrared safe" observables. An infrared safe observable does 
not change its value when an object in the final state (parton, particle, energy 
cluster, ...) is split collinearly into smaller units, or when a soft object is added. 
Examples for infrared safe observables are collective observables like energy flows and 
certain event shape and jet observables. A counter example is the total multiplicity. 
Infrared safety ensures that perturbative predictions can be made without the need 
to introduce a cut-off parameter against infrared divergencies. This property is also 
experimentally advantageous, as it does not require to know the particle composition 
of a certain energy deposition in the detector that is used in the analysis. For 
example, calorimeter clusters usually do not have a one-to-one correspondence to 
incident particles. 

Event generators 

The calculational techniques differ. Some specific final state observables can be cal- 
culated analytically (rarely), or numerically with a program. More ambitious are 



event generators, "Monte Carlo" programs (see section 4/7 ) . They attempt to model 
event by event the complete final state. Ideally, an ensemble of Monte Carlo gen- 
erated events would correspond in all aspects to an ensemble of events collected in 
nature. An intrinsic difficulty with this approach is that the probabilistic generators 
have to use probabilities, not probability amplitudes. Quantum mechanical interfer- 
ence effects therefore have to be implemented "by hand" , for example by enforcing 
an angular ordering for parton emissions (see section |4.4j ) . 

Hadronization 

One of the biggest problems of QCD is to make contact with the real world of 
hadrons. Perturbative QCD makes predictions for parton cross sections. Observ- 
able are hadrons though. The transition from partons to hadrons ( "hadronization" ) 
cannot be treated perturbatively, because it happens at a scale where the strong 
coupling constant becomes large. Non-perturbative QCD can be considered as one 
of the last frontiers, where new insight into nature can be gained, or as a nuisance, 



'Calculations in next-to next-to-leading order (NNLO) are not yet available for DIS. 
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because it hides the beauty and simplicity of the partonic world from direct obser- 
vation. 

In any case, much has been learnt about hadronization from the study of hadron 
production. This has led to some rather successful phenomenological models to de- 
scribe hadronization (see section |4.5|) . Being models, they are not uniquely defined 
by theory. Rather, they depend to some extent on ad hoc assumptions and param- 
eters which are being adapted to observation. At least the hadronization models 
are universal and can be used for different kinds of reactions. Probably the most 
sophisticated model is based on a colour string (a flux tube) that connects coloured 
partons. When stretched by the separating partons, the string breaks, pulling new 
qq pairs from the vacuum. When no more energy is left in the string, colour neutral 
hadrons are formed. 
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Chapter 3 

Inclusive Cross Sections 



3.1 The Structure Function Fo 



ZEUS H, H, H, ||] and HI 



1 



41, 421 have measured F 2 in a com- 



pletely new kinematic domain compared with fixed target experiments, most notably 
towards much smaller values of x, and towards much larger Q 2 . For kinematic rea- 
sons the cross section falls oc (1/Q) 4 (see eq. |2.1| ). Here we shall concentrate on 
the region of Q 2 less than a few thousand GeV 2 , because very little (though not 
uninteresting |43], ^4], |45|) data exist for higher Q 2 , especially as the hadronic final 
state is concerned. 

In Fig. p.l| F 2 is shown as a function of Q 2 for fixed x values. At x ~ 0.1 scaling 
is observed - F 2 does not depend on Q 2 . At x > 0.1 F 2 decreases with Q 2 due 
to parton splittings, the products of which are found then at smaller x. Therefore 
F 2 increases with Q 2 for x < 0.1 in accord with the qualitative discussion of section 
27T| . When plotted as a function of x for fixed Q 2 (fig. |3.2| ), F 2 exhibits a steep rise 
towards small x , which flattens at smaller Q 2 . The increase signals growing parton 
densities with decreasing x. 

From a DGLAP evolution of pre-HERA data this sharp rise could not be pre- 
dicted a priori, because input distributions at small x were not available. (Assuming 
however valence-like parton distributions at a small scale Qo ~ 0-3 GeV 2 , a sharp 
rise at larger Q 2 was predicted from DGLAP evolution p9| .) It was known though 
that asymptotically for Q 2 — > oo the small x behaviour is given by the DLL formula 
eq. |2.17| . Are the HERA data still consistent with DGLAP evolution, or is there a 
need for other effects, for example BFKL, which one may expect at very small x? 
It turns out that the data with Q 2 > 1 GeV 2 can be fit perfectly well with parton 
densities which obey the next-to-leading-order (NLO) DGLAP evolution equations 
(see Fig. |J) |4, 
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Standard QCD evolution appears to work over many 
orders of magnitude in both x and Q 2 \ 

When fitting F 2 ex (l/x) x , HI finds that the exponent A increases from A ~ 0.15 
toA« 0.40 between Q 2 = 1 GeV 2 and Q 2 = 1000 GeV 2 gOJ. F 2 rises faster than 
expected from the soft Pomeron model (A ~ 0.08), but less fast than expected from 
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Figure 3.1: The structure functions ^(x, Q 2 ) as a function of Q 2 for different x bins. For better 
visibility, a function C(a;) = 0.6(i — 0.4), with i the x bin number (i = 1 for x — 0.32), is added to 
F 2 in the plot. Shown are the HERA @ |§ and fixed target data ||, |47|, ^|[ , as well as a NLO 
QCD fit pj. 
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Figure 3.2: The structure functions F2{x,Q 2 ) as a function of a; for different Q 2 bins (Q 2 values 
are given in GeV 2 ). Shown are data from HI (yj and NMC |47), together with a NLO DGLAP 
QCD fit to the data with Q 2 > 1.5 GeV 2 from HI, NMC and BCDMS @. The starting point 
for the evolution was Qf, = 1 GeV 2 , and a s (m 2 z ) = 0.118 was set. 

the LO BFKL equations (A = 0.5 in LO; A is expected to decrease in NLO). In 
fact, both the x and the Q 2 dependence of F 2 can be attributed predominantly to 
the DLL formula eq. |2.17j , which can be displayed nicely with a suitable variable 
transformation ("double asymptotic scaling" |50[). A unified BFKL and DGLAP 
description of the F 2 data using the unintegrated gluon distribution (eq. |2.19 ) is 
also possible, with significant contributions from the lnl/x resummation ||51|| . The 
structure function data are thus compatible with pure DGLAP evolution, but cannot 
exclude significant contributions from BFKL evolution pl|. The structure function 



data are probably too inclusive to resolve the question of non-DGLAP evolution. 
One has to resort to less inclusive measurements on the hadronic final state, which 
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will be discussed in chapter |. 

The F 2 data can also be described by evolving flat or valence-like input quark 
and gluon distributions from a very low scale, Qq = 0.35 GeV 2 up in Q 2 as was 
shown by Gliick, Reya and Vogt (GRV) fJ9fl . The steep rise with decreasing x is 
achieved by the long evolution length from Ql to Q 2 (see eq. |2.17| ). The success 



of the GRV prediction came as a surprise for many, as perturbative QCD should 
only be applicable for Q not too close to Aq CD w 0.2 GeV , because otherwise 
a s (Q 2 ) = 12tt/[(33 - 2n f ) \n(Q 2 / A 2 QCD )] (LO) diverges. 

To conclude this section, it is quite satisfying that the evolution of the F 2 data 
can be described with parton densities following standard QCD evolution. In the 
next section the results of the NLO DGLAP QCD analyses will be given. However, 
the goal remains to calculate the measured magnitude of the growth (l/x) A from 
QCD, rather than tuning it with the starting point Q^ of the DLL evolution. 

3.2 QCD analysis of F 2 



From the DGLAP equations eq. [2.13| it is clear that the scaling violations of 



F 2 depend on both a s and the gluon density. In fact, for x < 0.01 and in low- 



est order one can derive the approximate formula |52 



dF 2 (x/2,Q 2 ) 10 a s (Q 



2\ 



dlnQ 2 27 vr 



xg(x,Q ), (3.1) 



because at small x the proton is dominated by gluons, and the scaling violations arise 
from quark pair creation from gluons. The full NLO QCD analyses now employed 
at HERA are of course more involved. Fig. [3.3| a shows the gluon density x ■ g(x, Q 2 ) 
extracted from NLO QCD fits to the F 2 data [53|. Previous data from NMC cover 
x > 0.01, and the HERA data extend down to x — 0.0001. In that region the gluon 
density increases sharply towards small x. 

Clearly, the density cannot increase forever; eventually saturation effects will 
set in |f)8| , |59[| . The effective transverse size of a gluon that is probed by a photon 
of virtuality Q is ~ 1/Q- Fluctuations below that scale happen, but cannot be 
resolved. When such gluons of transverse size ~ 1/Q fill up the whole transverse 
area offered by the proton, they will start to overlap and recombine. This would 
be a novel and very interesting situation indeed: high parton density, but Q 2 large 
enough for a s (Q 2 ) to be small! Given the size of the proton rj 1 fm, the critical 



condition for saturation effects to turn on can be estimated |59| as Xcritfi^crit, Q ) ~ 
Y7pj2 ~ 25 Q 2 ■ This value is by far not reached by the measured gluon density 
(Fig. |3.3| a). It could be however that saturation does not set in uniformly over 
the proton's transverse area, but starts locally in so-called hot spots |K|. In this 
case x cr it would be larger. The inclusive structure function data however do not 
require any saturation correction. Even a conspiracy of two new effects, a sharp rise 
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Figure 3.3: The gluon density xg(x,Q 2 ) at Q 2 = 20 GeV 2 f|, as extracted from NLO QCD 
analyses of the HI (4l| and ZEUS @ F 2 data together with the NMC result @. Shown are the 
±1ct error bands of the analyses. The values for the charm mass and for a s (m 2 z ) which were used 
in the analyses are given. The extracted gluon density is compared to the global analyses MRSR1 
H, CTEQ4m M and GRV94-HO jM. 



of F 2 from BFKL evolution dampened by saturation effects, cannot be ruled out. 
Probably saturation effects will be first seen in hadronic final state data 



From an analysis of the scaling violations dF 2 /dlnQ 2 



a. 



or equivalently, 



from a QCD fit to the F 2 data, the strong coupling constant can in principle be 
determined. From an analysis of the 1993 data a s (m 2 z ) = 0.120 ± 0.005(exp.) ± 
0.009(theor.) was obtained J6l| . It is estimated |62| that ultimately, with an in- 
tegrated HERA luminosity of 500 pb _1 , a s can be extracted by a NLO QCD 
analysis from the HERA structure function data with an experimental error of 
Aa s = 0.001 — 0.002, and a theoretical uncertainty of ~ 0.006. The theoretical 
errors could be reduced by higher than NLO calculations. 
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3.3 The Longitudinal Structure Function Fl 

One can express the longitudinal structure function Fl in terms of the cross section 
for the absorption of longitudinally polarized photons (see section |2.1|) , 

Q 2 (l-x) Q 2 . . 

F ^ = 1 2 ■ ^ W 7T'^' ( 3 - 2 ) 

Longitudinal photons have helicity and can exist only virtually. For transverse 
photons with helicity ±1 the spin is parallel to the direction of propagation, and the 
field vector perpendicular to it. In the QPM, helicity conservation at the electro- 
magnetic vertex yields the Callan-Gross relation (Fl = 0) for scattering on quarks 
with spin 1/2. This does not hold when the quarks acquire transverse momenta 



from QCD radiation. Instead, QCD yields the Altarelli-Martinelli equation |63 



i-X a O / \Al& 



F L (X,<?) = -t-X* -r 



4-7T Jx Z 



■:.\ 



^-F 2 (z,Q 2 ) + 8^1(1 --)-zg(z,Q 2 ) 

o Z 



(3.3) 



exposing the dependence of Fl on the strong coupling and the gluon density. At 
small x, the second term with the gluon density is the dominant one. In fact 
Fl(x, Q 2 ) = 0.3 • [Aa s /(Sir)] ■ xg(2.5x, Q 2 ) is not a bad approximation for x < 10 -3 



The extraction of F 2 from the cross section measurement (eq. |2.1|) so far had 
to make an assumption for Fl, because there existed no direct Fl measurements in 
the HERA regime. At large y, y w 0.7, this is a 10% correction. The argument can 
be turned around, and Fl can be extracted at large y from a measurement of the 
cross section, assuming that F 2 follows a QCD evolution and can be extrapolated 
from measurements at smaller y. This procedure has been carried out by HI |^2|, |41fl 
(Fig. p.4| ). The extracted Fl excludes the extreme possibilities Fl =F 2 and Fl =0 
and implies R = <JlJo t = F L /(F 2 — F L ) ~ 0.5, since F 2 m 1.5. This is self consistent 
with the gluon density extracted from the HI QCD fit. It should be noted however 
that the Fl data points at high y are somewhat above the expectation from the 
QCD fit. More precise measurements are necessary to clarify the situation. A 
measurement of Fl without theoretical assumptions on the evolution of F 2 will be 
possible by measuring the DIS cross section at HERA at two different centre of mass 
energies | 65fl , because these data will allow to vary y in eq. 2A while keeping x and 



Q 2 fixed. 

3.4 ct^o/ and the Transition between DIS and Pho- 
toproduction 

The total j*p cross section is shown as a function of W 2 in fig. |3.5| . The cross section 
increases with W 2 , and the slope increases with Q 2 . 
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Figure 3.4: The longitudinal structure functions Fl extracted from a QCD analysis of high y 
inclusive scattering data |42|, [Il[ at two different y values. The shaded bands give the range of 
expectations for Fl from a QCD analysis of structure function data with y < 0.35. The full line 
represents the solution Fl =F 2 . 



The data at Q 2 = are well described by the soft Pomeron parametrization 
by Donnachie-Landshoff (DL) PC|| , o~7ot P oc W 2 ' om . With increasing Q 2 , the cross 
section grows faster than the DL prediction. There appears to be a smooth but 
rather fast transition between the soft and the perturbative regions. Already at 
Q 2 ps 1 GeV 2 , a distinctively fast rise of the total cross section a^ oc W 2 ' x with 
A > 0.08 is observed. The HI fit result F 2 oc (l/x) A translates into o"7 ot p oc W 2X 
with A ps 0.2 - 0.4 for Q 2 > 1 GeV 2 . 

The perturbative GRV mechanism [|j| utilizing NLO DGLAP evolution provides 
an adequate description of the data down to Q 2 ~ 1 GeV 2 . It has not been 
possible to apply NLO DGLAP at lower Q 2 to describe the data. For example, 
below Q 2 ps 0.4 GeV 2 the GRV curves would turn over at small x, reflecting the 
valence-like behaviour of the GRV input partons at a low scale, by far failing to 



account for the data 35 . 



34 



■i 

2 10 

V 



10 



10 



10 



10 



• ZEUSBPC95 




* ZEUS 94 ° E6 
o HI 94 * ZI 
▼ H1SVTX95 * lo\ 


65 

:us, HI 7p 

v Wyp 

Q 2 (ZEUS) 
0.00 (x256) 
0.11 (x 128) 

0.15 (x 64) 

0.20 (x 32) 

0.25 (x 16) 

0.30 (x 8) 
0.40 (x 4) 

0.50 (x 2) 
0.65 (x 1) 
1.50 (x 1/2) 

3.00 (x 1/4) 
6.50 (x 1/8) 

1 


GKV(94) 

T-VT 


JJL/ 


H*M*<$^v* »fr --.... .^ *.».— ... 








♦ ♦ 


- 






♦ ♦♦♦ 


- 






♦ ,.♦♦♦ 








♦ •♦♦♦♦ .... 


0.23 (xl6) 
0.31 (x 8) ■-•, 
0.43 (x 4) o 


...(/ 








+ . ♦ ♦ ♦ ♦ 








, ..._• L 


'"6 


















♦ _.---- 


0.59 (x 1) -•0"°- 






.--#• 


-,-n- f '^"" 

~r~" T 


2.80 (xl/4) 0-^- 






, .*-*+ v * 


,^-'- (t 


7.16 (xl/8) --f 4 ^ * i? 

1 


, 







10 



10 



10 



10' 



W 2 (GeV 2 ) 



Figure 3.5: The total -f*p cross section as a function of W 2 for different values of Q 2 [J35[ . The 
data for virtual photon-proton scattering |3J| ^, 34, |5[ Q and photoproduction [[3(| [57], Q are 
compared to the predictions by the soft Pomeron model of DL |69| and the perturbative QCD 
model by GRV M. 
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Chapter 4 

Models for Hadron Production in 
DIS 

4.1 Reference Frames 

In this section the different reference frames and the variables used to describe the 
hadronic final state are introduced. Rough expectations for the event properties are 
derived mainly from phase space arguments. They will be refined by perturbative 
QCD and specific hadronization models. 

The laboratory frame 

The HERA experiments are performed in the "laboratory frame", in which the 
detector is at rest, and the incoming electron and proton beams are collinear. The 
z axis is defined by the proton beam direction. After the collision, the transverse 
momentum of the scattered electron is balanced by the one of the hadronic system. 
The hadronic final state is however better studied in a frame where the transverse 
boost is removed, such that the virtual photon and the incoming proton are collinear. 
Photon and proton 4- vectors are denoted by q = {E~* ,q) and P = (E p ,P). With 
E p 3> m p the proton mass can be neglected. 

The hadronic CMS 

The hadronic centre of mass system (CMS) is defined by the condition P + q — 0. 
The invariant mass of the hadronic system is given by W 2 = (P + q) 2 . The positive 
z-axis is usually defined by the direction of the virtual boson, q. All hadronic final 
state particles with 4-momentum p = (E,p x ,p y ,p z ) which have p z > are said 
to belong to the current hemisphere, and all particles with p z < are assigned to 
the target or proton remnant hemisphere. Target and current systems are back 
to back and carry energy W/2 each. They are not necessarily collinear with the 
incoming proton. In the QPM however, a quark in the proton with 4-momentum 
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xP absorbs the virtual photon and is scattered back with 4-momentum xP + q. 
The proton remnant retains the 4-momentum (1 — x)P. With these assumptions, 
the scattered quark and the proton remnant are both collinear with the incoming 
virtual photon and the proton (fig. |4.1| ). This is no longer the case if one considers 
intrinsic transverse momenta of the partons in the proton, or generates transverse 
momenta perturbatively by radiation. 



hadronic CMS 




' v. proton 




y 


W 




^ — 


x remnant quark 


2 



Figure 4.1: The hadronic centre of mass system (CMS). 



The Breit frame 

The Breit frame (BF) is defined by the conditions that proton and virtual photon 
are collinear, and that the virtual photon does not transfer energy, just momentum. 
Since q 2 = —Q 2 , q = (0, 0, 0, Q) with the same orientation as in the CMS, and 
P = (E p , 0, 0, —E p ). In the QPM, the incoming quark absorbing the virtual photon 
does not change its energy, but reverses its longitudinal momentum with magnitude 
Q/2 (fig. |4.2|) . The Breit system is also called the brick wall system, because in this 
picture the quark bounces back like a tennis ball off a brick wall. Neglecting the 
quark mass and transverse momenta, its 4-momentum is xP = {Q/2, 0, 0, —Q/2) 
before and {Q/2, 0, 0, Q/2) after the scattering. The incoming 3-momentum vector 
xP is merely reverted. 

Again, all particles with p z > are assigned to the Breit current hemisphere. 
The CMS and BF are connected via a longitudinal boost. The division into current 
and target hemisphere is frame dependent. In both systems, the remnant diquark 
is a mere spectator in the sense that its momentum remains unchanged. 
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Figure 4.2: The Breit frame. 
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4.2 Kinematics and Variables 

Let us consider a particle with 4-momentum p = (E,p) = (E,p x ,p y ,p z ) and mass 
m, which travels under an angle 9 with respect to the z axis. We define the trans- 



verse momentum pp = Jp 2 x + p 2 , the transverse mass rrip = Jm 2 + p\ and the 
transverse energy E T = E ■ p T /\p\. 



The rapidity j70| of a particle is defined as 



y :— - In — = In = coth —|. (4.1) 

y 2 E-p z m T E K J 

The rapidity transforms under a boost in the z direction with velocity (3 as 

y-*y , = y-coth/3. (4.2) 

Hence the shape of the rapidity distribution is invariant against longitudinal boosts. 
A useful relation is 

dy/dp z = lj E. (4.3) 

Often the mass of a measured particle is not known. One therefore defines the 
pseudorapidity 

n 

i] := — lntan-. (4.4) 

For E ^> in, r) « y . 

Energy and momentum of a particle are conveniently scaled by their kinemati- 
cally allowed maximum, leading to the definitions 

x F :=p z /pf^ x p :=\p\/p max x E :=E/E max . (4.5) 

For E 3> m and p z ^> px, these definitions converge: xp ~ x p ~x^ . We shall use 
z as a generic symbol for the variables defined in eq. |4.5| . 

The scaled longitudinal momentum xf is known as Feynman-x. In the CMS, 
xp ~ 2-p z /W, with — 1 < xf < 1. Similarly, in the Breit frame current hemisphere 



x 



v 



2\p\/Q with < x p < 1. None of the variables defined in eq. [4.5| is Lorentz 



invariant. A Lorentz invariant that is currently not being used at HERA is Zh := 
(Pp)/(Pq). In the proton rest frame it gives the fraction of the energy transfer taken 
by the produced hadron, Zh = Ejv. For large W and xp ,Zh > 0.1 the differences 
between z and xp are small. 

The rapidity can take any values with y min < y < y max (see fig. fO|), with 



?/max = In — 2/min = In — - — . (4.6) 

m m 

In the CMS, -y min = y max w In VT/m. At HERA with W < 300 GeV, pions can be 
produced with —7.7 < y < 7.7. Neglecting transverse momenta, the CMS rapidity 
y is related to xp (for xp > 0) by 

2/max - 2/ ~ ln(l/x F )• (4.7) 
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The difference of rapidities measured in the Breit and the CM systems is 

ln(Q/m) - \n(W/m) = \n(Q/W). (4.8) 
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Figure 4.3: The rapidity plateau. Indicated are the sizes of the target and current regions in the 
Breit and CM systems. The dimensions are chosen for typical HERA events with W = 120 GeV , 
Q 2 = 30 GeV 2 and for m — m,. 



4.3 Simple Mechanisms for Hadron Production 

Much can be learnt about the gross event features already from phase space consider- 
ations and a few simple assumptions. We start with the simple model of independent 
fragmentation, applied to the scatterd quark and the remnant. In the following sec- 
tions, this naive picture will be refined with perturbative QCD radiation and more 
sophisticated hadronization models. 

The Lorentz invariant cross section for hadron production may be written as 



E 



dV 



d 3 



a 



dp d(f>dyp T dp T ' 



(4.9) 



since dy/dp z = 
gration yields 



1/E. For a distribution that is isotropic in azimuth 0, the inte- 



^ d 3 a 



dV 



d'p dydpj^ 

We describe the inclusive densities to find a hadron h with energy fraction 
z := E/E q from the fragmentation of a quark q with fragmentation functions D^(z). 
(Dq(z) is not a probability density, because its integral is not 1.) They scale ap- 
proximately, that is they depend only on the fractional hadron energy z and not on 
the quark energy E q fflfl . The resulting hadron spectrum from quark fragmentation 
in the reaction ep — > e'qX is then 



1 dnh 



1 da E q e 2 Mx)D%(z 



N dz otot dz 



l^iq e qJq\ X ) 



(4.11) 
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Here cr to t is the total event cross section, N the number of produced events, and rih 
the number of produced hadrons of type h. f q (x) is the quark density function for 
flavour q in the proton. 

In the simple model of independent quark fragmentation [72|, a quark with energy 



E q fragments into a hadron with energy E according to a distribution function f(z), 
where z = E/E q , see fig. [4.4| . The process is iterated with a quark carrying the 
remaining energy (1 — z)E q until there is no energy left to produce a hadron. It 
is assumed that in the fragmentation process only limited transverse momenta are 
produced, usually parameterized with a falling exponential distribution in p\ with 
(p£) « 0.44 GeV 2 . 
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Figure 4.4: Independent fragmentation a la Feynman-Field W% . A u quark with initital energy 
E q branches into a 7r + meson with energy E = zE q , and a u quark with energy (1 — z)E q , and so 



The model describes the gross features of fragmentation |72| , namely energy 
independence, and for small z a behaviour 

D h q (z) oc l/z. (4.12) 

Since dy ~ zdz, this gives a uniform number density of hadrons in rapidity, the 
so-called central (in the CMS) rapidity plateau: 



N dy N dz q 



■ z oc D q (z) ■ z = const. (4.13) 



Fig. |4.5| a gives the relation between the longitudinal momentum p z = z ■ p™ ax ~ 
z ■ W/2 for a typical HERA situation. \xp | < 0.1 translates roughly to the rapidity 
region \y\ < 3 (fig. |4.5|a). 

A convenient parametrization for the fragmentation functions is 

D h q (z)=a- z (l-z) b . (4.14) 

Energy-momentum conservation requires 

J2[ 1 zD q \z)dz = l, (4.15) 

h Jo 
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Figure 4.5: a) Relation between rapidity y and longitudinal momentum p z for a particle with 
m = m„ and px = 0.6 GeV at W = 120 GeV , typical for HERA. The curves for pseudorapidity 
•q and rapidity y are indistinguishable. Differences between y and r\ become noticeable for smaller 
Pt ; for pt — 0.1 GeV it is at most 0.3 rapidity units. j/ max is given by p™ ax rs W/2. b) The 
fragmentation function D(z) = — (1 — z) 2 ' 5 . 



where the sum runs over all hadron species /i. If we do not distinguish between 
hadron species, we can define D q (z) := J2h Dg(z). The normalization a is then fixed 
by 

/ zDJz)dz = a/(b + 1) = 1. (4.16) 

Jo 

A useful parametrization for light quarks q = u,d, s, accurate to 20%, is 



D q (z) 



3.5-(l-z 

z 



,2.5 



(4.17) 



see fig. [4.5| b. The average number of hadrons per unit rapidity in the plateau region 
(z small) is then 

of which roughly 2/3 will be charged. The rapidity plateau is depicted in fig. [O . 
There will be a kinematic fall-off for y — > y max , y m >m ■ The available rapidity range 
|/max — 2/min increases logarithmically with the quark energy (longitudinal phase 
space), and so does the average total hadron multiplicity. 

The average total hadron multiplicity in a quark jet is given by 



(4.19) 



(n) = / D q (z)dz, 
where z m [ n = m/E q is determined by a typical hadron mass m. From 



(4.20) 
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we obtain a logarithmic scaling law for the multiplicity 

(n) = aln^ + c, (4.21) 

m 

where c is the integration constant. 

Though the independent fragmentation picture describes crudely the data on 
e.g. e + e~ — > qq — > hadrons, it has its limitations. It cannot be applied consistently 
for small z, z < 0.1, it is not Lorentz invariant (the constants a and c are frame 
dependent), and energy momentum conservation has to be enforced by hand after 
both quarks have fragmented. The lack of Lorentz invariance even leads to contra- 
dictions. Let W be the total CM energy. In the CMS, one calculates for the total 
multiplicity (n) = 2 ■ (c + a In — ) with E q = W/2. In a frame where one quark takes 
almost all of the energy, and the other one only very little, one gets a different value 
( n ) = c + a In ^ w ith E a = W 

The discussion so far covered only the lowest order processes, where hadron 
production is determined by the assumption of limited transverse momenta and 
longitudinal phase space as a simple hadronization model. Though simple, it pro- 
vides a good guideline for hadron production. QCD radiation modifys this simple 
picture and will be discussed next (section |4.4j) . Afterwards improved hadronization 
models are introduced that do not suffer from the problems with the independent 
fragmentation model (section |4.5|). 



4.4 Perturbative QCD Radiation 

The perturbative part of QCD radiation is treated either with fixed order matrix 
elements, or with parton showers, or a combination of both. It depends on the 
application which treatment is most appropriate. In general, the matrix element will 
become more important with increasing hardness of the interaction. Perturbative 
QCD makes predictions for partonic final states, but observed are hadronic final 
states. To make contact with the experimentally accessible world, hadronization has 
to be taken into account. Models for the non-perturbative effects of hadronization 



are discussed in section [4.5|. Excellent reviews on perturbative QCD evolution and 



on hadronization are [74, 75]. Further information specific to DIS can be found in 
|76| . We start the discussion of perturbative QCD with matrix elements and parton 
showers. An alternative technique which is not based on Feynman diagrams is the 
dipole radiation approach, to be discussed afterwards. 

Matrix elements 

In fixed order perturbation theory, the incoming parton flux is folded with the 
matrix element [[FT]] for electron-parton scattering which leads to some final state 
parton configuration. In LO possible final state configurations for eq scattering are 
(apart from the scattered electron) q (QPM), qg (QCDC), and for electron gluon 

42 



scattering (BGF) qq (fig. [4.6[). In NLO, one additional parton can be emitted, and 



so on. 




QPM 




BGF 



QCDC 



Figure 4.6: Matrix elements in DIS. In leading order, scattering on a parton i can proceed via 
the graph that corresponds to the simple quark parton model (QPM), or the QCD analogue of 
Compton scattering (QCDC), or the boson-gluon fusion graph (BGF) in case the incoming parton 
i is a gluon. The graphs on the right hand are the LO realizations of the "blob" in the left graph. 

The calculations are ususally done numerically. For state of the art NLO cal- 
culations the programs DISENT @, MEPJET J7| and DISASTER++ @ are 



available. Originally such programs could only be used for jet analysis. The present 
programs are more flexible and allow also calculations of for example event shape 
variables. 



Leading Log parton showers 

A complete calculation to an order higher than NLO appears prohibitive. For many 
applications, where higher orders are important, the parton shower ansatz is used. 
For example, higher orders can be summed to all orders in the leading log approx- 
imation. In the DGLAP (leading logQ 2 ) approximation, they result in splitting 
functions which give parton emission probabilities. 

When a parton with 4-momentum p radiates, it changes its virtuality T 2 := — p 2 . 
The evolution parameter t := ln(T 2 /A 2 ) increases in the initial state parton shower 
(space-like shower) by successive emissions, and decreases in the final state parton 
shower (time-like shower). The probability V that a branching a — > be will take 



place during a small change dt is given by the evolution equation [15 



dV a . 



•bc 



i\t 







1 dz ^flp^ bc{z) . 

Z71 



(4.22) 



The functions P a -+bc(z) are just the Altarelli-Parisi splitting functions in eq. |2.14j , 
p p 

1 a^rbc J ba- 
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Such calculations are usually done with Monte Carlo event generators, where 
the parton shower evolution is simulated step by step according to the emission 
probabilites, until the whole event is generated. Perturbative evolution is stopped 
at some small scale T 2 of 0(lGeV 2 ), when it becomes unsafe to apply perturbation 
theory due to the growth of a s . 

When combined with the exact fixed order matrix element to take care of very 
hard emissions that are not properly covered with leading logs, the Monte Carlo 
event generators can be expected to provide a good representation of what is actually 
happening in DIS eventsQ (fig. |4.7|) . Ambiguities result from the way the matrix 
element and parton shower are combined ("matching"), coherence effects are treated, 
divergencies of the matrix element are cut-off, and from other approximations in the 
implementation. Unfortunately todays generators implement the matrix element 
and parton showers only to LO. 

Examples of DIS Monte Carlo generators that are based on leading log DGLAP 
parton showers are LEPTO [§]], HERWIG g^ and RAPGAP g£§. They are ex- 



pected to be valid where the DGLAP approximation is valid. As DGLAP based 
models, the initial state parton shower is strongly ordered in kx , increasing from 
the proton towards the matrix element. 

The leading log approximation (LLA) can also be applied directly for specific 
observables, without utilizing an event generator. Most relevant is the modified 
leading log approximation (MLLA) which takes into account destructive interference 
between soft gluons, in conjunction with the assumption that the resulting parton 
spectra resemble the observable hadron spectra, up to a constant factor (local parton 
hadron duality, LPHD). 

The colour dipole model 

Another type of parton shower model is the colour dipole model (CDM) for QCD ra- 



diation |84[. The colour charges of the scattered quark and the remnant are assumed 
to form a colour dipole, from which gluons can be radiated (fig. |4.8| a). Subsequent 
gluon radiation emanates from dipoles spanned between the newly created colour 
charges and the others, and so on. To good approximation it can be assumed that 
these dipoles radiate independently. The CDM uses the LO cross section for the 
emission of a gluon with transverse momentum p? at rapidity y in the soft gluon 



approximation |85| 



da = '-^^fdy. (4.23) 



n c a s dpf. 

2"7T p\ 

The cross section is uniform in rapidity and \wp\. Kinematically, the phase space 
is bounded by \y\ < \ii(W/pt), where W is the total energy of the radiating system. 
Gluons with pt above a certain cut-off, pt > A, lie inside the triangle in fig. |4.8|b. 



1 Remembering the double slit experiment, we should feel a bit uneasy about a statement of 
what is actually happening. 



44 




Figure 4.7: Elements of an ep event generator. As an example an event is shown with the BGF 
matrix element (ME). The initial and final state parton showers (PS) produce additional emissions. 
Finally, the partonic final state (the parton level) is hadronized to yield the observable hadrons 
(the hadron level). 



quark 




K=ln(pj.A 2 )i 

L 




remnant 



Figure 4.8: a) CDM in DIS: gluons are emitted from successively built colour dipoles. b) The 
phase space for the emission of gluons with pr > A in the variables k = ln(p^/A 2 ) and y. The 
phase space is bounded by the triangle with height L = ln(W 2 /X 2 ). Subsequent emissions have 
decreasing pt , pn+i < Pn, but no ordering in y. If the emissions are reordered in y, they lose 
the pt ordering. 



In DIS one of the colour charges is not pointlike, the proton remnant. The 
suppression of radiation with short wavelength from an extended source is taken 
into account with a special parameter of 0(1 fm) that describes the extendedness of 
the remnant remnant^]. The CDM is implemented in the program ARIADNE [36 



2 In the recent version also the colour charge of the scattered quark is spread out, depending on 
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In ARIADNE the QCDC graphs are covered by dipole radiation (and corrected to 
yield the exact LO matrix element contribution), but for the BGF graph the matrix 
element is used. Further radiation is then according to the CDM. 

In contrast to the LL parton shower in the CDM it is not possible to distinguish 
between initial and final state radiation, or to reconstruct an "evolution path" for 
the parton in the proton that is hit by the virtual photon. Kinematically, the 
kr phase space for radiation is bounded by the k? of the previous emission. In 
rapidity the emission probability is uniform. Therefore the final gluon configuration 
is not ordered in kr, if one arranges them according to increasing rapidity (see fig. 
fL8|b). Rapidity is directly related to lnxj for fixed W and kn (see section |8jJ) . In 
that respect the CDM is similar to what can be expected from BFKL evolution |[7 , 



though perhaps in a somewhat unconventional fashion [p8| . Whereas the LL parton 
showers are connected with the intuitive picture of an evolution path, it is probably 
fair to say that the significance of the dipole model is not yet understood very well. 
We shall see that this model provides effortlessly overall the best description of most 
final state data. This is probably the main reason why it is being discussed, whether 
or not it has a deeper physical justification not yet appreciated by the community. 

The Linked Dipole Chain Model - LDC 

The CCFM approach for the hadronic final state has been reformulated in a picture 
similar to the CDM, giving rise to the Linked Dipole Chain model (LDC) p4j . This 



model is particularly interesting, as it should converge to the DGLAP and BFKL 
predictions in their respective regions of validity, and handle properly cancellations 
between real and virtual emissions for the final state. In the LDC also the case is 
treated where the largest pr of the process is not attached to the virtual photon 
vertex. Some promising results from a Monte Carlo implementation of the LDC 



have been obtained [39 , but the program is not yet publicly available. 



Coherence 

Quantum mechanical interference is an important issue for the hadronic final state. 
With matrix elements that sum and square amplitudes interference effects are taken 
into account naturally. Coherence effects in parton showers that rely on emission 
probabilities, not amplitudes, require special attention. One may distinguish a) 
interference between initial and final state radiation; and b) interference between 
partons emitted either in the initial state or in the final state parton shower. 

In the CDM there is no distinction between initial and final state radiation, so in 
the dipole approximation problem a) does not arise. Also soft gluon interference is 
automatically taken into account by the dipole mechanism. In other generators such 
interference effects can be realized by phase space restrictions for parton emission. 
For example, the effect of destructive interference between subsequent emissions in 



the "size" (ex 1/Q) of the final state quark that can be resolved with the photon of virtuality Q 2 
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a parton shower can be approximated by imposing a posteriori "angular ordering" : 
only emissions are allowed with a smaller opening angle than the previous one, 
9 i+ i < 9i. 

The physical reason for this condition is that large wavelength quanta cannot 
be emitted from dipole sources with small transverse dimensions. This can be seen 
easily for e + e~ cascades pD |. Consider the example fig. [4.9| a, where a virtual photon 
branches into an e + e~ pair, and the e~ radiates a photon. The formation time of the 
final eey state can be estimated as the lifetime of the intermediate e~ from its off- 
shellness M, r = 7 • r = E'/M ■ 1/M. We have M 2 = (p + kf « 2EE 1 {1 - cos6 2 ), 
and 1 — cos 9 2 ~ $1/2 for small angles, Since for soft photons, E 7 <C E' « E, 
t ph 1/(E^6 2 ). During this time the e + e~ pair have separated a transverse distance 
Sr fh rd\. The transverse wavelength of the emitted photon is Ay = 1/kx with 
k T = E 1 9 2 . If \t is larger than Sr, the photon cannot resolve the individual charge 
of the e~ , it rather sees the combined charge of the e + e~ pair, which is 0. For 
allowed emissions it follows Sr > At, or Q\ > 6 2 , angular ordering. 

In QCD the situation is a bit more complicated, because the incoming parton 
is colour charged (fig. |4.9| b) . Parton b cannot be radiated incoherently from parton 
v when 9 2 > Q\. In that case it would see the combined colour charge of a and v, 
which is the colour charge of the incoming parton i. Parton b with 9 2 > 9\ can 
therefore be treated as effectively being radiated from i before a branches off, thus 
restoring angular ordering. 






V, 



b) v 

Figure 4.9: a) Angular ordering in an electromagnetic cascade. The 4-momenta of the outgoing 
and intermediate electron are p — (E,p), p' = (E 1 ,p*), and of the outgoing photon k = (2? 7 , k). b) 
Angular ordering in a QCD cascade. 



4.5 Hadronization 

Experimentalists follow two different philosophies how to prepare their hadronic 
final state data. 

• In one approach, the data are corrected just for detector effects in order to 
remove any experimental bias and be able to compare the data to other ex- 
periments. The data are corrected "to the hadron level", they represent a 
measurement of observable hadrons. The data can also be compared to QCD 
predictions that include hadronization effects, but not directly to perturbative 
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QCD calculations (unless one is interested to study the difference between 
observable hadrons and partonic calculations). 

• In the other approach, a hadronization model is used to correct the data for 
hadronization effects back "to the parton level" . The data can then be com- 
pared directly to perturbative QCD calculations. A problem of principle arises: 
partons are not observable, so strictly speaking an observable "at the parton 
level" does not exist. With great care one may proceed nevertheless. One has 
to specify the parton level in perturbative QCD (like NLO, LLA, including 
cut-offs etc.), and correct the data for all other effects beyond that: higher 
order corrections plus hadronization effects. The observable thus becomes 
scheme dependent. There is another price to pay: the data have picked up 
a theoretical bias due to the assumptions on higher order and hadronization 
effects. In practice this bias is studied by comparing different Monte Carlo 
generators. The encountered differences are treated as systematic error of the 
measurement. 

We summarize the main approaches to hadronization, some of which have already 
been mentioned. It has to be emphasized that the boundary between perturbative 
QCD evolution and hadronization is not uniquely defined. Ideally, predictions for 
final state observables would not depend on how the boundary is defined, for example 
on the lower cut-off for virtualities in the parton shower. 

We start with hadronization models that are implemented in Monte Carlo gen- 
erators (independent fragmentation, string and cluster model). Other approaches 
that model only specific aspects of hadron production are discussed afterwards. 



Independent fragmentation |72[| 



From the fragmenting parton mesons are produced, carrying a certain fraction z 
of the original energy, see section [4.3| and fig. W!\ The branching is repeated with 



the remaining energy until all the energy is used up. In the physical picture for 
independent quark fragmentation the quark forms a meson with an antiquark from 
a qq pair, with the remaining quark continuing fragmentation. It is assumed that 
the distribution of z is energy independent, leading to scaling fragmentation func- 
tions. Furthermore, all partons of an event fragment independently. That leads to 
inconsistencies like non-conservation of energy-momentum, that have to be cured 
by hand. 

Transverse momentum components are assumed to be Gaussian distributed^ p2J . 
which results in an exponentially falling distribution in p\. From the parametriza- 

3 The transverse momentum pt which a primary hadron receives from the fragmentation process 
can be described by Gaussians in p x and p y with widths a, 



dp x dpy ct\/2tt \2a 2 J (j^/2tt \ 2er 
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tions used in fragmentation models p2| , |9T| , (p f ^ ag ) ~ (0.3 GeV ) 2 is obtained for 
primary hadrons. 

String fragmentation ||93|1 

In the Lund string model, a quark and an antiquark moving apart stretch a colour 
field between them. The field is thought to be string-like with constant energy den- 
sity per unit length of (9(1 GeV /fm), and with transverse dimension of (9(1 fm). 
When the stored energy becomes large enough, the string can break up and create 
a qq pair from the vacuum (fig. |4.10| a). The newly created q and q terminate the 
loose string ends, such that the new string pieces are themselves colour neutral. 
The process is iterated with the new strings until all the available energy is used 
up. Transverse momentum components result from a tunneling process; they are 
parametrized as in eq. |4.24| . The string model gives rise to scaling fragmentation 



functions, a rapidity plateau with uniform particle density, and a logarithmic in- 
crease of particle multiplicity (n) = c + a In W, where W is the string invariant 
mass. 

String fragmentation can be formulated in a covariant way, it does not suffer 
from the mentioned shortcomings of independent fragmentation. Furthermore, some 
quantum mechanical interference effects in gluon radiation are taken into account 
with string fragmentation, due to the way gluons are treated. Carrying two colour 
charges, gluons are always the endpoints of two strings. In a qqg configuration the 
gluon is realized as a kink in the colour connection between the qq pair. Destructive 
gluon interference leads to a suppression of radiation (and finally hadron production) 
in the region opposite to the gluon. Because in the string model there is no string 
spanned directly between the q and the q this "string effect" @j appears naturally 
in the Lund model. 



Cluster fragmentation [p5| 

According to the idea of "preconfinement" [p6|| , colour connected partons tend to 
be close in phase space towards the end of the perturbative phase. This property is 
exploited in the cluster fragmentation model. After the perturbative evolution all 

This leads to an exponentially falling distribution in p\ with (jpj) = 2<r 2 , 



dn 1 / —p\ 



2 - exp ( ) . (4.25) 



It follows, that the distribution in px is neither Gaussian nor exponential, 



- ■ ex P hr4 , (4-26) 



dpx o 2 \ 2a 2 

with (p T ) = y/Za = ^y/ijprj. The default parameter in the Monte Carlo program JETSET (9l) 



2 2 

is a = 0.36 GeV 
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Figure 4.10: a) String fragmentation, qq pairs are created when the stretched string between 
the separating qq pair breaks up. Mesons are formed when the string energy is too small to create 
further qq pairs, b) Cluster fragmentation. Gluons are represented by double lines, corresponding 
to the colour and anti-colour they carry. After the perturbative evolution gluons are split non- 
perturbatively into qq pairs, from which colour neutral clusters are formed. They decay into 
mesons. 



gluons are forced to split into qq pairs. Colour connected qq pairs are combined 
into low mass colour neutral clusters, which decay into mesons according to phase 
space (fig. |4.1U| ). Due to the small cluster mass, relatively few parameters are needed 
to describe the cluster decay. Nevertheless, large mass clusters may occur; they are 
split longitudinally in a fashion similar to the string model before they are allowed 
to decay. 



Local Parton Hadron Duality (LPHD) 

The hypothesis of local parton-hadron duality (LPHD) |)7] is made to allow pre 



dictions for hadron spectra from perturbatively calculated parton spectra. It states 
that any hadron cross section depending on a quantity ( is related to the corre- 
sponding partonic cross section simply by a constant factor, dcr^/d^ = c ■ da p /d( p . 
The idea behind this seemingly surprising ansatz is that if perturbative evolution is 
used for multi-parton emissions down to a very low scale, there is not much energy 
left for hadronization to change the shape of the distribution. This hypothesis could 
be applied successfully to a variety of phenomena, but one should keep its limita- 
tions in mind ||98|| . For example, when a quark and antiquark separate and do not 
radiate any gluons (unlikely, but possible), LPHD would predict a gap in rapidity in 
between, devoid of any hadrons. A more physical assumption is that the colour field 
stretched between quark and antiquark leads to the production of hadrons filling the 
gap. At HERA, the LPHD hypothesis has been applied to charged particle spectra, 



see sections K2\ and 6.1. 
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Fragmentation functions 



The number density to find a hadron h with energy fraction z := E/E q from the 
fragmentation of a quark q is described by a fragmentation function D^(z), and 
similarly for a gluon g. They scale approximately, that is they depend only on 
the fractional hadron energy and not on the quark energy |7I| . QCD effects like 
gluon radiation off the quark make the fragmentation functions scale dependent: 
Dq(z) — > Dy(z,t), where the scale t is provided by the hard process (fig. [4.11| ). 
While the fragmentation functions themselves are not calculable in perturbative 
QCD, their scale dependence can be calculated in the same way as the scaling 
violations of the structure functions in DIS 



3S 




Figure 4.11: Scaling violations of fragmentation functions In this case the scale is given by the 
virtuality of the photon, t = Q 2 . 



The fragmentation functions obey an evolution equation 



dD h/i (z,/j, 2 D ) 
<91n \i 2 D 



E 



a s {\x 



2n 



^Pji(u, a s (fM 2 D ))D h/j (^, fj, 2 D ). 



(4.27) 



t = fi 2 D is the factorization scale between the parton production process and the 
parton fragmentation process. The splitting functions Pji for branchings i — > j are 
in leading order the same as for the parton densities in DIS, 



eq. 



2.14, but differ in 



higher order [Q. They have been calculated up to next-to- leading order ||100|| . 
The fragmentation functions are universal, that is process independent. 



For 



example, quark fragmentation functions measured in e + e 
to predict the hadron distribution in the quark jet for ep - 



— > hadrons can be used 
e'qX. 



High energy and power corrections 

Due to the limited px generated in hadronization, it appears plausible that with 
increasing energy Q the relative importance of hadronization effects diminishes for 
suitably chosen observables. It is expected that for infrared safe quantities (see 
section [2.6|) like jet observables and some event shape variables their perturbative 
prediction approaches the value that is observable with hadrons. The difference, 
due to higher orders and hadronization, should decrease like an inverse power of 
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the energy oc (1/Q) a , a "power correction" |101 |. Power corrections have been 



investigated for thrust and other event shape variables at HERA, see section 6.2 



4.6 Rapidity Gaps 



102, 103 



At HERA there exists a class of events known as "rapidity gap events 
104 , |105|| . They have given rise to tremendous activity, which would deserve a review 
in its own right (recent summaries are ||). Here as much is covered as is relevant 
for our purpose. 

Rapidity gap events exhibit a sizeable region of rapidity without any hadronic ac- 



tivity - the rapidity gap, see fig. |4.12| . In normal events a colour field forms between 
the scattered quark and the proton remnant, leading to a homogeneous distribution 
of hadrons populating the available rapidity range. Large gap sizes between neigh- 
bouring hadrons from normal fragmentation are expected to be unlikely (though 
possible). Large gaps can be produced however in processes where two systems are 
produced that are not colour connected. 

ZEUS 1994 



W 



a) 
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Figure 4.12: a) Rapidity distribution for normal events (top) and for rapidity gap events 
(bottom). Mx and My are the invariant masses of the two systems X and Y separated by the 
largest gap in the event. The system Y is the proton or its remnant. The proton direction is to the 
right. r/ m ax is the forward rapidity edge of the system X. A large gap implies a small mass Mx- 
b) ?7max distribution from ZEUS (uncorrected) in the laboratory system. The data are compared 
to a normal DIS model (ARIADNE |^6| ) and to a model for rapidity gap events based on Pomeron 
exchange (POMPYT |l06| ). For most events 77 max ~ 3.5, given by the forward edge of the main 
calorimeter. Rapidity gap events are found in the tail at small rj mSiX . 



Events with large rapidity gaps had been anticipated for HERA from processes 
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though the large rate came as a surprise: about 10% 

a 



with Pomeron exchange |[107j 

of all DIS events have a large rapidity gapQ. In the Ingelman-Schlein model 
Pomeron is exchanged in the t channel between the proton and the virtual photon, 
see fig. |4.13| a. The proton may be left intact, or fragment into a low mass system 
My < 1.6 GeV . In the Ingelman-Schlein model, the virtual photon interacts with a 
parton in the Pomeron, thus probing the structure of the exchanged Pomeron. The 
resulting system with mass Mx is not colour connected with the proton remnant 
system, because the Pomeron carries no colour. A rapidity gap between the two 
systems is the consequence. 





a) 

Figure 4.13: a) DIS with Pomeron exchange. The proton momentum fraction carried by the 
Pomeron, xjp, can be reconstructed by xip = (M^ + Q 2 )x/Q 2 . The variable j3 := x/xjp = 
Q 2 /(Mv + Q 2 ) gives the Pomeron momentum fraction carried by the parton that is hit by the 
virtual photon, b) DIS with colour rotation. After the hard interaction (BGF, broken circle) a 
soft gluon is interchanged with the proton. 

Suppose the above picture to be valid, the partonic structure of the elusive 
Pomeron can be probed! In fact, the evaluation of the rate of rapidity gap events 
suggests that the exchanged object is predominantly gluonic, and that most of the 
Pomeron momentum is carried by one single gluon ||. The analyses suggest fur- 
thermore that the Pomeron intercept in DIS is somewhat larger than expected from 
the soft Pomeron, and possibly increasing with Q 2 |3[. 

Other models assume that rapidity gaps are due to changes in the colour struc- 
ture of the event after the hard interaction has taken place 



1091 , TO - 



At low x, 

most events are due to photon-gluon fusion, resulting in a colour octet - antioctet 
configuration in the final state. In the string model, there would be two fragment- 
ing strings spanned (fig. |4.14| a). It is assumed that while the produced quark and 
antiquark travel to leave the proton, soft gluons are exchanged between them and 



the proton remnant, see fig. |4.13| b. These gluons do not change the momenta sig- 
nificantly, but may rotate the colour configuration of the quark - antiquark pair. 
The final configuration may turn out to be a colour singlet, in which case there 



4 Often these events are called "diffractive" , or being due to a diffractive process. The meaning 
of these terms is not yet uniquely defined amongst the experiments and theorists. We therefore 
prefer the term "large rapidity gap" , because the existence of a gap of a certain size is well defined 
event by event. 
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would be no colour field connecting it with the remnant (fig. f4.14| c). The attractive 
feature of this model is an explanation of the absolute rate of rapidity gap events. 
If soft gluon exchange results in a random colour configuration, the ratio of colour 
singlet to colour octet configurations would be 1 : 8, and 1/9 ~ 10%! In a related 
semiclassical approach, partonic fluctuations of the virtual photon are scattered by 
the colour field of the proton [ |1 1 1 1 . 




Figure 4.14: String configurations for the normal case (a), and after soft colour interactions, 
leading to either longer strings (b) or a rapidity gap (c) 

Both models describe at least qualitatively the HERA data on rapidity gap 
events. Detailed studies will be necessary to discriminate them (for example in the 
final state of the system X [ |112|| ), or perhaps discover them to be not exclusive at 
al|]. For the discussion of the hadronic final state in DIS it will be important to keep 
the existence of rapidity gap events and possible production mechanisms for them 
in mind. Some analyses [|113| , |114| , |115| , |116| , |117| , |118| explicitly exclude rapidity 
gap events, others don't, and some study the difference between the two samples. 
The possible existence of a mechanism to change the colour configuration will have 
a bearing on the interpretation of certain final state data at small x (section |8.1| ). 
Colour reconnections are also being discussed in connection with WW fragmentation 
at LEP P|, [120 . 



4.7 Monte Carlo Generators 

We shortly summarize the DIS Monte Carlo generators. They incorporate the QCD 
evolution in different approximations and utilize phenomenological models for the 
non-perturbative hadronization phase. The events are generated according to the 
electroweak cross section d ^Q 2 with experimentally determined input parton densi- 



ties. QED radiative effects are treated with special programs [ |121| , |122 ], 

The Monte Carlo results depend on the model parameters, which can be adjusted 
to the data. Information on generator tuning can be found in ||123| , |124| , |125| , [76[ . 



Pomeron exchange where most of the momentum is carried by one gluon in the Pomeron is 
not so much different from boson-gluon fusion, where the hard dynamics is determined by the 
incoming gluon, and soft gluons are only needed for colour neutralization 
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In general, the generators now give an overall satisfactory description of the data. 
None of the generators however is able to describe all aspects of the data. 



LEPTO 

LEPTO |81] is based on the LO QCD matrix element with leading log DGLAP 



parton showers for soft emissions. Therefore the model is often called MEPS (Ma- 
trix Elements + Parton Showers). Angular ordering is imposed to model colour 
coherence. The rather sophisticated Lund string model as implemented in JET- 
SET |91|] is used for hadronization. A special feature are "soft colour interactions" 
(SCI) [|109| in the hadronization phase. They lead to string rearrangements. This 



feature is now default, because it describes at least qualitatively events with a large 
"rapidity gap" that are seen in the data. 

RAPGAP 



Originally, RAPGAP [Q was created for rapidity gap events. They are described 
by scattering on a colour neutral object, the Pomeron, which is emitted from the 
proton according to a flux factor, and whose internal structure is described by parton 
density functions. The present version includes also normal DIS event generation, 
very similar to LEPTO. In the latest version it is also possible to simulate a resolved 
(not point-like) component of the virtual photon. 

HERWIG 

HERWIG [[82] is based on LL parton showers, with additional corrections to describe 



properly the hard interaction. Final state radiation is angular ordered, initial state 
radiation is ordered in E ■ 6, where E and 9 are energy and angle of the radiated 
parton. The HERWIG philosophy is to model the perturbative phase as accurately 
as possible, and to hadronize with a relatively simple cluster fragmentation model 
Though there is no mechanism foreseen for rapidity gap events, it turns out 



that they are generated at a surprisingly large rate. Also with HERWIG a resolved 
component of the virtual photon can be simulated. 

ARIADNE 

In ARIADNE |SB| perturbative QCD radiation is modelled with radiating colour 
dipoles according to the CDM. The BGF graph is added by hand with its LO 
matrix element. Hadronization is performed with the string model as implemented 
in JETSET. ARIADNE allows to model rapidity gaps either by scattering on a 
Pomeron, or by colour reconnections, but these options are by default not activatedR. 



6 It was found that colour reconnections in combination with the colour dipole model do not 
give a good description of rapidity gap events at HERA |l26|. 
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Chapter 5 

General Event Properties 



5.1 Energy Flow 

The first analysis of the hadronic final state in DIS at HERA was based on 88 events 



recorded in 1992 ||127|| . The measured energy flows and particle spectra were in rough 
agreement with what was expected from some models including QCD radiation and 
fragmentation, and exluded certain other radiation scenarios. 

The flow of transverse energy as a function of pseudorapidity rj provides a 
very simple, global characterization of the hadronic final state. It is measured with 
the calorimeters and includes all produced particles, except the scattered electron. 
The E T flow as a function of pseudorapidity is defined as 

i^l = ±V [ ET ^^dE T (5 1) 

N dr, a^J drjdErp ' ^ 

where the sum extends over all particle species h with cross sections o^. The defi- 
nition is analogous for the E? flow as a function of azimuth <ft. The distribution is 
normalized to the total number of events N, corresponding to the total event cross 
section a. 

It is instructive to examine the flow of hadronic transverse energy as seen in 



the HERA detectors in the laboratory frame, fig. |5.1| In this frame one measures 



according to eq. |5.1| the energy flow transverse to the beam line as a function of 



pseudorapidity rj = — lntan(0/2) and azimuth (p. Here 9 is the angle with respect 
to the proton beam axis (laboratory +z), and is the angle between the energy 
deposition and the scattered electron in the plane transverse to the beam line. One 
sees the current jet from the scattered quark emerging opposite to the scattered 
electron in (f> (fig. ^T] right). The rapidity region from the current jet towards the 
remnant is filled with particles due to the colour field stretched between the scattered 
quark and the remnant (the main calorimeter acceptance ends at typically rj = 3.5), 
see fig. |5J] left. 

Significantly more transverse energy is observed between the current jet and the 
beam pipe than is expected from the QPM without QCD radiation, but including 



56 




Figure 5.1: Transverse energy flow in the laboratory frame for x < 1CP 3 . The plots are normalized 
to the number of observed events N. Left) Et as a function of pseudorapidity rj by HI [ 1 13 1 . 
The proton beam direction is to the right. Shown are statistical and systematic errors added 
in quadrature, except for a 6% scale error. The data are compared to the QPM without QCD 
radiation, and to models including QCD radiation: CDM (colour dipole model, ARIADNE 4.03), 
MEPS (matrix element plus parton showers, LEPTO 6.1), PSWQ (parton shower only, LEPTO 6.1, 
with maximal virtuality in the parton shower W ■ Q), and HERWIG (HERWIG 5.7). Right) Et as 
a function of azimuthal angle <fi with respect to the electron scattering angle in the transverse plane 
for clusters with 9 > 10° by ZEUS |l28fl . Comparisons are made to different model predictions: 
matrix element only (ME, dash-dotted), matrix element plus parton shower (MEPS, full line), 
parton showers only (PS) with either W 2 (dashed) of Q 2 (dotted) as maximum for the virtuality 
in the PS. The HI data are corrected for detector effects. The ZEUS data are uncorrected, and 
compared to the models including detector simulation. 



fragmentation. The QPM predicts an Et of roughly 1 GeV per unit pseudorapid- 
ity; one arrives at a similar number from the simple considerations in section |4.3| , 
assuming 3.5 hadrons per unit rapidity with an average pt of 0.35 GeV . The ex- 
tra Et can be attributed to hard and soft gluon radiation. The models for these 
processes varied a lot in their predictions, but have evolved since. 

In fig. 1573 the latest preliminary HI data on energy flows in the hadronic CMS 



115| are compared to the model LEPTO 6.4 with different treatments for per- 
turbative QCD radiation. Without QCD radiation, a flat rapidity plateau with 
Et ~ 0.7 GeV per unit rapidity is expected. QCD radiation as expected from the 
hard matrix element generates Et in the current hemipshere only. It is shown below 
that the Et in the current hemisphere depends on Q 2 . With parton showers only, 
a reasonable description of the data is achieved. For this observable, the parton 
shower covers already most of the relevant phase space, including the matrix ele- 
ment. The best description of the data is obtained by combining radiation according 
to the LO matrix element with parton showers. 

The energy flow has also been studied as a function of Q 2 (fig. |5.3| ). At low 
Q 2 the rapidity plateau is almost flat. With increasing Q 2 , the Et in the "photon 
fragmentation region" (loosely defined as far in the current hemisphere) increases 
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Figure 5.2: The E T flow from HI |yj at (x) w 0.002 (left) and 0.0004 (right) for (Q 2 ) w 
14GeV 2 . The data are in the hadronic CMS (proton direction to the left). Shown are statistical 
errors only. The data are compared to LEPTO 6.4 with different treatments of perturbative QCD 
effects: QPM: no radiation; ME: matrix element only; PS: parton showers only; ME+PS: matrix 
element plus parton shower (default). 



considerably. In that region, the virtuality of the photon governing hard interactions 
plays an important role. At central rapidity, the Q 2 dependence is comparatively 
small. This supports qualitatively the picture that particle production at central 
rapidity is independent of the type of the colliding particles, be they ir,K,p or a 
virtual photon 7* ||129| . 
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Figure 5.3: Transverse energy production in the CMS. Et as a function of pseudorapidity r] 
by HI [130, 115]. The DIS data for different Q 2 bins are compared to photoproduction data at 
Q 2 ss 0. The proton beam direction is to the left. Shown are statistical errors only, except for the 
two foremost data points measured with the plug calorimeter. 



In sect. K.2I the Et flow is discussed in more detail in connection with low x 



physics. Further data on hadronic energy production are available from HI, namely 
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energy-energy correlations | |113| , 115|, and the distribution of transverse energy at 
central rapdity ||1 15 ]. 



5.2 Charged Particle Multiplicities 

Average multiplicity 

The multiplicity of charged particles in the CMS current hemisphere (xp > 0) 
has been measured by HI ||116|| with their central tracking chamber for events^ at 



(Q 2 ) « 23 GeV 2 in the kinematic range 80 GeV < W < 220 GeV (fig. |Ja), 
thus reaching fragmentation energies not yet accessible at LEP. From phase space 
arguments alone, one would expect (n) oc In W, (see section |4~3l). Scale breaking due 
to QCD radiation should lead to a faster rise of the multiplicity, however dampened 
by colour coherence and the running of a s . A faster than logarithmic rise is observed 
when comparing the HI data to lower W lepton-nucleon scattering experiments. 
The QCD model LEPTO somewhat overestimates the multiplicity (see fig. |5.4| ). 
The DIS data are also compared to the JETSET prediction (light quarks only) 
for e + e~data, which is known to give a good representation of the e + e~data. The 
behaviour of the total multiplicity in the CMS current hemisphere in DIS is similar 
to e + e~ annihilation data. Differences are found however for the distribution in 
rapidity, see below. 

The growth of multiplicity with energy has been calculated via the QCD evolu- 
tion of fragmentation functions for running a s in NLO |)9, 131|| , yielding 



(n) 



a s (W 2 



,c/y/a a (W*) . 



1 + dJa s (W 2 



(5.2) 



The two-loop expression eq. [2.8| for a s is used with W 2 as scale. The constants b 
and c are given by the theory as 



, 1 10n/ 

b — — I L 

4 27A, 



and 



v/96 



7T 



A) 



(5.3) 



with rif = 3 active flavours and (3q = 11 — (2/3)n/. a and d are free parameters. 

Predictions for hadron spectra can also be derived from the parton multiplicity 
in a jet, which is calculable in QCD by summing up all the branchings in the 
cascade down to a cut-off Qq. The hypothesis of local parton-hadron duality (LPHD) 
(see sect. [4.5| ) is made for the transition from partons to hadrons: the hadron 
spectra should have the same shape as the perturbative prediction, merely differ by 
a normalization constant. The prediction will depend on the cutoff Qo, because the 
charged particle multiplicity is not an infrared safe observable, see section [2.6| . If one 
assumes that the shower stops at a scale given by a typical hadron mass Qo ~ m n , 



1 Rapidity gap events had been excluded from this analysis. 
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Figure 5.4: a) The average charged multiplicity (n), b) the ratio {n)/D, and c) the factorial 
moment i?2 for charged particles in the CMS current hemisphere as a function of W. The DIS 
data are compared to the DIS QCD model MEPS, to the e + e~modcl JETSET for u,d,s quarks, 
and to a MLLA fit, and to the LO and NLO predictions for R 2 . 
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the Qo dependence is eliminated. One can then compare the perturbative prediction 
directly to the data. 

Most important for the average multiplicity is the growth at small fractional 
momentum due to soft gluon radiation. Quantum mechanical interference leads to 
a suppression of soft gluon radiation (soft colour coherence). The shower evolution 
in the leading-log approximation (LLA) is modified (hence called MLLA) to take 
into account destructive interference of soft gluons ||90|| . The modification consists 
of changing the evolution variable such that angular ordering |132| | between 
subsequent emissions i and i + 1 is fulfilled. That is, their opening angles satisfy 
6i > 9i + 1 (see section 



The MLLA+LPHD prediction for the hadron multiplicity for running a s is [90, 



98 



with 



(n) = cx-Nla + c 2 , 
y 



N LA {Y) = T(B) (-) h- 



b{z) 



(5.4) 



(5.5) 



where z = sJ^Y, Y = In(W/2Q ), a = ll + 2n//27 and B = a//3 . T is the Gamma 
function and I u the modified Bessel function of order v. Q = 0.27 GeV is used 
| 98fl . Eq. |57J] reduces to eq. |57J| for large z. 

The data can be fit by either eq. 
faster growth than eq. |5\ 



74] or eq. |h2 
namely a power law [B 



For fixed a s one would expect a 



(n) = a(W/W ) 2b ' - c 



(5.6) 



which is still consistent with the data. In fact, all these fits would be indistin- 
guishable in fig. 0. The fit with eq. §3 yields a = 0.034 ± 0.005 and A = 
In this fit the correction term d^/a s has been neglected, be- 
0.2 ±0.3). It has to be pointed out that A 



0.190 ± 0.060 GeV. 
cause it had been found insignificant id 



extracted from eq. [572] need not be identical with A-^^r , though they are expected 
to be similar if d is small 111 3 111. The fit results are summarized in table |57T. 



ansatz 


fit parameters 


fixed parameters 


power law, eq. 5.6 
MLLA+LPHD, eq. 5.4 
NLO resummed, eq. 5.2 
same, but A free 


a = 1.40 ± 0.04 
a = 1.21 ±0.05 
a = 0.041 ±0.006 
a = 0.034 ±0.005 


b' = 0.20 ±0.01 
c 2 = 0.81 ±0.08 
d = 0.2 ±0.3 
A = 0.190 ±0.060 GeV 


c = 0.5, W = 1 GeV 
Qo = 0.270 GeV 
A = 0.263 GeV 

d = 



Table 5.1: The parameters of fits to the dependence of the average charged multiplicity 
the CMS current hemisphere as a function of W [H6|. 



,n) m 



61 



Multiplicity Distributions 



HI has also measured the multiplicity distribution P n , giving the probability to 
observe in a given event n charged particles, as a function of W, Q 2 and for varying 
pseudorapidity ranges. In general, the data can be parametrized with the Negative 
Binomial Distribution (NBD), 



P n {k,n) 



k(k + 1) • • • (k + n- 1) 



nl 



n 



n + k 



k 



n 



k 



(5.7) 



The two parameters k, n can be expressed in terms of the average multiplicity and 
the dispersion D: 



((n — (n)) 2 ) (n) = n 



D 2 



(n) 



k' 



(5.8) 



One arrives at the NBD in many phenomenological models ||133|| for hadron pro- 
duction. A large body of data can be parametrized by the NBD. In QCD the LLA 
also leads to the NBD for the gluon multiplicity in a quark jet [ |134|| . An alternative 
function is the Lognormal Distribution (LND), derived for multiplicative branching 
processes [135|. A variable is distributed according to the LND, if its logarithm 
is Gaussian distributed. The HI data (fig |5~5| ) are equally well described by both 
parametrizations. The MEPS model deviates from the data both for very small and 
large multiplicities. In particular the origin of the excess of 0-prong events in the 
data over the MEPS model is not clear. For fixed W, a Q 2 dependence of (n) or D 
in the range 10 GeV 2 < Q 2 < 1000 GeV 2 was not found, within errors. 

Hard gluon radiation would lead to a superposition of different NBDs, hence 
deviations from a single NBD for the hadron multiplicity distribution. A multi-jet 
induced shoulder structure as observed at LEP [|136j| is however not visible in the 
HI data. 

In fig. |5.6| the multiplicity distribution is plotted in the scaling KNO (Koba, 
Nielsen and Olesen) form ||137|| , (n)P n as a function of n/(n). In a broad class of 
production mechanisms, based upon scale invariant stochastic branching processes, 
(n)P n is expected to depend only on nj (n), and not on energy [ |138| |. The HI data 
for pseudorapidity 1 < i] < 5 for different W bins exhibit KNO scaling, and agree 
with the JETSET result for e + e~ annihilation. 



In fig. |5.4| b and c the energy dependence for the ratio (n) / D and the normalized 
second order factorial moment [ 139|| R2 ■= (n(n — 1))/ (n) 2 are shown (further data 
on the moments of the multiplicity distribution can be found in ||116|| ). For strict 
KNO scaling, they should be energy independent, which is the case for HERA 
energies at large W, but not for data at smaller W. The data are rather well 
reproduced by the JETSET e+e" model, and by the MEPS DIS model. A QCD 
calculation approaches the R2 measurements when NLO corrections are taken into 
account ||140|| , but deviates still significantly from the data. The LO result in the 
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Figure 5.5: The charged particle multiplicity distribution P n for 115 GeV < W < 150 GeV in 



indicated CMS pseudorapidity domains. The data [116] are compared to LEPTO (MEPS 6.4) and 
to the NBD and LND fit results. The data for 1 < rf < 5 are shown at the true scale, the others 
are scaled down by factors of 10 with respect to the previous one. Shown are statistical errors only. 



double logarithmic approximation (DLA), R2 = 7/4, is reduced in NLO[|l4Tj|] to 



R2 = (1 — xV^s) ' 7/4, where x ~ 0.55. Missing higher order effects, which are taken 
into account with parton showers in MEPS, may be responsible for the disagreement. 

Rapidity Distribution 



HI has measured the charged multiplicity as a function of pseudorapidity ||118|| for 
different kinematic regions at low Q 2 , see fig. fT7\ . In the plateau region about 
2.5 charged particles per 7] unit are observed. The QCD models describe roughly 
the data, but are far from perfect. In particular HERWIG overshoots the data at 
low x towards the central region. We mention already here that quite substantial 
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Figure 5.6: The multiplicity distribution [|ll6| for 1 < 77 < 5 in KNO form for different W 
values, and compared to the JETSET prediction for e + e~ annihilation. Shown are statistical 
errors only. 

differences between the models and the data emerge at small x when only hard 
particles (pr > 1 GeV) are considered (between 0.2 and 0.3 such particles are 
found per i] unit). This effect will be discussed in detail in connection with low x 
physics in section |8.3| . 



The charged particle density measured by HI [|116|| close to the central region 
1 < i] < 2, (n) ~ 2.5 per rj unit, is consistent with data from hadron-hadron 
interactions, see fig. |5.8| . When compared with other DIS experiments, a clear 
increase of the particle density with W is seen. The increase of the total multiplicity 
with W is not entirely due to the increasing longitudinal phase space (rapidity), but 
also due to an increase of the number of particles per unit rapidity, i.e. the height 
of the rapidity plateau increases with W. The average E T at central rapidity (see 
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Figure 5.7: The CMS pscudorapidity distribution of charged particles [US]. The proton direc- 
tion is to the left. Data are shown for nine different kinematic bins with the mean values of x and 
Q 2 as indicated, plus the combined sample (bin 0). The models ARIADNE 4.08, LEPTO 6.4 and 
HERWIG 5.8 are overlayed. 
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section |S.2| ) in DIS and in hadron-hadron collisions increases with a similar slope 

(fig- ED 

The particle density at HERA is significantly lower than for e + e data, repre- 
sented by the JETSET result. Also the W dependence in e + e~ is steeper than in 
hadron-hadron interactions. The DIS data appear to behave as one would expect 
from hadron-hadron data. That is also true for the transverse energy (see section 
T|) measured at central rapidity ||130|| , see fig. [5T8 . 



Why are there for the same CM energy W less charged particles in DIS than 



in e e~ at central rapidity? A possible explanation is the "antenna effect" ||141 
that leads to suppressed radiation from the smeared out colour charge of the proton 
remnant. Another possible explanation is that scaling violations leading to larger 
multiplicities are stronger in e + e~ interactions than in ep reactions due to the larger 

scale, Q 2 + _ = W 2 + _ = s„+„- 3> Q 2 . The xp spectra in the CMS and their 

e e e e e e ^ e P ^ 

scaling violations are presented in the next section. 



5.3 Charged Particle Momentum Spectra 

In this section inclusive charged particle spectra measured by HI and ZEUS ||113 



148|| with their central tracking chambers [|113j , |148| are presented. This limits the 



acceptance to the CMS current region, xf > 0. Measured are longitudinal and 
transverse momentum spectra with respect to the virtual photon direction (fig. |5.9| ). 
The kinematic region covered is 10 < Q 2 < 100 GeV 2 and 55 < W < 200 GeV for 
HI, and 10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV for ZEUS. 

xf spectra in the CMS 

The Feynman-x (xp = 2p z /W) spectra are steeply falling with xp , in fact more 
steeply than what would be expected from the naive QPM, that is quark fragmen- 
tation without QCD radiation (fig. [5.10| a). The data can be described when QCD 
radiation is taken into account. Due to gluon radiation, more than one parton share 
the available energy for fragmentation; therefore less fast particles are produced, 
and the multiplicity of soft particles increases (fig. |4.11| ). 

The spectra can in principle depend on the available phase space, determined by 
W, and on the virtuality Q 2 . A dependence on x ~ Q 2 /W 2 is implied and important, 
because the parton composition of the proton changes with x. The following table 
5T2| gives the scales involved. The DIS xp spectra measured at HERA ||113| , |148| 



are softer than the spectra measured at lower energy fixed target experiments by 



EMC [|149|| and E665 | 150 |, see fig. |5.11| b. From this comparison it is however not 
clear whether the difference is due to the different values of x, W or Q 2 . 



In fig. [5. llj a the HERA data are compared with LEP data from DELPHI ||151 



(divided by 2 to account for the two hemispheres) at similar, though somewhat 
smaller centre of mass energy ^Js e + e - = W = mz- For xf > 0.1 fragmentation 
univerality appears to hold, since the HERA and LEP spectra agree. At smaller 
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Figure 5.8: Solid symbols: W dependence of charged particle density ((n) = # charged parti- 
cles/unit pseudorapidity, left scale). Open symbols: W dependence of the mean central transverse 
energy (Et ) per unit pseudorapidity (right scale). The dotted line is the JETSET result for (n) 
in e+e~ reactions, and the full curve a parametrisation (n) = 0.35 + 0.74(W^ 2 ) 0105 [pi. The DIS 



multiplicity data [116, 143, 144 1 and transverse energy data [ 130 1 are compared to hadron-hadron 
collisions ||| p| |w|, H^. The DIS E T data are discussed in section B.2. 
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Figure 5.9: The longitudinal (p z ) and transverse (pt ) momentum components of a hadron h 
in the 7*p CMS. 



xf, the multiplicity is smaller in DIS than in e + e~ interactions. (The different 
energy dependence of multiplicities at central rapidity in DIS and e + e~is discussed 
in section f>.2[ see fig. |5J|.) The following circumstances may be held responsible for 
the difference: 1) mass effects due to the different flavour compositions; 2) the BGF 
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Figure 5.10: Charged particle distributions in the CMS current hemisphere from HI [ 127 1 and 
ZEUS 148]. a) The Feynman x (= Xp) spectrum, b) The pt spectrum for xp > 0.05. c) The 
(pj.) vs. xp (seagull plot for Xp > 0). The data are compared with the QPM (dotted line) without 
QCD radiation, and the QCD models MEPS (LEPTO 6.1, full line) and CDM (ARIADNE 4.0, 
dashed line). 



contribution to DIS, which is absent in e + e~ interactions; 3) the scale Q 2 , which 
is smaller at HERA than at LEP; 4) the antenna effect due to the extended colour 
charge of the proton remnant (though xf = 0.1 is already 2 units of rapidity away 
from the target region 77 < 0). Further studies are needed to identify the reason(s) 
for the difference. 

A recent analysis by E665 [ 152 ] of the W and Q 2 dependence of the xf spectra 
in the range 7.5 < W < 30 GeV and 0.15 < Q 2 < 20 GeV 2 reveals that at low 
Xf they depend mainly on W, and at large xp mainly on Q 2 . The increase with 





(W) (GeV) 


{Q 2 ) ( GcV^ ) 


(GeV) 


(x) 


LEP e+e-(s/s = m z ) 


91 


8300 


91 




HERA ep 


120 


28 


300 


10~ 3 


fixed target IN (EMC/E665) 


14/18 


10 


24/30 


5 • lO" 2 



Table 5.2: Approximate kinematics and scales for the e + e , ep and IN data which are discussed 
here. 
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Figure 5.11: a) The x F spectrum from ZEUS (148) compared with DELPHI data p51j . The 
DELPHI data are divided by 2. b) The xf spectra measured at HERA ]113| , |148[ , compared with 
the QPM (dotted line) without QCD radiation, the MEPS model (LEPTO 6.1, full line), and with 



fixed target IN data at lower W from EMC [149| and E665 |15 



W at small xf is essentially due to the increase of longitudinal phase space. The de- 
crease at large xf has been attributed to the scale dependence of the fragmentation 
function. 
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QCD predictions with fragmentation functions 



The hadron production cross section (see fig. |5.12| ) da/dz can be written as a con- 
volution of the parton production cross section with the fragmentation function 
Dh/j(z, which gives the number density to observe a hadron h with momentum 
fraction z from the fragmentation of parton j. The parton production cross section 
itself is a convolution of the parton density function f%/ p (£, [ip) for parton species 
i in the proton with momentum fraction £, with the hard scattering cross section 
&ij to produce parton j from electron scattering on that parton i, e + i — > j + X. 
Finally, the cross sections need to be summed over all parton species 



da 
dz~ 



fi/p(€> ®^ij®D h/j (z^ 2 D ). 



(5.9) 



hi 



The parton densities and the fragmentation functions depend upon their respec- 
tive factorization scales, \j? F and \j? D . The hard scattering cross section <jy = 
^ijif^Fy A*D; A*R; depends in addition upon the renormalization scale n R and a s . 
The scales fi R , fi 2 D , \j? F are in principle arbitrary; here they are taken to be Q 2 (not 
W 2 !), the physical scale related with the hard scattering in DIS [ 153 1 . The situation 



in DIS with Q < W is different to e + e", where Q = W = y/s = 2E q 





Figure 5.12: a) The process ep — > e'hX. In lowest order the blob labelled Sij represents 
electron-quark scattering via the exchange of a virtual photon (b). 



With experimentally determined parton densities and fragmentation function 
parametrizations from other, mainly e + e~, reactions [|154|| , the NLO QCD calcu- 
lation for the HERA xf spectrum is quite satisfactory (fig. [5.13|) |[L53|| . The LO 
calculation predicts a spectrum that is too hard compared with the data. When 
moving out of the current fragmentation region, i.e. towards small and negative xf, 
these calculations are expected to break down due to fragmentation effects that are 
not related to the scattered parton. 

Since the scaling violation of the parton density and fragmentation functions as 
well as the hard cross section depend on a s , a measurement of the scale dependence 
of da/dz allows in principle a measurement of a s . Due to the Q 2 range covered 
at HERA, it would be possible to measure ot s from the scaling violations of the 
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Figure 5.13: The calculation by Graudenz [153] for different parton density parametrization 
compared to the HERA xf data from a) HI and b) ZEUS. The spectrum is calculated in LO 
with the GRV LO parametrization (full line), and in NLO with the parametrizations GRV HO 
(dashed), MRSA' (dotted) and CTEQ3M (long dashed). 



fragmentation functions in a single experiment. For an integrated luminosity of 250 
pb _1 a competitive statistical error for a s (m 2 z ) of ±0.0007 is expected. The theo- 
retical uncertainties are presently estimated to be much larger ||153|| : ±0.005 from 
the parton densities and ±0.011 from the scale uncertainty. A better knowledge of 
parton densities would reduce the corresponding error. Possibly the theoretical situ- 
ation is more favourable in the Breit system, where the scale of the hard interaction 
also sets the energy of the fragmenting quark, Q/2 — E q , similar as in e + e~ . First 



indications of scaling violations in the Breit frame are presented in section |0 



Transverse momenta 

The transverse momentum spectrum of charged particles from ZEUS with xp > 
0.05 is shown in fig. 5.10| b. The data exhibit a high px tail over the expectation 
from the QPM without QCD radiation, which is well described by the QCD models. 
Pt can be generated in the QCD cascade, where partons are emitted before and after 
the hard scattering. Somewhat arbitrarily one may distinguish radiation from the 
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hard matrix element (BGF and QCDC) and soft gluon radiation (parton shower). 

In general, contributions to the px of primary hadrons are expected also from 
other sources than QCD radiation (p^ d ) ||155| , |5j : from an intrinsic px of the initial 



parton in the proton (p™ tr ), and from fragmentation (p^ ag ). These contributions 
can be summed up|, 

04) « xp(p^ 2 ) + (pp g 2 ) + (pf d 2 ). (5.10) 

The pt of the finally observed hadrons will be diluted due to decays of the primary 
hadrons. At HERA, the relative importance of p r p d is bigger than at the lower energy 
DIS experiments. This is best studied in the so-called seagull plot. 

In the seagull plot the mean p\ of charged hadrons is plotted as a function of 
x F . The data from EMC P5fl , HI fTTJ and ZEUS fn§ are shown in figs. |7Lj and 



5.10 c. The shape of this distribution is reminiscent of a seagull in the skvQ Due 
to the limited acceptance to xf > 0, at HERA only one wing of the seagull plot is 
measured. 

The HERA data are well described by QCD models. Without QCD radiation 
(the QPM) the prediction falls far below the data. While the QPM predictions are 
quite similar for EMC and HERA energies, QCD radiation effects are much larger 
at the higher energy experiments. The (p T ) increases with W [ |148| . 



The fixed target data also cover the target fragmentation region, xp < 0. In 
the target region QCD radiation is suppressed as compared to the current region. 
This suppression can be explained with the smeared out colour charge in the proton 



2 An intrinsic px of a parton in the proton leads to a tilt of the fragmentation axis with 
respect to the 7* axis. A daughter hadron with momentum fraction xf will therefore inherit a 



transverse momentum xf • Px ■ Models assume a Gaussian distribution as in eq. 4.24 with 



{Pr tr ) = (0-39 GeV ) 2 in accord with the uncertainty principle applied to a parton confined 



to the proton volume. The experimental value is (pf tv ) = [0.29i°;o7(stat.)i°;^(syst.) GeV 12 



[149 1 



Pt from fragmentation is also assumed to be Gaussian distributed (eq. 4.24) |92j|. A good 



choice for the description of the HERA and LEP data with the Lund string model is {p, 



(0.3 GeV ) 2 0, y|, [L2§. (p T ag ') = (0.41 ± 0.02 GeV ) 2 has been measured with fixed target 
data p9| . 

3 The dip at if =0 is mainly a phase space effect. Consider massless hadrons distributed 
according to the longitudinal phase space model, 

, d / 2 = a ■ exp(-p|/b). (5.11) 

Because dy/dp z = l/E, we have 

d 2 a _ W/2 d?a 
dxFdp T E dydp T ' 

For small longitudinal momentum, E = y 'p\ + p T w px- So for xp — * 0, the cross section is 
weighted with a factor that grows with 1/pt for small px, giving a large weight to small px 
hadrons in the average p\. 



2, 
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Figure 5.14: The seagull plot, (p|) vs. x F for charged hadrons in the CMS |157[ . x F > 
defines the current hemisphere, x F < the target hemisphere. The data from HI [ 1 13 1 , ZEUS 
[148[ and EMC pi are compared to a QCD model (CDM, ARIADNE 4.08) and the QPM. 



remnant, if one considers radiation from a colour dipole. Perturbative QCD evolu- 
tion as implemented in LEPTO yields a similar asymmetry. The description of the 
EMC data by the CDM (ARIADNE 4.08) is certainly not satisfactory. It would be 
interesting to try to tune the QCD models to both data samples simultaneously. 

At HERA energies an increase of (p T ) with Q 2 is observed at fixed xf, whereas 
at the lower energy fixed target data from EMC ||149|| almost no Q 2 dependence was 
found, see fig. |5.16| . With very precise data in the kinematic region 7.5 < W < 
30 GeV , 0.15 < Q 2 < 20 GeV 2 and 1.5-10" 4 < x < 0.6 the up experiment E665 



has recently seen a W as well as a Q 2 dependence for the mean p\ [|152 



We can use a QCD model to investigate the different contributions to the soften- 
ing of the xf spectrum and to the increasing p?. In the MEPS model (LEPTO 6.4) 
radiation from the LO matrix element at the photon vertex (BGF and QCDC) can 
be identified as the main cause; additional soft gluon radiation is relatively unim- 
portant (fig. |5.15| ). On the other hand, the QCD parton shower as implemented in 
LEPTO appears to cover already a large part of the LO matrix element contribution. 
Note that the parton shower changes when the matrix element cut-off changes (or 
is switched off) in order to cover the maximal phase space without double counting 
(this adjustment is sometimes called "matching"). 



5.4 Strangeness Production 

Strange particles may be produced in DIS when the scattered quark is a strange 
quark (either from the sea, or produced in BGF events), of from ss pairs created 
during parton showering and hadronization, from heavy quark decays, or from non- 
standard processes like instantons. Strange particles produced directly from the 
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Figure 5.15: QCD radiation effects according to the LEPTO 6.4 model JL58| . The data are 
compared to the LEPTO model without QCD radiation (QPM), QCD radiation from the LO 
matrix element alone (dotted, ME), QCD radiation from parton showers alone (PS, dash-dotted), 
and the full model (full line). 



scattered quark are expected to carry a relatively large fraction of the scattered 
quark's energy, whereas most of the strange particles are expected from hadroniza- 
tion with small fractional energy. In the Lund hadronization model, qq pairs are 
created in the colour field of the string from the vacuum. Because m s > m u « rrid, 
ss pair creation is suppressed with respect to the other light quarks with a ratio 
A s : 1 : 1. Typically, A s = 0.2 — 0.3, with 0.3 the Lund default value representing 



the experimental mean from [159 



HI ||117[ and ZEUS | |160|| identify and reconstruct K° and A particles via the in- 



variant mass of their decay products in the decays — > n + n , A 



A 



P7T 
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Figure 5.16: The Q 2 dependence of the charged particle (p 2 ^) for different values of xp . The 



ZEUS 1 148] data at (W) = 120 GeV are compared to EMC data [149| at (W) = 14 GeV . Also 



shown are the MEPS predictions for the two values of W. 



In the following, the notations K° and A include both particles and antiparticles. 
The decays are measured in the central drift chambers. This limits the acceptance 
for strange particles in the laboratory system to —1.3 < 77 < 1.3 and px > 0.5 GeV. 
The results for K° and A yields are given in table |5.3| . Neglecting the slightly differ- 
ent kinematic selections and averaging the ZEUS and HI results, for —1.3 < rj < 1.3 
and pt > 0.5 GeV the K° yield is 0.111 ± 0.005 per event and unit of pseudora- 
pidity, and for A baryons it is 0.042 ± 0.0005. 





HI 


ZEUS 


event selection 


no large rapidity gap 
10 GeV 2 < v\ < 70 GeV 2 
0.0001 < x < 0.01 
0.05 < y < 0.6 


10 GeV 2 < v\ < 640 Gey2 
0.0003 < x < 0.01 
0.04 < y 


strange 

particle selection 


0.5 GeV <p T < 2.12 GeV 
-1.3 <r) < 1.3 


0.5 GeV <p T < 3.5 GeV 
-1.3 < t] < 1.3 


# K v /event 
#X° / # charged * 
A event 


0.287 ±0.008 ±0.012 
0.058 ±0.002 ±0.003 
0.053 ±0.007 ±0.007 


0.289 ±0.015 ±0.014 
0.077 ±0.006 ±0.008 
0.038 ±0.006 ±0.002 



Table 5.3: Strange particle production at HERA. A lower p T cut of 0.15 GeV for HI and 0.2 
GeV for ZEUS is applied for the charged particles in the ratio. * For the ratio #if / # charged 
both experiments exclude events with a large rapidity gap. 
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The strange particle yield has been studied as a function of x, W, Q 2 , r], pr and 
xp. In general, the same features are observed as for charged particles. These 
are an almost flat distribution in pseudorapidity, steeply falling pt spectra, and no 
significant dependences on the event kinematics other than on W. The ratio of K° 
to charged particle yield was found to be the same in events with and without a 
large rapidity gap, within errors. The average K° multiplicity in the CMS current 
hemisphere increases logarithmically with W, see fig. |5.17| . This behaviour is similar 
to what has been seen for charged particles, though a faster than logarithmic rise is 
not seen with K° mesons (compare |5.2| ). 
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Figure 5.17: The average K° (top) and A (bottom) multiplicity in the CMS current hemisphere 



as a function of W 2 . The HI data [117] are compared to fixed target data from E665 [ 1 6 1 1 and 



WA21 [162| 



The data can be reasonably well described by standard QCD models invok- 
ing Lund string fragmentation, provided X s ~ 0.2 is chosen. As an example, the 
xf spectrum and the seagull plot for K° mesons are shown in fig. |5.18| . The larger 
ipj) observed when going from fixed target to HERA energies can be understood 
with QCD radiation in the same manner as for charged particles. About the same 
(pj) are observed for K°s and for charged particles. It is not yet clear why the 
supposedly scaling xp spectra of the fixed target data lie below the HERA data for 
xf > 0.1, but this is also expected from Monte Carlo | 163|| . Probably mass effects 
play a role at smaller energies, in contrast to the HERA data. The HI data agree 
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with the expectation for e + e annihilation to light quarks at the same CM energy, 
represented by the JETSET generator. 



5.5 Charm Production 

The charm yield 

Charm production in DIS has been investigated by detecting D* + — ► D°tt + — > 
K~7t + tt + 1 166, |167|| and D° — ► K~n + decays (and their charge conjugates). The 



decay products are detected in the central drift chambers. D° mesons are recon- 
structed by the Ktt invariant masses. D* mesons are identified by the kinematically 
tightly constrained mass difference Am = m(D°Tr^ ow ) — m(D°) between a D° — > Ktt 
candidate and the slow pion, and the D° candidate mass. 

Apart from possible anomalous sources of charm (new particles, instantons, ...), 
the main interest in charm production stems from its sensitivity to the gluon content 
of the proton at small fractional momentum x g . The two lowest order standard 
charm production processes in DIS are boson gluon fusion (BGF) and scattering 
off a charm sea quark in the proton (QPM)Q, see fig. |5.19[ At high x also an 



intrinsic (valence) charm component to the proton at the few permil level has been 
discussed [ |168fl . At small x the BGF contribution is expected to dominate due to 
the increasing gluon content. When the invariant mass \/§ of the qq system is large 
enough (I ^> (2m c ) 2 ) so that mass effects can be neglected, the quark flavours are 
expected to be produced according to u : d : s : c = 4/9 : 1/9: 1/9:4/9 due to their 
electric charges, provided one is below the b threshold, bb production is suppressed 



by the large bottom mass. It has been estimated ||169|| that in DIS b : c production 
is »s 0.02, integrated over a large range in x, Q 2 , but can be as large as 0.1 at small 
x and large Q 2 , that is large W. There are no results on open b production from 
HERA yet. 

The charm cross section^ is inferred from the D production rate by taking into 
account the charm fragmentation function and the appropriate D branching ratios, 
and by extrapolating from the px and r\ regions accessible to the measurement. The 
HI and ZEUS measurements cover 5 GeV 2 < Q 2 < 100 GeV 2 and are summarized 
in tab. [5.4| . They are in reasonable agreement with each other, and with NLO 
calculations [ |171|| based on current parton density parametrizations. About 4% of 



the DIS events contain a D*^ ||167|| , and about 25% of the DIS events are due 



to charm production [ 166 , 167 . The asymptotic value of 40% expected from the 
quark charges is not yet reached. HI finds for the ratio of D* : D° production 
0.38 ± 0.07 ± 0.06, in agreement with other experiments (see compilation in [ |166|| ). 
Note that the D°s include non-primary D°s from D* decays. 

4 The boundary between QPM and BGF is subject to definition. 

5 The elastic (diffractive) cross section for ep — > epJ/tp is comparably small, 100 ± 20 ± 20 pb 



for Q 2 > 8 GeV 2 and 30 GeV < W < 150 GeV [170| 
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Figure 5.18: Top) xf spectrum of K° mesons. The HI data [117] are compared to fixed target 



data from EMC and E665 |164], |161[ , and to the DIS event generators CDM, LEPTO and HERWIG. 
For CDM and LEPTO the DELPHI fragmentation parameter set |L65| with A s =0.23 were used. 
The curve JETSET is a simulation of e + e~ annihilation (u,d,s quarks only) at ^/s = 138 GeV. 



Bottom) The average -p\ as a function of xf for K° [117] and charged particles [113] from HI, 



compared to fixed target data [164, 161 and QCD models 
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Figure 5.19: Charm production in LO: a) boson gluon fusion (BGF) and b) from sea or valence 
quarks (QPM). In the BGF process the fractional gluon momentum x g is related to Bjorken x by 
x g = x(l + s/Q 2 ). 





HI 


ZEUS 


NLO 


er C c(nb) 


0.01 < y < 0.7 


< y < 0.7 




5 GeV^ < Q 2 < 10 GeV^ 
10 GcV 2 < Q 2 < 100 GeV 2 


17.4 ± 1.6 ± 1.7 ± 1.4 


13.5 ± 5.2 ±1.8iJ;!j 
12.5 ± 3.1 ±1.8iJ-j 


9-13 
8-14 




8-10- 4 < x < 8-10~ a 


2-10- 4 < x < 5-10~ 3 




{Ff/F 2 ) 


0.237 ±0.021+^« 


m 25% 





Table 5.4: Charm production at HERA [166, 167 . The extrapolation of the HI measurement 
from the range 0.01 < y < 0.7 to < y < 0.7 would modify the measured cross section by 
approximately +6%. The NLO calculations [171|, depending on the assumed gluon distribution 
function, and on the value of m c used, vary by the range given. 



Momentum spectra of D mesons 



In order to investigate further the charm production mechanism, momentum spec- 
tra of the D mesons in the hadronic CMS have been measured ||166| , |167|| . The 



Pt spectrum is well described by the Monte Carlo generator AROMA ||1 72|| based 
upon the BGF process (see fig. |5.20| a), and also by NLO calculations [ |167|| . The 
spectrum of the scaled momentum variable := 2\p\/W falls with increasing 
(see fig. |5.20| b), in agreement with the BGF based charm generator, and also with 
the NLO calculations [ 167 |. In contrast, QPM-like events from charm sea quark 
interactions would lead to an xd spectrum that is peaked at xd ~ 0.6, reflecting 
the hard charm fragmentation function. HI excludes at 95% C.L. a 5% contribution 
from charm sea quarks to the total charm production. 



The charm structure function F£ 



The cross section measurements in bins of x and Q 2 can be expressed in the form of 
a charm structure function F% c (fig. |5.21|) analogous to F 2 . When compared with 



79 



2 4 6 8 o 10 

p t D °[ GeV ] 

a) 

Figure 5.20: a) p T spectrum of D° and D* mesons fl66f , compared to the AROMA 2.1 [Jl72| 
expectation for m c = 1.3 GeV (full line) and m c — 1.7 GeV (dashed line). Only statistical errors 
are shown, b) xd distribution for D° (open points) and D* + (full points) mesons in the laboratory 
pseudorapidity range —1.5 < 1)d < 1-5. The data are compared to the expectation for BGF 
processes (AROMA, shaded), for QPM-like contributions from a charm sea (dashed histogram), 
and an extrapolation from charged current vN scattering off strange sea quarks |173}] (full line). 

EMC data at x — 0.02 — 0.3, F| c rises towards small x. The data are consistent with 
NLO calculations based exclusively on the BGF process with gluon densities that 
describe the current HERA Fi data. In the calculation both the renormalization and 
factorization scales are chosen to be /j, — \JQ 2 + 4m 2 . The dominant uncertainty in 
the calculation arises from the charm quark mass. The predicted F^ changes by 
15% when m c is changed by 0.2 GeV. 

The measurements establish the boson gluon fusion process as the dominant 
one for charm production in DIS at Q 2 < 100 GeV 2 , and thus the sensitivity to 
the gluon distribution function in the proton. The data are consistent with parton 
density parametrizations which were derived from the inclusive F2 measurements, 
providing an independent cross check of the interpretation of the F2 data in terms 
of parton densities. 

The gluon density from charm production 

In a new preliminary analysis [ |175|| to determine the NLO gluon density from charm 
production HI measures the visible D* cross section <t^ s (corrected for detector 
effects) in the range 2 < Q 2 < 100 GeV 2 , 0.01 < y < 0.7 for D* mesons with 
laboratory pt > 1.5 GeV and \r)\ < 1.5. The measured cross section is a®* = 
(5.63±0.66to;gg) n b- The measured differential D* cross sections as functions of the 
kinematic variables x , Q 2 , pr, rj (not shown) are well described by the AROMA 
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172 charm generator and by a NLO calculation (program HVQDIS 1.1 [ |176|| with 
the Peterson charm fragmentation function ||177|| ). 

Xg/ I} is the proton momentum fraction carried by the gluon entering the BGF 
process. x g / p can in leading order be reconstructed from 



x 



g/p 



x(l + — ) with 



Pfc + ™ 2 C 

z(l-z) 



and 



P-c _ (E-p z 
P ■ q 2yE e 



\lab 



(5.13) 



Here c is the charm quark's 4-vector and p^ c it's CMS transverse momentum. P 
is the proton 4-momentum, and E e is the electron beam energy in the laboratory 
frame. The Lorentz invariant z can be calculated from the charm quark's energy 
E and longitudinal momentum p z in the lab. frame. For BGF events it can be 



correlated with the quark scattering angle in the photon-gluon CMS, see section \7A 
and eq. f77f| . 

Measured are charm mesons though, not charm quarks. Therefore the observable 
is defined by replacing p* Tc with 1.2 • p* TD * and [E — p z ) l ^ h with (E — p z ) l '$l in 



x 



obs 



eq. |5.13| . x° bs is well correlated with x g / p ||175|| . 



The measured cross section as a function of x° bs is shown in fig. [5.22| a. It is 
well described by the NLO calculation. The NLO calculation takes into account the 
gluon initiated processes j*g — > cc, ccg and the quark initiated processes 7*g — > ccq, 
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Figure 5.21: The charm contibution F^ c to the proton structure function F2 from HI [16£] and 



ZEUS [|T63 . Results from EMC [|174| are shown as crosses. The shaded band gives the NLO QCD 
prediction from the GRV NLO gluon density for charm masses between 1.35 GeV (lower bound) 
and 1.7 GeV (upper bound). 
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7*<? — > ccg, where the cc pair is produced from a radiated gluon. Using the NLO 
calculation for the x° hs cross section, the gluon density can be unfolded from the 
data. In the calculation m c = 1.5 GeV is used, and the factorization scale is set 
to fMp = \/^ m c + Q 2 - The resulting NLO gluon density xg(x,fi 2 ) (fig. |5.22| b) at an 
average scale fi 2 = 24 GeV 2 reaches down to x g / p ~ 10~ 3 , and agrees well with a 
recent indirect extraction from the F 2 data HTJ. Smaller values of x g / p could be 
probed with a larger angular acceptance for D* decays in the backward region, to be 
provided with silicon trackers. The ultimate limit is given by x g / p > s/(sy) ~ 10~ 4 . 

By reconstructing event-by-event x g / p , this measurement probes the gluon den- 
sity locally in x g / p . In contrast, with F 2 CC the charm contribution is measured as a 
function of Bjorken x, integrating over all x g / p > x. The comparison of the directly 
measured gluon density with the one extracted from the inclusive F 2 , and the agree- 
ment of F^ with the prediction based upon such gluon distributions constitutes an 
important test of QCD with universal parton (here: gluon) densities. A similar test 
was obtained with dijet events in LO, see fig. |7.5| . 
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a) The visible D* cross section as a function of x° ha (see text), corrected for 



a) 

Figure 5.22: 

detector effects. The preliminary HI data |l75| (data points) are compared to a NLO QCD 
calculation [176] using the GRV 94 HO |l9| parton density function (histogram; the shaded bands 
represents the variation of the cross section when m c is varied between 1.3 and 1.7 GeV). b) The 
NLO gluon density x-g(x) at a scale = 25 GeV 2 , unfolded from the D* cross section |175[ |. The 
statistical errors are shown as error bars, the systematic errors as shaded histograms. The line 
shows the CTEQ4F3 parametrization p7|, the shaded band the gluon density from a NLO QCD 



analysis of the HI Fi data |41 
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5.6 Bose-Einstein Correlations 



Due to Bose-Einstein statistics, identical bosons prefer to occupy the same quantum 
state. In particle physics, the Bose-Einstein effect leads to an enhanced probability 
for identical bosons - like sign pions - to have similar momenta ||178||. The shape 



and strength of their correlation function provide information on the production 
process. 

As a measure of distance in momentum space for two particles with 4-momenta 



Pi and P2 with invariant mass M = J (pi + P2) 2 one definesf] 

T 2 := -(pi -P2) 2 = M 2 -4m 2 , (5.14) 

assuming the particles to be pions. One measures the number of like sign charged 
particle pairs as a function of T. This correlation function p(T) is compared to a 
function p Te f (T) obtained similarly, but from a reference sample where Bose-Einstein 
correlations are absent: 

R(T):=p(T)/p Ici (T). (5.15) 

As reference sample may serve a) unlike sign charged particle pairs, or b) pairs from 
different events (event mixed method). 

If the particle emitting sources are distributed in space-time according to p s (£,)y 
a ratio 

R(T) = R (l + \\p s (T)\ 2 ), (5.16) 

is expected for the Bose-Einstein enhancement. p s (T) is the Fourier transform of 
p s (£) and Rq a normalization constant. The correlation strength A is for completely 
coherent emitters, and 1 for completely incoherent emitters. 

The measurement at small T, where the Bose-Einstein effect is expected, is 
notoriously difficult due to the finite double-track resolution of the tracking devices 
and due to e + e~ pairs from photon conversions. The construction of a reference 
sample from the data introduces a residual bias. This effect is estimated by Monte 
Carlo, and cancelled by building the "double ratio" 

RR{T) = R(T) dat jR(T) MC . (5.17) 



The HI data ||179|| are shown in fig. |5.23| . A clear enhancement at small T, at- 
tributable to the Bose-Einstein effect, is visible. The data cover 6 < Q 2 < 100 GeV 2 , 
1 • 10~ 4 < x < 1 • 10~ 2 and 65 < W < 240 GeV. No dependence of the Bose-Einstein 
effect on the kinematic variables x, Q 2 or W was found, nor a difference between 
the diffractive and non-diffractive data samples. 

Different parametrizations of the effect have been fitted to the data. The tradi- 



tional Goldhaber parametrization [180] 



R[T) = R [1 + A exp(-r 2 T 2 )J (5.18) 

6 This variable is traditionally denoted as Q 2 . Here we use T 2 instead, because in DIS Q 2 is 
already used for the virtuality of the exchanged photon. 
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Figure 5.23: Bose-Einstein correlations measured by HI ]179[ . The double ratio RR(T) is plotted 
as a function of T and M . For the reference sample, tracks from different events are paired (top), 
or unlike sign pairs are used (bottom). The fits of the Gaussian, exponential and power law 
paramctrizations (with extra linear terms to accomodate the sloping background) are shown. The 
open points in the bottom plot are excluded from the fits, because they are in a region which is 
affected by correlations due to K® — » tt + tt~ and p a — > 7r + 7r~ decays. 



is derived from the the assumption of a static, spherically symmetric Gaussian source 
density p s (£) = p s (0) exp(— £ 2 / (2r 2 )). The parameter r gives the size of the produc- 
tion volume. In high energy reactions, the sources are moving relativistically ||181 



and the assumption of a static source becomes invalid [|182|| . In the relativistic string 



fragmentation picture |183| , an approximately exponential shape is predicted, 



R(T) = R [l + Aexp(-rT)], 



(5.19) 



where the parameter r is related to the string tension |184|| . Alternatively, assuming 
a self-similar, scale-invariant pattern for perturbative QCD cascades, a power law is 



derived [185, 139 



R(M) = A + B ■ (1/My. (5.20) 

Compared to the data (fig. |5.23|) , all parametrizations are acceptable, though the 
exponential and the power law ansatze give somewhat better results at low T. A 
power law behaviour had been observed previously ||186|| . 

In order to compare with other experiments, the following results are obtained 
with the Gaussian (Goldhaber) parametrization, eq. |5.18| (table |Q| ). With the 
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reference sample 


r (fm) 


A 


event mixed 
unlike sign pairs 


0.54±0.03±£$ 
0.68 ± 0.04±°;°j? 


0.32 ±0.02±£$| 
0.52 ±0.031^21 



Table 5.5: The fit parameters from a Gaussian parametrization of the Bose-Einstein effect seen 



in the HI data fl79]. 



event mixed reference sample, a lower r value is obtained than with unlike sign 
pairs. Such a systematic effect has been observed also in other analyses. Other 
measurements in e + e~ and lepton-nucleon interactions scatter between r = 0.39 fm 
and r = 0.97 fm, and A = 0.27 and A = 1.08 (see the compilation in ||179|| ), but 
no pattern emerges. At least the extracted values of r, including the HERA ones, 
can be said to be in the ball park of the length scale at which hadronization takes 
place, 1 fm. HI finds marginal evidence for a dependence of r on the charged 
particle multiplicity density n, given as number of particles per unit pseudorapidity, 
n = dn/drj. A strong multiplicity dependence of r, roughly r = 0.4 + 0.075 • n 
(with 20% accuracy) for 2 < n < 20 was seen in pp collisions 
confirmed by the HI data with a limited lever arm, 2 < n < 5. 



187, 1881. This is 
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Chapter 6 

The Quark Fragmentation Region 



6.1 Charged Particle Spectra 

Comparisons with e + e reactions 

The physics issues that have been addressed with charged particle spectra in the 
Breit frame ||189| , |190| , |191| , |192| , |193|| are the study of fragmentation properties, 



QCD coherence effects, and scaling violations with the goal of determining a s . The 
advantage of ep experiments over e + e~ experiments at fixed beam energies is that 
the scale variable Q can be varied continuously Complications arise however from 
initial state parton radiation and from boson-gluon fusion events, which are absent 
in e + e~. 

With the Breit current hemisphere a relatively small subsystem of the whole 
event is considered, (for most events Q <C W), which should be not much affected 
by the dynamics of the proton remnant (see sections 41 and L2 with fig. |4.3| ). The 



Breit frame current hemisphere of an ep event, the quark fragmentation region, is 
similar to one hemisphere of an e + e~ event. Considering just the lowest order 
process, in e + e~ interactions the outgoing quarks have opposite momenta of equal 
magnitude. The same is true for the incoming quark and the outgoing quark in 
the Breit frame of DIS. Furthermore, as in e + e~, the scale of the hard interaction, 
Q, also determines the phase space for QCD effects, given by the quark energy 
E q = Q/2. In contrast, in the CMS E q = W/2 > Q/2. In the CMS it is therefore 
more difficult to disentangle a kinematic effect from a scale dependent QCD effect. 

The scaled momentum variables of a hadron with momentum p and energy E 
are defined as (compare sect. POI) 



^=JSr e-ln 1 (6.1; 



p"^ x. 



p 



with < x p < 1. In practice, one often approximates p max pa E q = Q/2 (massless 
kinematics). With the variable £ the soft part of the momentum spectrum at small 
x p , where the particle density is the largest, is expanded. 



8(3 



Fig. [TT] shows the measured | |189| | x p and £ distributions of charged particles in 



the Breit current hemisphere for large and for small Q 2 . D(x p ) is steeply falling, 
and F(£) is approximately Gaussian[[ When Q increases, soft particles (large £) are 
produced much more copiously. The hard part of the spectrum (large x p , small £) is 
approximately invariant of Q, "scaling". This can be understood naively. If hadron 
production stems from a branching process with energy independent branchings, 
scaling at high x p is expected. At low x p , the reduction of the phase space forces a 
turnover of the spectrum, which happens at lower x p for larger initial energy. 

The evolution of F(£) can be effectively summarized by the Q dependence of 
the peak position £ pea k and width £ w idth of the Gaussian, shown in fig. |672|. The 
average charged multiplicity (n) (the integral over the fragmentation function) is 
shown in fig |6.3| . An approximately logarithmic dependence of these variables on Q 
is observed: 

(n) « -1.7 + 4.9 • log 10 (Q/ GeV ) (6.6) 
W ~ 0.6 + 1.5 • log 10 (Q/ GeV ) (6.7) 
£width « 0.4 + 0.45 • log 10 (Q/ GeV ). (6.8) 



1 Qualitatively the Gaussian shape can be understood assuming that hadrons result from a 
chain of branching processes i — l,...k Jl94{ ] (gluon coherence giving rise to a hump backed, 
approximately Gaussian, shape is discussed later). The branching stops, when there is not enough 
energy left to produce a hadron of mass m. On average, 

(6.2) 

The Zi give the energy fraction retained in each branching, and are assumed to be randomly 
distributed according to some splitting function Pj(zj). We expect that the hadron multiplicity n is 
proportional to the number of branchings k, hence a logarithmic growth of the average multiplicity 
with energy follows, 

(n) = c-]n(Eg/m). (6.3) 

The energy left after r branchings is 

r 

E r q =E q -z l ---z r = E q T\z i . (6.4) 

»=i 

A hadron from this stage has energy E r = x r ■ E q with x r — z' r E' q /E q , where z' r = 1 — z r+ i. We 
define the random variable £ r := ln(l/x r ), which is 

^ r ■= ln t = ln { v ' n t) = ln v + £ ln h (6 - 5) 

Xr Zi J z r Zi 

t; r can be written as a sum of random variables, and is thus Gaussian distributed for r — > oo 
according to the central limit theorem, regardless of Pi(zi). 

The final hadron spectrum of all hadrons from all k branching stages would then be the result 
of the sum of Gaussians G r (£) with r > k, where k itself is a random variable. For not too large k, 
this would give an approximate Gaussian, resulting from fall-offs at the kinematic boundaries at 
small and large £. For large k (large energy), a plateau F(£) =const can develop, which translates 
into a rapidity plateau dn/dy =const. 
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Figure 6.1: The distributions D(x p ) := l/NAn/Ax p (a) and := l/Ndn/d£ (b) for charged 
particles from the Breit current hemisphere for large and for small Q 2 from HI |19(| . The curves 
are Gaussian fits to the data. 
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Figure 6.2: The evolution of of the peak position £ pea k (left) and width £ w idth (right) of F(£) 
with Q. Data from ep [ 190| , 195 and e + e~ reactions are compared. Also shown is a MLLA fit to 



the HI data |190 



£ pea k and ^ w idth measured in e + e~ annihilation at ^/s e + e - = E* = Q are within 
errors consistent with the ep data. For Q > 10 GeV , also the total multiplicities 
(n) agree, but at smaller Q, (n) is significantly smaller in ep than in e + e~ . 

The deviation for Q < 10 GeV can be attributed largely to LO QCD processes 
in ep, in particular to the boson-gluon fusion process which occur at small x (hence 
small Q) where the gluon density is large. For these events, the scattered quark 
may not be collinear with the photon, or may not even point into the Breit cur- 
rent hemisphere ||196|| . Its fragments are then emitted partly into the Breit target 
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Figure 6.3: Average charged particle multiplicity in the Breit current hemisphere as a function 
of Q. The HI [ 190 1 and ZEUS [ 1 9 1 1 data are compared to a curve representing half of the charged 
multiplicity of e + e~ annihilation events at CM energy E* . 



hemisphere^. For Q > 10 GeV charged particle production in the ep Breit current 
hemisphere appears to be similar to one hemisphere in e + e~annihilation, supporting 
the notion of universal quark fragmentation 

Predictions for the hadron spectra are available from a) perturbative QCD in the 
modified leading log approxiation (MLLA), where the transition to hadrons rests on 
the hypothesis of local parton hadron duality (LPHD); b) perturbative QCD in 
NLO, based on the QCD factorization theorem with universal, but scale dependent 
fragmentation functions; c) event generators with their sophisticated schemes to 
model perturbative QCD evolution and hadronization. They will be compared to 
the data in the following. 



MLLA + LPHD predictions 



The MLLA+LHPD derivation of the spectrum F(£) is discussed in JH], [T57|, [DJB 



The spectrum of partons in a parton cascade initiated by a quark of energy E q 
(and related virtuality) is calculated perturbatively in the modified leading log ap- 
proximation (MLLA), including colour coherence. The cascade is calculated down 
to a low cut-off Qo ~ A QCD for the parton virtuality. It is then assumed (the 
hypothesis of local parton hadron duality, LPHD) that the final state hadrons fol- 
low the same spectrum, up to a constant factor. The predicted spectrum F(£) has 
a "hump backed" form of approximately Gaussian shape around the peak. On 

2 Consider the subsystem of the hard scattering subprocess with invariant mass s (QPM events 
with massless quarks have s=0). One can easily verify that for s > Q 2 the z momentum component 
of the hard subsystem is negative, pointing into the target hemisphere. 

3 In e + e~ data the tails of the distribution deviate from a Gaussian. In the ep data, no 
significant deviation from a Gaussian has been found yet. 
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the other hand, when coherence effects are neglected, and E q — ► oo a plateau 
F(£)=const. should form in the soft region, until the kinematic limit is reached 
at a scale £ = ln(l/x p ) ~ ln(E q /m). Coherence suppresses soft parton radiation, 
thus depleting the plateau at large £ - a "hump backed" plateau is formed?]. 

Analytical formulae in their various approximations for the spectrum are in gen- 
eral quite lengthy |5t], |197| , |198| . However, for large energies E q and concentrating 
on the main part of the spectrum F(£), without the tails, the peak and width of the 
Gaussian are predicted to evolve with E q as ||201 



£ peak = 0.5 • ln(£ g /A efr )+c 2 J\n(E q /A eS ) + « (6.9) 



^width = J[ln(E q /A eS )] 3 / 2 /(2ci). (6.10) 



The constants are c\ = y36n c /6, c 2 = BJb/(16n c ) with b = yn c — |n/ and 
B = -^(y^c + |^// n c)j they depend on the number of colours, n c , and the number 
of flavours, rif (here n/ = 3). k accounts for higher order corrections and is expected 
to be of order 1. A eff is an effective scale for the evolution. 

It is not a priori clear that the above MLLA+LPHD calculations for showering 
quarks can be applied to ep data, where BGF events play a role at small Q. A good 
combined fit to the HI data is obtained though, yielding A e ff = 0.21 ± 0.02 GeV 
and k = —0.43 ± 0.06, in agreement with an analysis of combined e + e~ data, for 
which A eff = 0.21 ± 0.02 GeV and k = -0.32 ± 0.06 was obtained p01|] . 



It should be noted that the existence of a hump backed plateau, fig. |67X], does 
not prove colour coherence effects. At finite energy, a hump backed, approximate 
Gaussian shape for the hadron spectrum can also be obtained in Monte Carlo simu- 



lations without coherence and with an independent fragmentation model |197| . One 



may argue that for an incoherent shower the turn over of the plateau happens at 
a scale x p ~ m/E q . This point would define £ pea k, hence d^ pea ^/d\n.E q = 1, in 
contradiction to the slower growth of the data and of the MLLA prediction. How- 
ever, in an incoherent shower model with string fragmentation, also a slower growth 



close to the data can be obtained |191|| . Furthermore, in ||197|| it is demonstrated 
how string fragmentation can turn a flat plateau for partons into a hump backed 
one for hadrons. String fragmentation appears to compensate for negligence in the 
perturbative evolution. On the other hand, the growth of the average charged mul- 
tiplicity with Q is less well described by string fragmentation with an incoherent 



parton shower 191 



4 A flat plateau dn/d£=const. implies also a flat rapidity plateau, dn/dj/=const., since dy/dp z = 
1/E. A hump backed distribution dn/d£ therefore results in a depict ion f or y — » 0. A dip in the 
rapidity distribution has been measured also in e + e~ annihilation ]199| , [200| . The observation 
of a dip for 2-jet like events in e + e~ annihilation provided evidence for soft and collinear gluon 
radiation 19£]. 
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Invariant energy spectra and the soft limit 



Based upon MLLA+LPHD, the spectra from fragmenting quarks have been calcu- 
lated also in the Lorentz invariant form 



1^ (6.11) 
N d 3 p 



to study the small momentum limit [2U2|. The calculation agrees well with e + e data. 



One salient feature is that for p — > the calculated spectra do not depend on 
the energy of the fragmenting parton. The invariant spectra from HI ||190| are 



shown in fig. |6.4J. The particle momentum p is measured in the drift chamber. 



Its energy is set to E = y Qo + \p\ 2 > following the presciption in [202], where the 
"particle mass" is given by the mass cut-off of the parton shower, taken to be 
Qo = 0.27 GeV (« A QCD ). 

Comparison of the MLLA+LPHD prediction for fragmenting quarks, which has 
been shown to agree with e + e~ data [|202|| , with the ep data is hampered by the 



presence of BGF events. At low Q, the prediction is far off the HI data [[L90|| , see 



fig. |6T4| a. That effect can be reduced by a special anti-BGF event selection ||190|| , 
which is based on the calorimetric energy seen in the Breit current hemisphere. That 
selection reduces the discrepancy at low Q. At large Q, where the discrepancy is 
smaller, there is little change after that cut, see fig. |6.4j b. This study shows that 
in the ep Breit frame spectra there are fewer fast particles than in e + e~ events at 
the same energy. The difference is presumably due to BGF events that are absent 
in e + e~. For larger Q, where the BGF contribution decreases, also the difference 
decreases. 

The prediction of an energy (here: E q = Q/2) independent soft spectrum is 
confirmed - the measured spectra converge for \p\ — > 0, independent of Q (fig. |S~4]c). 
It can also be noted that a calculation with fixed a s , instead of a running a s , is clearly 
untenable (figs. |6.4| a,b). A meaningful extraction of a s is however not possible yet, 
because the calculations are available only in LO. 



In the same HI analysis ||190 |, the KNO mulitplicity distribution was studied in 



the Breit current hemisphere. KNO scaling was found to be broken; at low Q the dis- 
tributions broaden, and deviate the most from the approximately scaling e + e~ data. 
Again, BGF events and the limited rapidity range can be held responsible for this 
observation. It would be interesting to test also KNO-G scaling ||203|| , which is a 



reformulation of the KNO form, valid for discrete distributions (like multiplicities) 
at finite energies. 



Scaling behaviour 

The scaling behaviour of the x p spectra in the Breit current hemisphere are best 
studied by displaying the Q 2 dependence of the normalized hadron cross section 
— da/dxp for fixed x p , where cr to t is the event cross section. With the large 
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Figure 6.4: Invariant charged particle spectra in the Breit current hemisphere 19C]. In a) and 



b) the data at small and large Q are compared to the MLLA+LPHD prediction [ 202 1 with running 
a s (solid line), and with constant a s (dashed line). The open and full points are before and after the 
anti-BGF selection. In c) the specta are shown for different Q 2 ranges. For better visibility, they 
are offset by an incremental spacing of 0.5 GeV . The full lines are the MLLA/LPHD expectations. 



range of Q 2 accessible at HERA, the scaling behaviour can be studied in a single 



experiment. Ultimately, a s could be extracted from such measurements [ 153|| . 

The ZEUS data P2fl , covering 10 < Q 2 < 1280 GeV 2 are shown in fig. The 
spectra were found to be almost independent of x for fixed Q 2 ||189| , |192| . Approx- 



imate scale independence is observed for x p ~ 0.3: the cross section — da/dx 



O"tot 

does not depend on Q 2 . For smaller x p , the cross section increases strongly with 
Q 2 (compare also fig. |6.1| a). These effects are well described by the ARIADNE 4.08 
p6| event generator; in LEPTO 6.5 |^TJ the scale dependencies are too strong, see 
fig. |6.5| . Possibly the strong coupling parameter for the final state parton shower in 
LEPTO can be adjusted to describe the data. 

In fig. |6.6| new preliminary ZEUS data (1995) ||193|| are shown together with 
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Figure 6.5: Charged particle cross section as a function of Q 2 for different x p bins. The data 
[ 192| are compared to a) ARIADNE 4.08 f|] and LEPTO 6.5 @; and b) to the NLO calculation 



with CYCLOPS [204]. 



lower statistics HI data (1994) fl90| for Q 2 > 100 GeV 2 . They agree very well. 
The ep data agree also very well with data from e + e~ annihilation (corrected for 
K° and A decay products, and divided by 2), supporting the notion of universal 
quark fragmentation in this kinematic region. The measured cross sections at small 
x p exhibit clear scaling violations with a positive slope. It is not clear however that 
this effect is due to the expected scale dependence of the fragmentation function in 
perturbative QCD. These expectations will be discussed now. 

The spectra have been calculated [|209|| according to eq. |5l| in NLO pertur- 
bative QCD (program CYCLOPS |[204| ) by folding parton density functions (here 
MRSA' [ plO|| ) with the NLO matrix element for parton production and a NLO 
fragmentation function | |211[ | for parton fragmentation into hadrons. All of these 
ingredients are scale dependent. Aq CD = 0.23 GeV was used. 

The calculation agrees well with the data in the kinematic region shown in fig. |6.6| . 
x p > 0.1 and Q 2 > 100 GeV 2 . In this regime, the expected scaling violations are 
small, and can hardly be established with the present precision of the data. With 
higher precision data, and with a larger lever arm both at higher and smaller Q 2 , 
the predicted scaling violations will be measurable. It will be very interesting to 
study to what extent this effect is due to the scale dependence of the fragmentation 
function, the matrix element, or the parton densities. 

The fragmentation function picture has its limitations [ 209 |. The fragments must 
not be too far separated in rapidity from the parent parton, assumed to fragment 
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Figure 6.6: Charged particle cross section as a function of Q 2 for different x p bins. HI data [|190| 



and preliminary ZEUS data 193 are compared to e + e data [205, 20£, 207, 20S], and to a NLO 



calculation with the program CYCLOPS [153, 204 



independently. It is estimated that the hadron rapidity in the Breit system should 
be above some minimal value y', with The condition y > y', translates into 

a lower bound on x p , 

2m T 1 



Q 



1 — tanh y D 



(6.12) 



with the transverse mass m-r, typically of O(0.5 GeV ). Furthermore, hadron masses 
are not taken into account with fragmentation functions. The calculation thus must 
break down when x p • Q/2 approaches m^. Indeed, the NLO calculation would by 
far overestimate the data in a region where Q 2 < 100 GeV 2 and x p « 0.1 p09| , 
see fig. |B~5| . We note however deviations also at Q 2 w 15 GeV 2 and x p > 0.5, where 



the NLO calculation should be reliable according to eq. |6.12. 
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6.2 Event Shapes 



Definition of event shape variables 

Information on the shape of events, - pencil-like, spherical, cigar-like, planar etc. - 
is conveniently given by simple functions of the hadron momenta pi, F = F{pj). 
Measurements of such event shape variables provide information about both per- 
turbative and non-perturbative aspects of QCD. Due to the large kinematic range 
covered, HERA is particularly well suited to study their scale dependence, for exam- 
ple on energy, and by that means disentangle the two contributions. Furthermore, 
provided the non-perturbative part can be estimated, the strong coupling constant 
a s can be extracted. 

First measurements at HERA of the event shape variables thrust T, jet mass 



parameter pc, and jet broadening parameter Be are reported by HI |[212| . The 
measurements are restricted to the Breit frame current hemisphere in order to avoid 
complications from the not so well understood region towards the proton remnant. 
The energy of the scattered quark in the Breit frame (in the QPM), Q/2, provides 
the scale for the process. Q values between 7 and 100 GeV are covered. 
The variables are defined by 

T:=max ^7^ Bc:= ^e7W Pc:= 1? = — • (6 ' 13) 

For the thrust calculation, the unit vector Ht, defining the direction of a thrust axis, 
is variedF] in order to maximize T, the normalized longitudinal momentum sum of 
all particles i. T = 1 for a collinear event configuration, and T = 1/2 for a spherical 
configuration. Thrust measurements are usually expressed in 1 — T, so that low 
values correspond to pencil-like events, as for the other event shape variables. 



In the HI measurement calorimeter clusters are used in the sums eq. |6.13| , which 
do not have a one-to-one correspondence to incident particles. However, these event 
shape variables are infrared safe quantities (see section |2.6| ), they do not change 
when an object in the sum is split: T{pi, ...pi, ...p n ) = F{vii ■■■ Z V%-, (1 ~ z )Pi-> ■■■Pn)- 

Thrust measurements 

Here we shall concentrate on the thrust measurements, see fig. |6.7|; the conclusions 



from the other shape variable measurements are similar. The measured mean value 
(1 — T) decreases with increasing Q, the events become more collimated. The thrust 
distribution as well as the mean thrust and their energy dependences are well de- 
scribed by the LEPTO model (version 6.1) which incorporates a perturbative 
phase with the LO QCD matrix elements and leading log parton showers, and the 
Lund string model for hadronization. 

5 H1 has also measured the thrust for hxed thrust axis, defined by the virtual photon. 



95 



Qualitatively, data from e + e annihilation [ j213| , |207| , |205| , |208| , |206| | follow the 



same trend as the ep data, but are systematically higher at small energies, Q < 
15 GeV. One would actually not expect exactly the same distributions, partially 
because in e + e~ the full event with two hemispheres is usedf]. Furthermore, the 
flavour compositions differ, and in ep there are 0(a s ) processes contributing, namely 
initial state QCD radiation and boson-gluon fusion, which are absent in e + e~. In 
fact, the other shape variables (pc, Be) show bigger differences than thrust between 
e + e~ and ep data. 

Power corrections 



It has been argued [101] that for any infrared safe variable the mean value (and 



higher moments) can be written as 

(F)= (F) pert - + (F) pow \ (6.14) 

The perturbative part (_F) pert ' is calculable in fixed order perturbation theory, de- 
pending only on a s and the choice of renormalization scaleQ The only energy de- 
pendence comes from the running a s , since there is no other intrinsic scale in pertur- 
bative QCD than Aqcd • QCD radiation patterns would be identical at all energies, 
if a s were constant. If the shape variable is invariant against momentum scaling 
with a constant c (like T), Tipi) = T{c ■ Pi), the distribution of the event shape 
variable depends on Q only through the running a s : 

d diV d dN ( d dN\ da s 

+ b— 1TB --ST, (6-15) 



dQdF OQdF \da s dF J dQ 

because the first term on the right hand side vanishes. In ep reactions a Q depen- 
dence may arise however via the Q dependent parton densities. 

Higher orders and hadronization effects are collected in the "power correction 
term" (F) pow ', which according to ||101|| scales approximately oc 1/Q (in general 
oc (1/QY with the power p). A plausibility argument for such a term is given 
below. In contrast to the perturbative term, this term does depend explicitly on 
the energy, because hadronization introduces a scale, given by hadron masses. The 
higher order perturbative effects are calculated with running a s down to a scale /if, 
with Aqcd "C Hi to avoid the divergence of ot s for Q — > Aqcd , and fii <S Q to 
maximize the perturbative evolution. A conventional choice is \ii = 2 GeV . Below 
the scale /i/, the calculation uses a constant effective coupling o7o, a free parameter, 
which has to be determined by experiment, and which depends on the choice of /i/. 

The NLO expressions with Q as renormalization scale are [|101|| 

(F) pcrt - = Cl a s (Q 2 ) + c 2 a 2 s (Q 2 ) (6.16) 

6 Note that the thrust for a system of two identical back to back jets equals the thrust of a single 

jet. 

7 For the purpose of extracting a meaningful, process independent a s , the calculation has to be 
done in fixed order perturbation theory to at least 2nd order (NLO). 
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Figure 6.7: Top: thrust distributions pi2fl . The differential thrust distribution (l/iV)diV/d(l - 
T), compared to LEPTO 6.1 (top left) and a NLO calculation (top right). The spectra for (Q) = 
8.3 — 68 GeV are multiplied by factors of 10™, n — 0, ...4. Bottom: the mean (1 — T) as a function 
of Q. The HI data are compared to (bottom left) e + e~ data for the full event (two hemispheres) 
at y/s — Q, and the LEPTO result for hadrons and partons; and (bottom right) to the NLO 
calculation, and the fit result including the power correction. 
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pow. 



16 Hi 

a F ~ 

3tt Q 



27T V jJLl (3 



1 



(6.17) 



where (3 = 11 - (2/3)n/, # = 67/6 - tt 2 - (5/9)n/ and n/ = 5 flavours^. The 
constant ap is predicted by the theory and depends on the event shape variable F 
under consideration (see tab. |6.1|). In the HI analysis, the constants C\ and c 2 are 



calculated in the MS scheme with the program DISENT [|7|]. For e.g. (1 — T) they 
are c x = 0.384 ± 0.033 and c 2 = 0.57 ± 0.21. 

With the ansatz eq. |6.14| a good fit to the HI (1 —T) data (and (pc), (Be) 
as well) is obtained (see fig. |6.7|) with a s and a as the only free parameters. The 
perturbative contribution decreases with energy Q due to the decreasing a s . The 
power correction is substantial, but decreases from « 50% at Q = 8 GeV to ~ 20% 
at Q = 70 GeV. The prediction from LEPTO for partons is higher and closer to 
the data than the NLO calculation (see fig. |6.7| ), because the LEPTO parton shower 
includes QCD radiation in principle to all orders in the leading log approximation. 
In the ansatz |6.14| , the higher orders are part of the power correction term. Whereas 
the differential thrust distribution dN/dT is well described by the LEPTO model, 
the NLO calculation (covering states with up to 3 final state partons) is insufficient, 
producing too narrow events at small Q, but approaches the data for larger Q. This 
is not surprising, as the large NLO correction to the LO in eq. |6.16| (NLO/LO = 
«s cz/ci) suggests important higher order corrections, which should decrease with 
increasing energy due to the running a s . 

From the thrust fit one obtains for the two free parameters ao(A*j = 4 GeV 2 ) = 



0.497 ± 0.005lffi, 



and a Ami 



0.123 ± 0.002l8;885 



for the strong coupling ex- 



trapolated to the Z massQ. Similar values are obtained from the other event shapes, 
see tab. |6.1| 1U | . Assuming the power corrections to be universal (same «o for all 
shape variables), a combined fit to the thrust and jet mass data yields a s 



o.ii8±o.ooil8:886 



212 ; one has to admit however that the small combined statistical 



error disregards the fact that the data samples that enter the different distributions 
are the same, not at all statistically independent. To be on the safe side, one should 
quote 0.003 as statistical error instead. The obtained value agrees with the world 
WM with an error of similar size as in other determinations (see tab. |6TT 



average 



and 215 



It is quite remarkable that the ansatz |6.17| gives consistent values for a s , and 
within 20% the same constant Uq for all investigated shape variables. Further- 



8 It is not clear what should be used for nj in this situation. The lowest values ofQ/2 are above 
the charm threshold, the larger ones are above the bottom threshold. On the other hand, for the 
parton shower with many soft emiss ion s onl y li ght flavours are expected to be active. In fact, in 
the multiplicity analyses in sections [5.2| and |6.l| n / = 3 was used. 

9 Note that eq. 6.17 is not any longer a pure power law oc 1/Q, but is modified by corrections 
up to 0(a 2 s ). A pure po wer law fit of (I — t^p ow = 2\/Q would give a worse % 2 than the "QCD 
improved" expres sion eq. 6.17 [214]. 

10 According to [ 101 j , the power correction term for the jet broadening parameter Be should be 
oc (1/Q)ln(5. From a comparis on o f the HI data with the DISENT NLO calculation the extra 
lnQ term cannot be supported ]212[. 
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more, similar values are being obtained in shape analyses of e + e data, see table 



Q| . These encouraging results support the concept of universal power corrections, 
which deserves further investigation of its origin in QCD, and of its limitations. 
Besides new insights into hadronization, a better understanding would also improve 
the a s measurements. Recently also "power corrections" for the event shape distri- 
butions, not just the mean values, have become available [2TB ]. It would be very 
interesting to compare these calculations to the data. 



Observable 


ap 


a (ti 2 = 4 GeV 2 ) 


a s (m|) 


HI DIS 

(l-T) 
(Be) 
( Pc ) = (M 2 /Q 2 ) 


1 

2 
1/2 


0.497 ± 0.005^36 
0.408 ± 0.006±° ° lt 
0.519 ± 0.009t°;o2o 


0.123 ±0.002±° °^ 
0.119 ±0.003t£$2 
0.130 ±0.003±g ^5 


e + e data 

(l-T) 
(M 2 H /Q 2 ) 


1 
1 


0.519 ±0.009±g^ 
0.431 ± 0.020±°;°^ 


0.123 iO.OOlto'SoI 
0.115 ±0.002t° ° l 



Table 6.1: Results of the event shape analyses [212|. Following 101], different values for ap are 
used for the jet mass parameter in ep and e + e~. 



Naive model for power corrections 

The power behaviour oc 1/Q can be made plausible [[75[] with the simple tube or 
longitudinal phase space model ||217| | for hadronization. In this model a colour 



connected pair of partons with total energy Q produces two back to back jets of 
hadrons with energy Ej = Q/2 each, where the hadrons are distributed uniformly 
in rapidity y and according to a probability density function p(pr) in transverse 
momentum. The resulting hadron distribution function is then $(?/,pr) = n p(pT), 
where Uq gives the number of hadrons per unit rapidity. The energy and longitudinal 
momentum of one jet are given by 

Ej = / n p(p T )Edp T dy = Asinhy max (6.18) 

Jo Jo 

fUraax rco 

Pjz= / n p(p T )p z dp T dy = A(coshy max - 1) w Ej - A, (6.19) 



io Jo 

where A is related to the average hadron pr by 



A 



POO 

■= / n p(p T )p T dp T = n (p T ). (6.20) 
Jo 



Used were E = ttlt coshy, p z = itlt sinhy, ttlt ~ Pt for light hadrons, and the quick 
convergence of cosh?/ towards sinhy. For a pencil like jet pj Z = Ej would have been 
expected; hadronization yields a negative correction of relative size X/Ej, a simple 
power of the energy. 
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In this model we can calculate hadronization corrections to the otherwise van- 
ishing 1 — Tj and jet mass Mj of the jet. According to the definitions eq. |6.13 , 



M 2 



T* = 1 - ^ 
3 Ei 



En ~ A 



E; 



E 2 -p 2 

3 F 3 Z 



2A 

~Q 



2Ej\ 



A_ 

Q' 



(6.21) 



(6.22) 



Q 2 Q 2 Q 2 Q 2 

The corrections are proportional to 1/Q, and twice as large for thrust as for the 
scaled jet mass squared, in accord with the constants ap for ep in tab. BTT. The 
different values for ap in e + e~derived in ||101|| appear somewhat counterintuitive. 
One has to take into account however that the structure of DIS and e + e~ events is 
more complicated than just two separating colour charges. 

We can estimate the order of magnitude of the parameter A = noipp) from the 
data. From the seagull plot at low energies, where QCD effects are small, and for 
xp — > where intrinsic kp effects vanish, we estimate (pp) ~ 0.3 GeV , see fig. |5.14 . 
From measured hadron rapidity distributions one gets for the hadron density per 
unit rapidity n w 2.0 ||, compare also fig. That yields A ~ 0.6 GeV. In fact, a 
good description of e + e~ event shape variables is obtained with power corrections 
parameterized with A = 0.5 GeV ||75| . The fitted a from the thrust analysis yields 
for the 1/Q term in eq. |6.17| a coefficient of similar order, 2A ~ 1.7. 
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Chapter 7 
Jets 



7.1 Introduction 



Amongst the spray of particles emerging from a high energy reaction the human 
eye can recognize collimated subsystems of hadrons, so-called jets (fig. |L2[) . Jets 
are fascinating objects; they allow to view high energy quarks and gluons which are 
not observable as free particles due to confinement. As such they can be considered 
as the parton images imprinted on the hadronic final state. At short distances, 
partons can be treated as "asymptotically free"; their interactions can be calcu- 
lated perturbatively. When the partons emerge from the confinement volume, they 
fragment into observable hadrons. At high energies, the hadrons will be collimated 
around the original parton direction. These hadron jets still carry information on 
the underlying partonic interactions. 

Jets in DIS at HERA result from the scattered quark, and from additional QCD 
radiation either in the initial or the final state. Clean jet structures can develop due 
to the large available phase space (W up to 300 GeV). The main experimental chal- 
lenge in jet studies is twofold: 1) to measure jet quantities like energy-momentum, 
charge etc., and 2) to relate them to the corresponding parton quantities. A good 
correlation between partonic and hadronic jet quantities is vital to allow perturba- 
tive QCD to be tested. 

Jets at HERA have been studied with different intentions. The early measure- 
ments established clear jet structures in DIS |P18| , |219|| . Fragmentation properties 

E2TL Mea- 



of quark and gluon jets have been studied by measuring jet shapes 
surements of jet rates allow to study the partonic QCD mechanisms and to measure 
the parameters by which they are determined, like the strong coupling strength 
222[ |223|| , or the density of gluons in the proton ||224|| . It turns out that parts of 



a 



the phase space for jet production, in particular at small x and Q 2 , is not yet well 
understood [f225|, E2!| [227], [22J. The dedicated search for BFKL "footprints" with 
forward jets is covered in a chapter on low x physics (section 
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7.2 Jet Algorithms 



For quantitative measurements jets need to be defined: a jet is the output of a jet 
algorithm. A jet algorithm prescribes how to combine objects (partons, hadrons, 
energy deposits, ...) close in phase space to jets. It hinges on a resolution parameter 
below which jets cannot be resolved. The jet algorithm needs to be implemented 
not only experimentally, but also in the theoretical calculation for meaningful pre- 
dictions. There exist several jet algorithms for different applications. 

Clustering algorithms, like the JADE ||229|| or the k? algorithm 1230 1, combine 



affinal state particles into few jets, based upon a measure of distance. In the JADE 
algorithm, all particles i, j with invariant mass 



m% = 2E l E j ( 1 - cos % ) < y cut ■ W 2 (7.1] 



are merged to form new objects. (Experimentally, at HERA a pseudoparticle carry- 
ing energy and momentum of the unobserved proton remnant is introduced.) Here 
y cn t is the resolution parameter, and W serves as reference scale. The final set of 
jets is obtained, when no further merging is possible. 
In the hp algorithm two objects are merged, if 

2 • min(£?, E])(l - cos%) < k 2 cut . (7.2) 

Particles can be merged also with the remnant. At HERA either fc cut = Q is chosen 
||231|| , or k cut is set to a fixed value, for example k cnt = 3 GeV ||225|| . For small 
angles, the merging requirement becomes kj- < k cnt , where k? is the transverse 
momentum of the less energetic object with respect to the other one (hence the 
name &t algorithm). 

For the cone algorithm | |232| | a distance measure in r] (pseudorapidity) - <fi (az- 
imuth) is introduced: AR = J (Ar/) 2 + (A0) 2 . All particles inside a cone with radius 
R (at HERA usually R = 1) are combined to form a jet, if the resulting jet with 
respect to the proton axis is above a certain cut-off (usually a few GeV). The cone 
axis is chosen such that the Et inside the cone is maximized. 

The result of a jet algorithm depends on the way masses are treated when com- 
bining objects (recombination scheme). Measured jets consisting of more than one 
particle always have a finite invariant mass, while most jets in NLO calculations are 
massless (they consist of just one massless parton). For this reason it is often found 
that the final result (e.g. a s ) depends on the chosen recombination scheme. A lot 
of effort is being spent on minimizing such effects. For the cone algorithm it also 
matters in which order objects are combined. 



7.3 Jet Shapes 

ZEUS has measured the internal structure of jets in DIS and has compared them 



to jets in other reactions |220|. One expects gluon jets to be broader than quark 
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jets due to the larger colour charge of the gluon. In the laboratory frame, jets with 
Et > 14 GeV and rj e [— 1,2] were selected with the cone algorithm (R = 1) for 
events with Q 2 > 100 GeV 2 . In most cases the current quark jet recoiling from 
the scattered electron will be selected. 

One measures the average fraction ip(r) of jet transverse energy Et{t) inside a 
cone of radius r around the jet function r, defined by 

iVjets is the number of jets in the event sample. ip(R) = 1 by definition. The jet 
shapes are corrected for detector effects to the hadron level. 

It is found that the jet shapes narrow as Et increases, as expected kinemati- 
cally. (A jet will look narrower with increasing boost.) They show no significant 
7] dependence. (In a preliminary analysis of 2+1 jet events in the Breit frame, HI 
found that these jets become broader towards the remnant |[221|| .) The gross shape 
features are reproduced by the models investigated, PYTHIA 5.7 ||233|| , LEPTO |8l[ 



and ARIADNE |86] , but none of them provides an accurate description of the data 
over the entire phase space. 

On average, DIS jets are seen to be narrower than jets from photoproduction. 
This can possibly be explained by the larger fraction of gluon jets in photopro- 
duction, mostly from resolved processes, where a quark from the photon scatters 
with a gluon from the proton, qg —>■ qg. (However, multiple interactions between 
the proton and photon remnants could also lead to a a jet broadening in photo- 
production ||234||). A similar observation is made when comparing DIS jets to jets 



from e + e~ annihilation [|235|| , and to jets from pp interactions ||236|| with similar 
Et, see fig. [77I|. DIS jets and e + e~ jets are predominantly quark jets and have the 
same universal shape within errors. They are narrower than the jets measured in 
pp interactions which contain a significant fraction of gluon jets. 



7.4 The Strong Coupling a s 

Analysis method 

Considering the LO graphs for dijet production (fig. 7.2), the QCDC Compton 



(QCDC) and boson-gluon fusion (BGF) graphs, it is clear that the rate of events 
with 2+1 jets (the +1 refers to the remnant jet which escapes largely unobserved 
down the beam pipe) depends on the strong coupling constant a s , and on the parton 
densities in the proton. Inspired by the great potential to measure the behaviour of 
a s over the large span of Q 2 accessible at HERA with high statistical precision, jet 
studies have focussed on the determination of a s right from the beginning of HERA 
data analysis. For this purpose one assumes parton densities that are constrained 
by other data. 
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Figure 7.1: The average jet shapes ip(r) for DIS charged (CC, ep — > vX via W exchange) and 



neutral current (NC, ep — > eX via 7, Z exchange) events ]220[, compared to e + e jets [ 235 1 and 





DIS variables: 
Q 2 = -q 2 
x = Q 2 /(2P-q) 
W 2 = Q 2 (l-x)/x 

s = (ji + J 2 ) 2 
MMMP £ = X (Q 2 + s) IQ 2 
z 3 = {P. 3 )/{P.q) 



Figure 7.2: Feynman diagrams for the production of 2 + 1 jet events to first order a s in ep - 
collisions, q and P denote the four-momenta of the photon and proton, ji and ji are the four- 
momenta of the jets associated to the hard subprocess with invariant mass squared s. x g / p and x g / p 
denote the proton fractional momenta carried by the gluon or quark entering the hard subprocess. 



The kinematics of the hard subprocess in 2+1 jet evenst is depicted in fig. |7T2 . 
One defines the invariant 




s 



(7.4) 



where s := (j\ + J2) 2 is the invariant mass squared of the 2-jet system. In LO £ can 
be identified with the momentum fraction of the proton carried by the parton that 
enters the hard scattering, x q / p or x g / p for a quark or gluon. 
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At present the strategy is to measure a s with dijet data at relatively large Q 2 , 
where the QCDC process dominates with well known quark densities at large £. 
Dijet data at small Q 2 , where the BGF process is the dominant one, are used to 
determine the less constrained gluon density at small £, assuming a value for a s , see 



section |7.5| . Best use of the information in the data will in the future be made from 
a combined determination of a s and the gluon density. 

In the measurements the ratio of the number of 2+1 jet events iV 2 +i to the 
number of all events N tot , 

R W := ^ (7.5) 

-'not 

is determined as a function of Q 2 . (HI defines N to t '■= ATi+i + ^2+1 •) In the 
a s analyses the JADE algorithm is applied in the laboratory frame with resolution 
parameter y cnt =0.02. Another approach uses the dijet rate differential in the JADE 



resolution parameter y cnt , but no a s value has been quoted yet ||237|| . R2+1 is being 
corrected for a) detector effects and b) hadronization effects to the parton levelQ, 
and is then compared to the QCD prediction in NLO: 

R 2+1 = A ■ a s (Q 2 ) + B ■ (a s (Q 2 )) 2 + 0(a 3 s ). (7.6) 

A and B are calculable coefficients that depend on the kinematic conditions, the 
QCD renormalization and factorization scales, and on the parton densities. a s is 
determined from a comparison of the measured jet rate, corrected to the parton 
level, with the NLO prediction of eq. |7T6. 



The analyses cuts are designed such that a) the observed jet rates are well de- 
scribed by the Monte Carlo used for corrections; b) the correction factors are reason- 
ably small; c) the NLO prediction for partonic jets resembles the partonic jets found 
in the Monte Carlo generator. Jets from initial state parton showers complicate 
matters; to a large extent they are not covered by the NLO calculation. 

Data analysis 

The corrections for detector and hadronization effects rely mostly on the MEPS 
model (program LEPTO) |8l|], where the hard scattering at the photon vertex is 
calculated with the LO matrix element; higher orders are taken into account with 
leading log DGLAP parton showers. Hadronization is performed with the Lund 
string model |J3|. Other programs (ARIADNE Jg§, HERWIG @) are used as 



cross checks. Fig. |7T3| a shows R2+1 ||231|| as a function of Q 2 in comparison to the 
models. The data are corrected for detector effects to the hadron level. LEPTO 
describes this rate reasonably well. It is not understood why the other models do 
not describe this data at relatively large Q 2 , where little theoretical ambiguity is 
expected. 

The QCD calculations in fixed order perturbation theory are performed with the 
programs PRO JET |!38|, DISJET |239j, MEPJET and DISENT J7§. PRO JET 



1 The parton level needs to be carefully defined, see section 4.5 
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Figure 7.3: a) The jet rate -R2+1 as a function of Q 2 for hadron jets (corrected for detector 
effects), reconstructed with the JADE algorithm with y cut — 0.02. The preliminary HI data [ (23 1| 
are compared to the models ARIADNE 4.08 @, LEPTO 6.5 §l| and HERWIG 5.9 §§. b) 

223l 



The strong coupling constant a s as a function of scale Q 2 . Shown are the HERA data 1 231 



together with the world average for a s {m 2 z ) |l4| ]. The lines give the expectation for the running 
a s with the QCD scales (for n f = 5 flavours) Ai^!_ = 0.2 ± 0.1 GeV. 

1V1 o 



and DISJET are no longer in use, because parts of the phase space at small Q 2 and 
large W were not fully included^]. 

In the following two paragraphs, typical features encountered in the HERA 
a s analyses are reported. One finds 3 — 20% 2+1 jet events in the total event sample 

EH 



327|| , depending strongly on y CVLt and also on the event and jet selection. The 
dependence of R2+1 on y CVLt for Q 2 > 100 GeV 2 is well described by Monte Carlo 
(LEPTO 6.1) |n| and by NLO calculations (PRO JET) |223 



With increasing y cut , 
fewer structures are resolved as jets. According to the Monte Carlo (LEPTO 6.1), 
around 75% of matrix element 2+1 jet events are due to BGF for Q 2 < 100 GeV 2 
(50% for Q 2 > 100 GeV 2 ), the rest is due to QCDC [|T|. In LEPTO, around 
50% of the 2+1 jet events originate not from the LO matrix element, but from the 
parton shower, that is mostly initial state radiation ||222| , |241|| . The parton shower 
contribution is reduced at high Q 2 [242| . 

The NLO correction to the LO prediction is roughly 20%, depending on the 
analysis details, in particular on y cut ||223|| . Migrations between jet classes due to 
hadronization effects can be large: ~ 60% of 2+1 parton jet events and ~ 90% of 
partonic 1+1 events are correctly classified after hadronization [[21911 . The result- 
ing corrections for hadronization are sizeable, around 30% |[231 |. Similar numbers 



2 The published ZEUS analysis [223 1 presented here was based on PROJET. It has been verified 
thou gh th at due to the cuts applied the effect on the a s extraction is negligible (less than 2 per 



mil) m 
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are obtained for migrations due to detector effects. In general, smearing effects 
become less severe for larger Q 2 . Therefore the a s analysis apply lower Q 2 cuts, 
120 GeV 2 for ZEUS |2| and 40 GeV 2 for HI f23T] . 



Two Lorentz invariants can be introduced to describe the jet kinematics, for 
example s = (ji + j2) 2 ; the invariant mass squared of the two jets, and 

z< ■,:= « -(1 -cos0A (7.7) 

3 p_ q 2 V J) \ ) 

The index j denotes the jet from a particular parton. Zj G [0, 1] and Zj x = 1—Zj 2 . The 
approximate relation with the scattering angle 6j in the parton-photon scattering 
frame (the rest frame of the two jet system, fig. |7.4j a) holds for small jet masses and 
small Q 2 . In the experiment one defines z := mm(zj 1 , Zj 2 ) with z G [0,0.5]. z can 
be reconstructed from the measured jet angles 9j and energies Ej in the laboratory 
frame, 

Ej(l- cos 6^ 



z « min — ^ . - — — (7.1 
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Figure 7.4: a) Two jet production in the centre of mass system of the incoming parton and 
the virtual photon. The scattering angle 9j of jet j is defined with respect to the incoming parton 
direction, b) The z distribution of one of the two non-remnant jets from the ZEUS analysis [ p4Cfl . 
The ZEUS data (223) are compared to the NLO calculations from PROJET p38| ] (dashed line) 
and from MEPJET (79) (full line). The histograms show the NLO results with the same binning 
as the data. 



The LO matrix element for the QCDC process diverges cx 
with the familiar soft (£ — > x) and collinear (zj 



-> 0, 1) divergencies. 



BGF oc -^-^ 

It turns out that close to the collinear divergency (z < 0.1), more 2+1 jets are 
observed than expected from either the NLO calculations (see fig. |7.4j ) or the LEPTO 
simulation ||227| , [22311 . In this part of the phase space the NLO corrections are largest, 
and according to the Monte Carlo, parton shower effects are most severe (apart from 
the LO BGF and QCDC events, QPM events with initial state parton showers can 
contribute up to 50% to the total 2+1 jet rate). The region z < 0.1, corresponding 
to small angles 9 jet in the lab system with respect to the proton remnant direction, 
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is therefore cut out for the determination of a s ||223|, |23 1|| . HI applies an additional 



cut 9 j et > 10°, and requires 9- ]e t < 145° in order to exclude the backward detector 
region which was not instrumented with hadronic calorimetry. 



Results 

From a comparison of the jet rate -R2+1, corrected to the parton level, with the NLO 
is determined function of Q 2 



a. 



expression of |7T6 

three Q 2 bins from 120 to 3600 GeV 
to 4000 GeV 2 , see fig.O 



ZEUS has measured a s in 
1 and HI (preliminary) in four bins from 40 
The data are consistent with the decrease of a s with 
increasing scale according to the renormalization group equation. A fit to the data 



with the NLO expression for the running a s (Q ), eq. with a s (m z ) as free 

parameter yields the results in table 7.1 , These a s measurement from jets at HERA 



agree with the world average . 



ZEUS 


a s 


>!) = 


0.117 ±0.005 (stat) 


+0.004 
-0.005 


0y s exp) 


±0.007 (sys th ) 


HI prel. 


a s 


>!) = 


0.115 ±0.003 (stat) 


+0.008 
-0.011 


(sys) 


±0.006 (rec) 


world avg. 


a s 


>!) = 


0.118 ±0.003 








e + e~ 


a s 


>!) = 


0.122 ±0.006 









Table 7.1: a s measurements by HI (prel.) |23l| and ZEUS |223| from 2+1 jet events. ZEUS 
quotes experimental and theoretical systematic errors separately. HI finds that the use of alter- 
native recombination schemes would increase the measured a s value by the amount given (rec). 
For comparison, the world average from the PDG Jl4|] is quoted. The combined result from LEP 
and SLC hadronic final state analyses at the Z mass is also given ||215f . Sec [ 215 for a recent 
compilation of other a s measurements. 

For the ZEUS result, experimental and theoretical systematic errors are given 
separately. HI found that alternative recombination schemes than the ones used 
would increase a s by the amount quoted. The dominating systematic errors are the 
corrections for hadronization effects and the dependence on the choice of renormal- 
ization and factorization scales. The largest experimental systematic error is due 
to the uncertainty of the absolute calorimeter calibration, which is known to about 
±5%. 



It is estimated |[243|| that with a large amount of luminosity to be collected at 
HERA in the future, 250 pb" 1 , the errors on (^(m^) can be reduced to ±0.0013 
(stat.) and ±0.005 (syst.). Other algorithms than JADE may turn out to be ad- 
vantageous in the future |]231j] . In any case, a good understanding of perturbative 
QCD beyond LO and hadronization will be mandatory to achieve a good precision 
on a s . For example, the puzzle why two otherwise respectable QCD models fail to 
describe the jet rate at the hadron level (fig. f7T3|a) needs to be resolved. 
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7.5 The Gluon Density 



The gluon density at small fractional momentum x g / p {x g / p < 0.04) had been ex- 
tracted previously only indirectly from inclusive structure function measurements. 
A direct measurement of the gluon density is in principle possible via the measure- 
ment of 2+1 jet events from the boson-gluon fusion (BGF) process 7*5 — > qq. The 
2+1 jet cross section can in leading order QCD be described as the sum of quark 
initiated and gluon initiated processes, QCDC and BGF. Schematically, 

CX2+1 = a s {Q 2 ) [Aqcdc • q(x q/p , Q 2 ) + A BGF ■ g(x g/p , Q 2 )\ . (7.9) 

Here x q / p and the fractional proton momenta carried by the incoming 

quark and gluon. They are given in LO by £ = x(l + j4r) ~ pfi , see fig. [77|. A BGF 

and Aqcdc are coefficients that can be calculated perturbatively. The experimental 
challenge is to reconstruct £ by measuring s. That requires to separate in the 
hadronic final state the hard subprocess from initial state parton showers and the 
proton remnant. 

For the extraction of the gluon density g(x g / p , Q 2 ), a s (Q 2 ) is assumed to be 
known and taken from the PDG world average. The contributions from the com- 
peting QCD Compton process (QCDC) 7*g — > qg, which is calculable using known 
quark density functions q(x q / p ,Q 2 ), and from wronlgy classified QPM events have 
to be subtracted statistically from the observed 2+1 jet rate. 

Of particular interest is the gluon density at small x g / p , where no direct mea- 
surements exist. In order to access small values of £ (= x g / p in LO), small invari- 
ant jet-jet mass need to be resolved. The JADE jet algorithm resolves 2-jet 
events with mf^/W 2 > y cut . That defines the accessible phase space, and implies 
£ « s/W 2 > y cut . 

To reach lower values of £ than given by the canonical y cut = 0.02, the HI 
analysis [|224j| uses the cone jet algorithm (R = 1) applied in the hadronic CMS. For 



the jets pt > 3.5 GeV is required. Since for massless objects 

V% = sz(l - z), (7.10) 

values of > %Pt can be obtained. With the large W 2 > 4400 GeV 2 , much lower 
£ values can be probed than with the JADE algorithm. 

From DIS events at small Q 2 (12.5 < Q 2 < 80 GeV 2 ), 2+1 jet events are 
selected with 10° < 6- iCt < 150° in the laboratory frame. Their distance in (lab) 
pseudorapidity had to be smaller than |A?y| =2. In the hadronic CMS, that cut 
translates roughly to z > 0.1. These angular cuts remove events affected by parton 
showers. \/I > 10 GeV is required, s is reconstructed 1) as the invariant mass 
squared of all particles belonging to the two jets, and 2) from the jet directions 
(pseudorapidities r]j) measured in the hadronic CMS, 

§ = W 2 exp(-r ]l -r] 2 ). (7.11) 
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It is required that the two methods agree to |AVI| < 10 GeV. 

The selected 2+1 jet event sample covers 0.002 < £ < 0.2 and 0.0003 < x < 
0.0015. According to the Monte Carlo (LEPTO), BGF and QCDC events are se- 
lected in a ratio 3:1. By unfolding the observed £ distribution, the (LO) gluon 
distribution x g / p ■ g(x g / p ) at (Q 2 ) = 30 GeV 2 is obtained (fig. [T.5|) 



The data ex- 
tend the region previously covered by NMC |244| , x g / p > 0.04 down to x g / p = 0.0019. 
This direct measurement of the gluon density confirms the steep rise of the gluon 
density at small x which had been deduced from the structure function data (see 
fig. |7. 5|) . A preliminary ZEUS analysis using a similar method gives consistent re- 
sults [[241 1, within large errors. We note that sizeable NLO corrections for 2 + 1 jet 



production have been calculated [245 



Recent analyses focus on the direct determination of the gluon density in NLO 
46| . Using the measured partonic jet rate R2+1 (JADE algorithm) for Q 2 > 
40 GeV 2 from the HI a s analysis [|231|| as input (see section |7.4j ), the extracted 
gluon density for x g / p > 0.01 is consistent within errors with results from structure 
function analyses |[247|| . Due to the phase space restriction implied by the JADE 
algorithm, not very small probed, x g / p > 0.01. Smaller x g / p va\ues can be 

accessed by tagging gluon initiated BGF events with identified charm decays ||175 



see section 5.5 



Recent attempts to extract the NLO gluon density at small x g / p 
production at low x and Q 2 are facing a problem [|228| , |248 



from dijet 
There are about 20- 



30% more dijets in the data than expected from NLO calculations or from e.g. 
LEPTO with gluon densities as input that are compatible with the F 2 data. At face 
value, that would indicate a directly measured gluon density that is incompatible 
with the indirect extraction from the F 2 measurements. It could also be that there 
are effects contributing to jet production at small x and Q 2 that had not been taken 
properly into account. To clarify the situation, data on jet production are essential 
that do not contain any theoretical bias from the correction to the parton level. 
Such measurements are presented in the next section. 



7.6 Jet Rates at Low 



In a preliminary ZEUS analysis ||228|| that addresses the gluon density, the dijet cross 
section (cone algorithm in the laboratory frame with R — 1) has been determined 
and corrected for detector and hadronization effects to the parton level using LEPTO 
6.3. In the laboratory frame, r]j Ct < 2 was required to exclude the forward region. A 
cut prjet > 4 GeV was applied in the laboratory frame and in the CMS. 

The shapes of the measured cross sections as a function of x , Q 2 , PTjet and 
7] in the CMS, and £ = x(l + s/Q 2 ) are well described by LEPTO and by NLO 
calculations |79| . However, the absolute dijet cross section is ss 34% larger than 



predicted by NLO QCD (see fig. [7T^| ). Systematic effects cannot account for such a 
large discrepancy ||228|| . 
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Figure 7.5: The unfolded LO gluon distribution as a function of fractional gluon momentum 
Xg/ p at (Q 2 )=30 GeV 2 [ 224 1 . The direct measurement from 2+1 jet events is compared to the 
indirect LO determinations from structure function measurements at (Q 2 )=20 GeV 2 by HI [ [249| 
and ZEUS 1 250 , and to data at large x from NMC 244] obtained from J/tp production. The other 
data are evolved to the (Q 2 )=30 GeV 2 of the jet analysis. Also shown are the LO gluon density 
parametrizations from GRV pi and CTEQ p51|. 



The previous jet studies were aimed at the measurement of observables defined 
in QCD - the strong coupling constant and the gluon density. To that end regions 
of phase space where the data deviate from the theory - cast either into NLO calcu- 
lations or Monte Carlo generators - had been carefully avoided. These phase space 
regions were found at low x and low Q 2 , and in the forward region close to the rem- 
nant |[252| ) [228| , p48|| . The forward region at low x is discussed in a special chapter 
on low x physics, section iO. Here we shall discuss jet production at low Q 2 , with 
Q 2 ranging from the DIS regime down to photoproduction, < Q 2 < 100 GeV 2 . 
The object of these studies is a test of QCD and to explore the validity of QCD 
models for jet production. The data extend into a new phase space region in DIS, 
where the transverse momenta squared of the jets may exceed Q 2 , p% > Q 2 . In such 
a situation it may be possible to probe constituents of the virtual photon at a scale 
given by the of the jet ||253| . 

In a preliminary HI analysis [ |226|| , low Q 2 DIS events are selected with y > 0.05, 
where the scattered electron angle and energy are 156° < 9 e < 173° and E' e > 
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Figure 7.6: Differential dijet cross section (ZEUS prel. ||228|] ) as a function of £, corrected to 
the parton level (LEPTO 6.3). Shown are statistical errors only. The data are compared to QCD 
calculations in LO and NLO. 



11 GeV . The data sample covers 5 < Q 2 < 100 GeV 2 and 10" 4 < x < 10~ 2 . 
Exactly two jets are required, reconstructed with the cone algorithm (AR = 1) in 
the hadronic CMS and with transverse momenta pt > 5 GeV . The remnant is 
not counted. The pseudorapidity difference \Ar]\ between the two jets has to be less 
than 2. In the CMS of the scattering parton and the virtual photon, the jet angles 
with respect to the photon axis are given by 

0i,2 = 2arctanexp(±Ar?/2). (7.12) 

The requirement |A^| < 2 thus excludes very forward jets, because | cos6*| < 0.76, 
or z > 0.12. Measured is the dijet rate R2, the fraction of DIS events that fulfil 
the dijet requirements. The rate is corrected for detector effects to the hadron level 
(fig. |7.7| top). The dijet rate increases from 4 to 10% for Q 2 between 10 and 100 
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Figure 7.7: Top: the dijet rate Ri as a function of Q 2 [ 226 1 . The data are corrected to 
the hadron level and compared to the models ARIADNE 4.08 [§(J and RAPGAP 2.06 |§| with 
(DIR+RES) and without (DIR) a resolved component of the virtual photon. Also shown is the 
partonic prediction from the NLO program DISENT jf8| . Bottom: the single inclusive jet Et cross 
section, corrected to the hadron level, for jets with CMS 2.5 < rj* < 0.5 (photon direction at +z) 
p25j . The events are selected with 9 < Q 2 < 20 GeV 2 and 0.3 < y < 0.6. The inclusive DIS cross 
section in this kinematic range is w 9.8 nb . The data are compared to the HERWIG model with 
only direct interactions (Direct), and with resolved contributions added (DG, where the resolved 



photon structure is given by [ 254 ). The shaded bands show an additional systematic error, mainly 
due to the uncertainty of the calorimeter calibration. 
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The data are compared to different QCD models, see fig. [7.8| . In RAPGAP DIR 
(direct) fl83f , the virtual photon interacts directly with a proton constituent. The 
hardest interaction (with the highest px) takes place at the photon vertex and is 
modelled with the LO QCD matrix element. Softer emissions are generated with 
DGLAP parton showers. In RAPGAP DIR+RES it is assumed that there exists 
an additional contribution from so-called resolved processes, where the structure of 
the virtual photon is resolved by a hard scale, in this case the jet px with p\ > Q 2 
|255| . The hard interaction takes place further down the ladder between (evolved) 
partons from the photon and from the proton. Virtual photon parton densities were 
taken from SaS-2D [[25611 . The resolved contribution breaks the kx ordering in the 
parton evolution. 



a) direct 



b) resolved 



c) CDM 




Jet 
Jet 




Figure 7.8: Parton evolution scenarios, a) "Direct": the hardest emission (with the largest px ), 
given by the QCD matrix element, occurs at the top of the ladder. Further emissions down the 
ladder are ordered with decreasing kx . b) "Resolved" : the hardest emission may occur anywhere 
in the ladder. Softer emissions decrease in kx towards both ends, c) CDM: gluons are emitted 
from succesively built colour dipoles. There is no restriction on kx ordering. 

In the colour dipole model (CDM, program ARIADNE [|86]]), gluon emission 
stems from a chain of independenly radiating colour dipoles. The amount of ra- 
diation depends on the parameters that describe the extendedness of the proton 
remnant colour charge, and on the extendedness of the colour charge of the scat- 
tered quark, which is determined by Q 2 . Emitted partons are not restricted to obey 
kx ordering. The radiation pattern thus shows similarity with BFKL evolution 

The data cannot be described by direct interactions only (fig. |7]7] top), when con- 
ventional DGLAP models are applied. At low Q 2 (and also small x), the measured 
dijet rate is more than a factor two above a model which contains direct interactions 
only. At low Q 2 around 50% of the dijet rate is "missing", that amounts to ~ 2% 
of the total event rate. Also a NLO QCD calculation for parton jets with DISENT 
|78| falls short of the data. Hadronization is unlikely to account for the difference, 
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because the difference between hadronic and partonic dijet rates in the Monte Carlo 
was found to be < 20% for px > 5 GeV. The discrepancy diminishes with increasing 
Q 2 (and also inreasing x). The models where the hard interaction is not tied to the 
photon vertex, implementing either a resolved component of the virtual photon, or 
colour dipole radiation, provide a good description of the data. 

In a similar analysis the single inclusive jet cross section {h? algorithm with 
&cut = 3 GeV ) has been measured as a function of jet E T and pseudorapidity for 
events with < Q 2 < 49 GeV 2 and 0.3 < y < 0.6 |22§. The dat a are corrected to 



the hadron level. It was found that when the transverse energy squared of the jet 
exceeds Q 2 , conventional models (LEPTO 6.5 gl|, HERWIG 5.9 Direct ||) with 



direct photon-parton interactions alone undershoot the data (fig. [7/7| bottom) by a 
large amount. A better description of the data is obtained when contributions are 
taken into account where the hard scale provided by the jet Et serves to resolve par- 
tons inside the virtual photon (HERWIG 5.9 DG, with resolved photon component), 
or by the colour dipole model (ARIADNE 4.08 |§). 

In summary, the mechanism for jet production at small x and Q 2 is not yet very 
well understood. It appears that for regions of phase space where the jet p\ is larger 
than Q 2 the conventional DIS framework, where the virtual photon resolves partons 
inside the proton, ceases to be valid. Instead, the high pt of the jet could resolve 
both structures of the proton and of the virtual photon. We note that the linked 



dipole chain model |24], g9[ should be able to treat both situations consistently where 
the jet p\ is smaller or larger than Q 2 . This field of research will continue with more 
detailed measurements. It has to be seen also in connection with "small-x" physics, 
to be discussed in the next chapter ||. Are these models alternatives to genuine 
"small- re" effects, as expected from BFKL evolution, or are they just a different 
language for the same effect? It remains to be seen whether the competing models 
are actually consistent with the large body of other hadronic final state data, or 
whether they can describe just one facette of the hadronic final state. 
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Chapter 8 
Low-x Physics 



8.1 Introduction 

Amongst the most interesting issues of HERA physics is QCD in the newly accessible 



regime of small x, see section |2.5| . The observed rise of the structure function 
F2 towards small x suggests a strong increase of the parton density in the proton, 
but what is its dynamical origin? Is DGLAP evolution sufficient to account for the 
data, or is BFKL dynamics at work, at least to some extent? The structure function 
measurements are probably too inclusive to answer such questions unambiguously 
Complementary measurements of the hadronic final state provide more detailed 
information on the reaction, which may help to uncover the underlying dynamics. 

In the simple quark parton model, a quark is scattered out of the proton by 
the virtual boson emitted from the scattering lepton. QCD modifies this picture. 
Partons may be radiated before and after the boson-quark vertex, and the boson 
may also fuse with a gluon inside the proton by producing a quark-antiquark pair. 
In fact, the parton which is probed by the boson may be the end point in a whole 
cascade of parton branchings. This parton shower materializes in the hadronic final 
state, allowing experimental access to the dynamics governing the cascade. 

Though there may be other dynamical features to be discovered that leave their 



footprint in the hadronic final state (see for example section (776] on the virtual photon 
structure or chapter |D| on instantons), most dedicated small x measurements have 
concentrated on the predicted signals for BFKL evolution. The leading log DGLAP 
resummation corresponds to a strong ordering of the transverse momenta kx (w.r.t. 
the proton beam) in the parton cascade (Qjj -C k\ x -C ■ ■■k\ i <C ...Q 2 ). In the BFKL 



regime, the transverse momenta follow a kind of random walk (k^ ~ k^ i+1 ) ||257 



The longitudinal momentum fractions are ordered according to Xj ^> for BFKL 
and Xi > Xi + \ for DGLAP. 

A radiated parton (or any other object) with CMS momentum fraction Xi and 
longitudinal and transverse momentum components p Z i,pn is to be found at CMS 
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Figure 8.1: Generic diagram for parton evolution at small x. A gluon ladder evolves between 
the quark box attached to the virtual photon and the proton. The gluon longitudinal momentum 
fractions and transverse momenta are labelled Xi and fcyj. 
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For fixed W and fixed pn, the CMS rapidity yi is determined by lnxj. Here the 
relations Ei = XiW/2 and p Z i ~ — XjW/2 in the CMS have been used (with the 
proton direction along the —z axis). Finally, 
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For DGLAP evolution the phase space for parton radiation between the remnant 
and the current is restricted by kr ordering. A generic signal for deviations from 
DGLAP evolution is therefore enhanced activity in the central and forward CMS ra- 
pidity region, between the current system and the proton remnant. The observables 
which have so far been exploited are the following: 



S j flow: Increased parton activity should result in increased 
transverse energy at central (CMS) rapidity |258| . 



Transverse energy 



Forward jets: High energy jets with Prjet 2 ~ Q 2 (kinematically bound to be 
measured in the forward calorimetric systems) could tag events with BFKL 
evolution, because DGLAP evolution is not allowed 1 12591, 12601 . 
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Hadrons at large px'. High kx partons, disfavoured by the strong kx ordering 



in DGLAP, are signalled by measureable high px hadrons [261 



Other observables: Other observables in the DIS hadronic final state like 
correlations and de-correlations | 262 , [263|| , multijet production |[264| , [265|| and 
photon production [[266|| are being discussed. 



Apart from a few theoretical calculations in the resummed DGLAP and BFKL 
schemes, and calculations in fixed order perturbation theory (mostly NLO), pre- 
dictions for the final states observables can be derived from Monte Carlo models 
only. Theoretical calculations have the advantage that their theoretical input is well 
defined. However, most of them (an exception is presented in section |8l^) neglect 
hadronization, making it difficult to compare them with data. On the other hand, 
the draw back with Monte Carlo models is that often their complexity and their 
flexibility to model hadronization makes it difficult to pin down which feature of 
their theoretical input is actually being tested when comparing to data. 

HERWIG @ and LEPTO which are based upon leading log DGLAP parton 
showers are used as representatives for DGLAP evolution with strong kx ordering. 
Unfortunately a Monte Carlo program based upon BFKL (or CCFM, which in the 
limits of small and large x approaches the BFKL and DGLAP evolutions) dynamics 
is not yet available for ep reactions, though there exist some promising develop- 

|268| , 22]. Rather, the DGLAP model predictions are 

|. In certain aspects 



ments in that direction ]E6 



84 



contrasted with the colour dipole model ARIADNE 
the CDM description of gluon emission is similar to that of the BFKL evolution 
~|. In particular the gluons emitted by the dipoles do not obey kx ordering 



87 



along rapidity (see section |4.4|). 



8.2 Energy Flows 



The measurements of the transverse energy flow ||127| |113| , |114| , |115| , |130| , [269| , |128| , 
270 , 271 , |272|| have had a great influence on our understanding and modelling of 
the hadron production mechanism in DIS. It had been noted already in the early 
papers ||113| , |128|| that some standard QCD models for the hadronic final state failed 



to describe the measured energy flows, in particular at low x (see also section |5.1| ). 
The MEPS model (LEPTO 6.1), based upon the LO QCD matrix element with 
DGLAP parton showers and string fragmentation, undershoots the HI data |113|] 
by a factor two in the "forward region" , towards the proton remnant (fig. |3.2| a) . 
In the following years, the energy flow and it's dependence on kinematic conditions 
has been studied in great detail. 

The "forward" measurement close to the beam hole at the edge of the calorimeter 
acceptance at r/i ab « 3.5, corresponding to 9\ ah w 4°, is very difficult, as particles 
are lost in the beam pipe, and backscatter from lower angle particles hitting dead 
material produce background in the forward calorimeter region ][273|, [274j, |275||. An 
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independent measurement is provided by ZEUS, where the instrumentation and 
dead material distribution in the forward region are quite different. Fortunately, the 
data from ZEUS [[27C|| agree with the HI data (see figs. |8.2| b and |3.4|) , confirming 
the "forward energy crisis" |[276| . 
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Figure 8.2: The transverse energy flow Et as a function of pseudorapidity 77 in the laboratory 
frame. The proton direction is to the right, a) HI data compared to the models MEPS (LEPTO 



6.1) and CDM (ARIADNE 4.03) JllJ. b) The measured forward E T flow from HI |113[ and 
ZEUS |27(| compared to the models MEPS (LEPTO 6.1) and CDM (ARIADNE 4.03, here labelled 
CDMBGF), and to BFKL and DGLAP based theoretical calculations (parton level) ]258fl . The 
data are corrected for detector effects. 

The discrepancy between data and LEPTO, considered as the reference model 
for hadron production, triggered further investigations, both theoretically [|258| , |277| 
and experimentally [|114| , |115| , |270| , |271|| whether low x effects, for example BFKL dy- 
namics, could be held responsible for the excess. Shown in fig. [B.2| b are also results 
of calculations [2"58"| based upon either BFKL or DGLAP parton dynamics. Qualita- 
tively, much more Et is expected from BFKL than from DGLAP evolution, roughly 
at the level of the data. However, these calculations do not include hadronization, 
and are therefore not directly comparable to the data. 

The Et flows have now been measured over a wide range of x and Q 2 fll 1| . 



115|, |271|| , and are presented in the hadronic CMS to eliminate the transverse boost 

The latest preliminary HI data 



115, 278 



introduced by the scattered electron 
(fig. |8.3|) also employ the plug calorimeter [579] installed at small angle behind the 



liquid argon calorimeter, extending the acceptance far into the target region. The 
plug data have large errors due to the large amount of dead material along the line 
of sight to the interaction vertex, between 2 and 4 nuclear absorption lengths. 

The Et flow is plateau-like (~ 2 GeV per unit if) at small x and Q 2 , and be- 
comes more peaked in the current region at larger x and Q 2 . The average Et ((Et)) 
measured centrally (—0.5 < i] < 0.5) in the CMS rises with decreasing x (see 
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fig. This behaviour is predicted for the partonic (Et) from the BFKL cal- 

culation, whereas for DGLAP evolution the opposite is expected ||258|| . CDM gives 



still the best overall description of the data. The DGLAP based models, MEPS and 
HERWIG, have developed and now give also a reasonable description of the data 

In order to draw conclusions from comparing the data with partonic QCD calcu- 
lations, hadronization effects need to be estimated from Monte Carlo models. The 
partons produced in the CDM give about twice as much Et in the central rapid- 
ity region as the ones emitted from the ordered cascades in MEPS and HERWIG, 
see fig. BJ3 a [ 261 |, which can be traced to the fact that hard gluon radiation is 
much more abundant in CDM, see fig. |S.5|b [ 281 1 . However, the observable parti- 



cles emerging after hadronization give rise to very similar Et flows. The Et flow 
predictions of the different models cannot be well discriminated with the present 
data. While hadronziation adds relatively little (~ 0.5 GeV /unit rapidity) to the 
partonic Et for CDM, most of the Et (~ 1.4 GeV /unit rapidity) is generated by 
hadronization for MEPS and HERWIG. It is noteworthy that CDM and MEPS use 
the same Lund string fragmentation model for hadronization. In addition, the par- 
tonic and the hadronic Et are well correlated event-by-event in CDM, as opposed 
to the other models 
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With the Et data alone, the hadronization ambiguity cannot be resolved. The 
situation is summarized in fig. S.6 . The (Et ) data rise with falling x and agree with 
both the ordered and the unordered Monte Carlo models for hadrons. BFKL predicts 
also a rise, as is seen in the CDM partons, but at a somewhat lower level. Apart from 
the normalization, the data are consistent with a BFKL interpretation, assuming 
hadronization effects as given by CDM. The DGLAP model partons produce much 
less Et than CDM, and have the opposite x behaviour, as expected for DGLAP 
evolution [258]. However, the data are also consistent with DGLAP evolution, when 
large hadronization effects are assumed. How do they arise? 

In LEPTO 6.4 changes to the non-perturbative part were introduced to simul- 
taneously describe the large Et flow seen in the data, and to explain the rapidity 
gap events ||104| , |105| , |102j , |103|| without invoking the explicit exchange of a Pomeron 
1 09 1 . In the new concept of soft colour interactions (SCI), colour quantum numbers 



may be exchanged between the outgoing partons after the hard interaction, leading 
to a reconfiguration of the fragmenting strings (see section [4.6|) . This may result in 
colour neutral subsystems, separated by a rapidity gap from the rest of the event 
after hadronization, or it may enhance the Et flow, when strings are spanned back 
and forth (fig. 4.14 ). Also the treatment of the remnant fragmentation has been 
changed for events where a sea quark is hit (new sea quark treatment - SQT). In- 



1 The x dependence for fixed Q 2 implies also a W dependence. In fig. 5.8 it is shown that the 
average Et at central rapidity increases with W in DIS, and that this behaviour is in agreement 
with data from hadron-hadron collisions. A deeper understanding of Et production in DIS might 
also be benefitial for the understanding of the hadron-hadron Et data, which so far could be just 
parametrized phenomenologically. 



120 




H1 94 Preliminary 
ARIADNE 
HERWIG 
LEPTO 



-5 5-5 5 

CMS ti* 



Figure 8.3: The transverse energy flow Et as a function of pseudorapidity 77 in the hadronic CMS. 



The proton direction is to the left. The data [115 are compared to the models CDM (ARIADNE 
4.08), MEPS (LEPTO 6.4), and HERWIG 5.8. Shown are statistical errors only, except for the 
two foremost data points from the plug. 



stead of recombining the sea quark partner (for example a s quark in case a sea 
quark s was hit) within the remnant to form a hadron, a fragmenting string is now 
stretched between the hit sea quark and the remnant. At large x these effects are 



minor, but at small x they are significant, see fig. |8.7| . As long as the new SCI 
concept has not been ruled out by other data - from detailed studies of rapidity 
gap events, for example - it has to be taken seriously. The SCI does not only offer 
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Figure 8.4: The mean E T in the CMS pseudorapidity interval —0.5 < r\ < 0.5 as a function 
of x for different Q 2 ranges [28C]. Besides the preliminary data from HI and ZEUS | |27f| 

for hadrons, the BFKL expectation for partons is shown [25S]. Hadronization effects have to be 
taken into account when comparing this data with the partonic calculation. For the HI data an 
additional 8% overall normalization uncertainty is not shown. 



an explanation for the rapidity gap events and the large Ex flow, but has also been 
advocated 



to explain the large production rate of heavy quarkonia in pp inter- 
actions at the Tevatron [[283)1 ■ Presently the SCI model has difficulties describing the 
thrust of the hadronic final state in rapidity gap events, but so do the other models 
based upon Pomeron exchange ||112|| . Further data will be necessary to discriminate 
the Pomeron exchange and SCI models. 

In the CDM hadronization effects are much smaller, because gluon radiation is 
more abundant, and less energy remains for hadronization. The large hadronization 
effects in HERWIG can possibly be explained by its very simple cluster fragmenta- 
tion scheme. When there are few partons, the masses of parton pairs (clusters) to 
be fragmented tend to be large. In the cluster fragmentation model, they undergo a 
two-body decay according to phase space, leading to large transverse energies. The 
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Figure 8.5: a) Transverse energy flows as a function of the pseudorapidity rj for partons and for 
hadrons. b) The multiplicity of hard gluons with pr > 2 GeV vs. 77. The events are generated 
with the models CDM (ARIADNE 4.08), MEPS (LEPTO 6.4) and HERWIG 5.8 in a "low x " 
kinematic bin with (x )—0. 00037 and (Q 2 ) = 13.1 GeV 2 . The proton direction is to the left. 
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Figure 8.6: The mean Et in the CMS pseudorapidity interval —0.5 < 77 < 0.5 as a function of 
x for Q 2 — 14 GeV 2 . Besides the data from HI [115] (which have an additional 8% normalization 



uncertainty) and ZEUS J271| , the BFKL expectation for partons is shown [258|, together with the 
results from the models CDM (ARIADNE 4.08), MEPS (LEPTO 6.4), and HERWIG 5.8 for 
hadrons and for partons. The LEPTO predictions for hadrons without the features soft colour 
interaction (SCI) and the new sea quark treatment (SQT) are also shown. 



same mechanism leads to a significant fraction of rapidity gap events which are "not 
exponentially suppressed" in HERWIG without Pomeron exchange or soft colour 
interactions. In a later version (HERWIG 5.9) this problem is being addressed by 
allowing large mass clusters to split longitudinally before they decay. 
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Figure 8.7: The E T flow from HI gl| at (x) « 0.002 (left) and 0.0004 (right) for (Q 2 ) w 14 
GeV 2 . The data (shown are statistical errors only) are compared to the MEPS model (LEPTO 
6.4, structure function MRSH) in the standard version, without soft colour interactions (SCI), and 
with neither soft colour interactions nor the new sea quark treatment (SQT). 

8.3 Transverse Momentum Spectra 

Not yet well understood hadronization effects, however interesting they may be, pre- 
cluded strong conclusions on the underlying parton dynamics from the transverse 
energy flow measurements. It has been shown that single particle transverse mo- 
mentum (jprp ) spectra represent a more direct measure of the partonic activity, and 
corresponding measurements have been suggested [284, |261[ |. 

To answer the question whether the Et observed in the data is generated predom- 
inantly by parton radiation or by hadronization, inclusive pt spectra are considered. 
Hadronization should produce typical spectra which are limited in pt , while parton 
radiation should manifest itself in a hard tail of the pt distribution. To test this idea 
particles from a "central" (in the CMS) r\ interval < rj < 2 are examined. Simu- 
lated events are compared which have similar hadronic Et (E T a,d ) in that interval, 
with either a large or small contribution E T Ar from partons. The parton dominated 
events indeed exhibit a harder pt spectrum than the hadronization dominated ones 
(s. Fig. |8.8| a), regardless of the mechanisms employed for parton evolution and 
hadronization. Different parton evolution scenarios which differ in the production 
of hard gluons (for example BFKL vs. DGLAP, or CDM vs. MEPS/HERWIG, 
see fig. |8.2| ) can thus be discriminated by the number of high-p^ hadrons which are 
produced. In a similar fashion the distribution of the Et measured centrally could 
be used for discrimination ||261|| (fig. |8.8| c). Preliminary HI data [ |115|| favour the 
CDM over the other models, but the study of systematic errors has not yet been 
completed. 

HI has measured 



118] the charged particle Pt spectra as much central in the 
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Figure 8.8: a) pt spectra of charged particles with < rj < 2 for events which are hadronization 
dominated (left, ££ ar < 0.2 GeV/unit rapidity) or parton dominated (right, E^ 1 /E^ ad > 0.5) 
]26l| ]. The events are required to have E^ ad between 1 and 2 GeV/unit rapidity, c) The 
Et distribution (measured in < rj < 2) at (x) = 0.00037 ]26l| . The events are generated 
with CDM (ARIADNE 4.08), MEPS (LEPTO 6.4) and HERWIG 5.8. 



CMS (0.5 < rj < 1.5) as the detector acceptance allows^. Measured is 

Idn f , dW(e P -> e'Xh) / d 2 a(ep ^ e'X) 

' ' d V—r^ — 3 ^ / > ( 8 - 3 ) 



iV dpr -'Ar? d^dprda;dQ 2 / dxdQ 2 

where the differential cross section to produce a hadron with transverse momentum 
Pt in a pseudorapidity interval Arj is normalized to the differential event cross section 
for given x, Q 2 . In the measurement, the differential cross sections are averaged over 
certain small x, Q 2 bins. 

The data are compared in fig. 1813 to models with suppressed (LEPTO 6.4, 



including SCI and new SQT; HERWIG 5.8) and unsuppressed (ARIADNE 4.08) 
radiation patterns. At "large" x and Q 2 (here "large" means Q 2 ~ 35 GeV 2 and 
x ~ 0.004) all models provide a good description of the data. At smaller x and Q 2 , 
LEPTO and HERWIG fall significantly below the data for p T > 1 GeV. ARIADNE 
gives a good description of the data over the full kinematic range. The shortfall 
of the models with suppressed gluon radiation indicates that at small x there is 
more high kj> parton radiation present than is produced by the models based upon 
leading log DGLAP parton showers. These models are closer to the data in the 
current region, or when all charged particles are considered, including the soft ones 
TTBI , see fig. £7[ 

In a preliminary analysis of the HI forward tracker, p? data have been obtained 



also in the "truly central region", —0.5 < rj < 0.5. In fig. 8.10 the charged particle 



2 The measurement relies mostly on the central tracking device; when systematic effects presently 
hampering the forward tracker are mastered, the measurable rapidity region could be enlarged by 
one unit of rapidity towards the remnant. 
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Figure 8.9: The pr spectra of charged particles, measured in the CMS in the pseudorapidity 



interval 0.5 < r\ < 1.5 118 . Data are shown for nine different kinematic bins with the mean values 
of x and Q 2 as indicated, plus the combined sample (bin 0). The models ARIADNE 4.08, LEPTO 
6.4 and HERWIG 5.8 are overlayed. 



126 



Pt spectra are shown for three different r\ ranges. At central rapidity, —0.5 < rj < 
0.5, the spectra from photoproduction and DIS are similar. Towards the photon 
fragmentation region, the DIS spectra become harder than the photoproduction 
data. Apparently, in that region, the virtuality Q 2 influences the hardness of the 
Pt spectrum, while there is little influence in the central region. A similar conlusion 
has been arrived at from the study of the flow, see section EO. The DIS data 



are well described by the CDM. LEPTO becomes too soft away from the current 
region. The prediction by the photoproduction model PHOJET | |285| | describes the 
photoproduction spectra well. 



HI Preliminary 



0.5<r| <1.5 




3 4 5 

p T */GeV 



190 GeV 



HI <Q 2 >=16GeV 2 
HI <Q 2 >=0 GeV 2 
_ CDM <Q 2 >=16 GeV 2 
. . . LEPTO <Q 2 >= 1 6 GeV 2 
PHOJET <Q 2 >=0 GeV 2 



Figure 8.10: The pt spectra of charged particles for different CMS r\ ranges (Hf preliminary, 
||286f ). DIS data with (Q 2 ) = 16 GeV 2 are compared with photoproduction data with (Q 2 ) = 0. 
Also shown are the DIS models CDM (ARIADNE 4.08) and LEPTO 6.5, and the photoproduction 
model PHOJET 1.3. The data samples are selected such that they have the same (W) = 190 GeV . 

Theoretical calculations for the charged particle pt spectra have become avail- 
able |[287|| , where a cross section <jj to produce a parton j is folded with an exper- 
imentally known fragmentation function Dh/j(z, /J 2 ) |[2 1 1|] to produce a hadron h 
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with momentum fraction z from the parton j, see fig. RTTTIa. Monte Carlo models 
for hadronization are thus avoided. Symbolically, 



a h = Oj ® D h/j . (8.4) 

The parton cross section <jj is divided into a conventional part where gluon 
emissions in the ladder are ordered, and a part with unordered emissions, which is 



calculated with the BFKL equation (see fig |8.11| a). The HI px data at low x are well 



described by this calculation, see fig |8.11| b. The absolute normalization for the BFKL 



part is obtained by requiring that the calculated parton cross section would match 



the measured HI forward jets j288[ at the hadron level, see section |8.4| . We note 
that this normalization neglects the effects of the jet algorithm and hadronization 
effects, both of which are not included in the single parton cross section calculation. 
Without the BFKL part the calculation is by a factor « 2 below the data. Both 
predictions show little dependence on the choice of the factorization scale /i 2 for the 
fragmentation function. 

For the "no BFKL" calculation the LO cross section for three partons was used, 
where the main contribution is the diagram with the qq from the quark box plus 



one additional gluon, see fig. |8.11| a. A complete NLO calculation is still missing. 
DGLAP evolution was neglected, because the kinematic conditions severely dis- 
favour DGLAP evolution. DGLAP evolution is not possible between the parton j 
and the photon if kxj ~ Q due to strong kx ordering. For the curves shown in fig. 
8. 1 l| b this condition is approximately fulfilled, Q 2 /2 < k T j < 2Q 2 , because a lead- 



ing hadron carries on average about half of the fragmenting parton's momentum, 
Pt ~ k Tj /2. 

An even better signal is expected when the single particle detection capability 
could be extended further towards the proton remnant [ |281| , |266| | . Typically at 
HERA, the central trackers are limited to 8^ > 20°, and the forward trackers to 
#iab > 7°. An upgrade with a forward silicon tracker could cover angles 9y ah > 3° 
28 1|| . Alternatively, one can detect single photons or tc°s in the calorimeters, which 



have acceptance for typically #i a b > 4°. 

HI has measured "forward" n° production |[289|| , where the 7r° meson is detected 
via the electromagnetic showers of their decay photons in the HI forward liquid 
argon calorimeter [25S]. The 7r°'s have been selected in the laboratory frame with 



5° < #iab < 25°, E > 8 GeV , p T > 1 GeV and x n := E/E p > 0.015 (data with 
other x n cuts exist as well) in order to enhance the phase space for BFKL evolution. 
The number of such 7r°'s increases with decreasing Bjorken x, see fig. |8.12j . The 



models studied, MEPS and CDM, agree at large x with the data. The rise towards 
small x is only followed by CDM. Calculations for forward 7r° production exist in 
principle ||266| , |287|| , but need to be repeated for the kinematic conditions of this HI 



measurement. 

Coherence effects, which can be taken into account by angular ordering, are 
expected to modify the BFKL behaviour. They are included in the CCFM equation 
|20|| . Predictions for the hadronic final state at low x based upon the CCFM ansatz 
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Figure 8.11: a) The production of a hadron h from a parton j through a fragmentation function 
D. A gluon ladder connects the quark box attached to the virtual photon with the proton. The 
encircled part of the graph is calculated either with or without BFKL evolution, b) The result of 



the calculation with and without BFKL calculation 287 , compared to the HI data at low x [118 
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Figure 8.12: The number of "forward" ir° mesons as a function of Bjorken x. N is the number of 
events in the selected kinematic region. The preliminary HI data [ 289 1 are compared to the MEPS 
model (LEPTO 6.5 §3, including SCI and the new SQT) and the CDM (ARIADNE 4.08). 



are now emerging 



22 



In general, there is less hr diffusion in CCFM than in 
BFKL. The number of soft emissions is much reduced in CCFM when compared to 
BFKL, due to angular ordering (fig. |8.13| ). For appreciable transverse momenta, the 
differences are not so big, affecting mainly the high multiplicity tail. 
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Figure 8.13: The probability distribution P(n) for the number n of gluon emissions with 
transverse momentum a) kx > 1 GeV and b) fc^ > 0.007 GeV. The calculations in the DGLAP, 
CCFM and BFKL schemes §| are for events with x = 5 • 10~ 5 . 



The CCFM equation has been reformulated in the linked dipole chain model 
(LDC) P^fl . First results from a Monte Carlo generator implementing the LDC in 
the ARIADNE framework are promising [89]. The critical distributions, Et flows 



and pt spectra at small x, are much better described than with LEPTO, though 
the agreement with the data is not quite as good as with the CDM |8£ 



8.4 Forward Jets 



The classic signature for BFKL evolution in the haSronic final state are so-called 
"forward jets" [|259| , |260|| close to the remnant, that is "forward" for the HERA 
detectors. The selection criteria are designed such as to suppress the phase space 
for DGLAP evolution, and to maximize it for BFKL evolution. Xj et = E- ]et /E p , the 
ratio between the jet energy £j e tand the proton beam energy E p , is required to be as 
large as possible, since BFKL evolution between the jet and the quark box re ires 
Xjet ^ x tQ allow for strong ordering x i+ i 3> Xi along the ladder (see fig. |8.1| , where 
the forward jet is thought to emerge from a radiated gluon close to remnant). The 
transverse momentum k^et has to be close to Q (£/rjet ~ Q)- The phase space for 
kr ordered DGLAP evolution (£>ri+i 3> kxi) is thus much reduced, whereas it is still 
open for unordered BFKL evolution. An enhanced rate of events with such jets is 
thus expected in the BFKL scheme ||259|| . The experimental difficulty is to detect 



jets close to the beam hole in the proton direction in the forward calorimeter. 

Previous measurements by HI [|1 14|| were limited in statistics, but were consistent 
with BFKL calculations ||290[ |. In the preliminary analysis of the larger statistics 



sample from 1994 ||288|| forward jets with Xj et > 0.035 are selected with the cone 
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algorithm (cone radius R = V At]' 2 + A0 2 = 1) in the HERA laboratory frame within 
7° < 9j e t < 20°. The transverse jet momentum requirements are 0.5 < PTjct 2 /Q 2 < 2 
and a lower cut-off prjet > 3.5 GeV. The forward jet rate, corrected for detector 
effects to the hadron level, increases sharply with falling x, see fig. ETT3 . This is 
expected from BFKL calculations, in contrast to calculations without the BFKL 
ladder |291| . The behaviour of the data is well described by the CDM with the 
unsuppressed gluon radiation pattern. The DGLAP based model MEPS cannot 
describe the rise towards small x when soft colour interactions are switched off. 
When soft colour interactions are activated however, without changing the parton 
dynamics, MEPS comes up rather close to the data. A good description of the 
forward jet rate is also provided by RAPGAP including a resolved (not point-like) 
component of the virtual photon [|255|1 . The physical significance of this fact is still 
being debated. 
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Figure 8.14: a) The forward jet cross section as a function of x from HI f288fl (1994 data). The 
data arc corrected to the hadron level and compared to the models CDM (ARIADNE 4.08) and 
MEPS (LEPTO 6.4 with and without soft colour interactions), b) The forward jet cross section 
from a NLO calculation and for parton jets found in CDM and MEPS, compared to a BFKL 



parton cross section calculation [291]. c) The ZEUS forward jet cross section |292j (1994 data) vs. 
x , corrected to the parton level |278fl . x- ]ct > 0.035, 0.5 < pTjet 2 /Q 2 < 4 and prjet > 5 GeV were 
required. The data are compared to a NLO jet calculation ]293[ and to parton jets from CDM 
(ARIADNE 4.08). Also shown are parton cross sections (no jet algorithm) with ("BFKL") and 



without ("Born") BFKL evolution [291]. The systematic errors do not include uncertainties due 
to hadronization. 



Hadronization effects may play an important role for small prjet, where it is 
possible to form a hadron jet from hadronization fluctuations which is uncorrelated 
with any parton. Only for p T j Ct > 7 GeV hadronization effects were found to be at 
the 10% level |[294|| . However, as one is leaving the BFKL condition p T j et Q when 
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increasing the lower cut-off, the BFKL signal diminishes. This could be compen- 
sated by extending the measurements further into the forward region with upgraded 
detectors ||295|| . The hadronization effect for jets depends on the Et produced in 



hadronization. In the CDM, where most of the is produced by parton radiation, 
hadronization effects are much smaller. The MEPS forward jet rate for partons 
agrees with a NLO calculation with the program DISENT [7S (see fig. |8.14|). The 



BFKL calculation for partons shows an even larger increase towards small x than 
the CDM partonic forward jets. Similar conclusions were reached in a preliminary 
ZEUS analysis of 1994 data p92] , see fig. |8TT4l c. 

First results from a larger data sample from 1995 have been obtained, allowing 
for a larger PTjet cut. In the preliminary ZEUS analysis ||296|| events with 0.00045 < 
x < 0.045, y > 0.1 and E' e > 10 GeV are selected with (Q 2 ) « 15 GeV 2 . The jets 
are required to lie in the target region of the Breit frame and to satisfy Prjet > 5 GeV, 
Xj e t > 0.036 and 0.5 < PT}et 2 /Q 2 < 2. In the laboratory frame 0j et > 8.5° is required. 
Here we discuss the analysis with the cone algorithm (radius R = 1). An analysis 
with the kx algorithm gives larger jet rates, but one arrives at similar conclusions. In 
a first step the data are corrected for detector effects. At the hadron level (fig. |S.15| a), 
the increase towards small x of the jet rate is well described by ARIADNE 4.08. 
LEPTO 6.5 falls short of the data by a large amount at small x. 

In a second step, the data are corrected for hadronization effects to the parton 
level with ARIADNE, see fig. |8.15[ b. Uncertainties due to this correction have not yet 
been studied. The data follow the expectation for jets formed from the ARIADNE 
partons. The partonic jet rate at low x is far above a NLO calculation. The data 
are also compared to a calculation with and without BFKL evolution, though these 
calculations do not include a jet algorithm. However, the NLO result is close to 
the Born graph calculation, which can be interpreted as a hint that the differences 
are not so large. The calculation without BFKL evolution yields a low rate close to 
the NLO calculation. With BFKL effects included, a much larger rate at small x 
is expected, that even overshoots the data. However, in view of the possibly large 
hadronization uncertainty and the missing jet algorithm in the calculation, no firm 
conclusion could be drawn from this comparison |[296| 1. 



8.5 Summary and Outlook 

The hadronic final state has been searched for an unsuppressed parton radiation 
pattern, which could for example be expected from BFKL contributions, in contrast 
to pure DGLAP evolution. The high level of transverse energy is consistently de- 
scribed when BFKL evolution is included, but possibly large hadronization effects 
could also mimick the effect. In order to describe the data with DGLAP evolution, 
hadronization phenomenology has introduced new concepts like soft colour interac- 
tions in LEPTO, or it involuntarily produces "not exponentially suppressed rapidity 
gaps" by cluster fragmentation in HERWIG. 
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Figure 8.15: The forward jet cross section as a function of x from ZEUS |296| . Statistical and 
systematic errors are added in quadrature. The shaded bands give the additional systematic error 
due to the uncertainty of the hadronic energy scale, a) The data are corrected to the hadron 
level, and compared to the expectations from ARIADNE and LEPTO. b) The data are corrected 
to the parton level. Systematic errors due to the hadronziation correction are not included. Also 
shown are the ARIADNE and the NLO predictions for jets [jT9|, and the BFKL and Born graph 



calculations for parton cross sections [291] 



The forward jet measurements at low p? are in rough agreement with BFKL ex- 
pectations, but are also plagued by hadronization uncertainties. With larger statis- 
tics data sample, the measurements could be made at higher pt where a closer 
relation between hadron jets and partons is expected |j294|| - From the measured 
Et flows and forward jets, a suppressed radiation scenario is presently only tenable 
when large hadronization effects are assumed. 

The measured single particle px spectra require more parton radiation than is 
expected from pure DGLAP evolution, regardless of the hadronziation model. The 
spectra have been calculated via known fragmentation functions in the DGLAP and 
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the BFKL scheme. In the DGLAP scenario the calculation by far underestimates 
the data. 

Whereas pure DGLAP evolution fails, BFKL evolution offers a plausible expla- 
nation of the measurements. This does of course not preclude other explanations. 
For example, it has been suggested that contributions from resolved photons, where 
the hard interaction happens not at the photon vertex but further down the ladder, 
could be responsible for the additional parton activity ||255| , [297|| . We have seen 
possible signs for such effects also in the inclusive jet rates, see section |7.6| . More 
refined measurements to unravel the details of parton dynamics at low x, including 
correlations between jets and hadrons, have been proposed |58|, 263 1. It will also be 
important to check the consistency of the competing models by confronting them 
with the entire data on the hadronic final state. 

BFKL effects are also being looked for in pp collisions at the Tevatron. A jet-jet 
decorrelation is seen with increasing rapidity difference, though the effect is larger 
than expected from BFKL, and is described by the HERWIG and PYTHIA models 
without BFKL dynamics ||298| . Recently, BFKL effects have been brought forward 
to explain the discrepancy between measured jet rates and NLO calculations at 
small jet transverse energies for different CM energies at the Tevatron |299 
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Chapter 9 
Instantons 



9.1 Introduction 

For a long time it has been recognized that the standard model contains processes 
which cannot be described by perturbation theory, and which violate classical con- 
servation laws like baryon number (B) and lepton number (L) in the case of the 



electroweak interaction, and chirality (Q$) in the case of the strong interaction ||300 



Such anomalous processes are induced by so-called instantons [ |301|| . The name indi- 
cates that these are non-perturbative fluctuations that are confined to "an instant" 
in space-time, with no corresponding free particle solutions for t — > ±00. The 
interest in instantons remained somewhat academic, as observable effects were pre- 
dicted to exist only at extremely high energies of O(10 5 TeV), until it was discovered 
that their exponential suppression is much reduced by the emission of gauge bosons 
|302| . In electroweak theory with massive gauge bosons still rather high energies of 
0(> 10 TeV) would be required, but not so in QCD with massless gluons and strong 
coupling. Instanton effects could play a role in QCD reactions already at present 
day colliders. Deep inelastic ep scattering at HERA is particularly interesting, be- 
cause the virtuality of the photon probe Q 2 provides a hard scale for the instanton 
subprocess, which is needed for theoretically sound predictions ||303| , [304], |305|| . In- 
stanton effects have not yet been observed in nature. Their experimental discovery 
would be of fundamental significance for particle physics. 

Here a short introduction to the basic theoretical ideas will be given (a pedagogi- 
cal treatment of instantons can be found in [ |30(j| |). Instanton phenomenology in deep 
inelastic scattering (DIS) will be discussed, covering cross sections and event topolo- 
gies. Finally, prospects for instanton searches and first results from the analysis of 
HERA data will be presented. 
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9.2 Instanton Theory 



Instantons originate from the non-trivial topological structure of the vacuum in non- 
Abelian gauge field theories, where the vacuum is degenerate in the Chern-Simons 
number iVcs- Ncs is defined as an integral over the gauge fields A a with coupling g, 



N, 



cs 



16vr 2 



d 6 x£ijk ( Aj djA k — —eabcAj A^AI 



(9.i; 




Neighbouring vacua have the same (minimal) 
potential energy, but differ in their topological 
winding number Ncs , and are separated by a 
potential barrier of height Eb (fig- |9.1|) . The 
usual perturbative expansion of the scattering 
amplitudes in the coupling constant a around 
one minimum (fig. |9.1| ) , conveniently represented 
by a series of Feynman graphs, does not allow for 
transitions between neighbouring minima. They 
may however occur classically when the energy E 
is large enough E > Eb, or by quantum mechan- 
ical tunnelling when E < Eb, corresponding 
to so-called instanton solutions of the classical 
field equations. The transition amplitude for the 

instanton-induced tunnelling process is exponentially suppressed oc exp(— An /a) 
very small number. 



Fig. 9.1. The structure of 
the vacuum. Instanton solu- 
tions represent tunnelling tran- 
sitions between topological in- 
equivalent minima, which cannot 
be reached perturbatively. 




Figure 9.2: a) the electroweak interaction with A(B + L) = —6 and in b) the strong interaction 
with AQ 5 = 8. 

In the electroweak theory, the minimal barrier height is Eb ~ m^ja w = O(lOTeV). 
Instanton transitions between vacua separated by ANqs (see fig. |9.2j a for an exam- 
ple) would violate baryon (B) and lepton numbers (L = L e + + L T ) according 
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to 

A{B + L) = -2 n gener . • AN C s , (9.2) 

but respect 

A(B-L) = AL e = AL^ = AL T = AB/3. (9.3) 

Ngener. = 3 is the number of fermion generations. 

In instanton induced QCD reactions (see fig. |9.2[ b) chirality is violated. The 
chirality Q5 is the difference between the number of left- and right-handed fermions, 
Q 5 = j^L — #R. For nj active quark flavours, the selection rule is 

AQ 5 = 2 n r AN CS . (9.4) 

The minimal barrier height is given by the hard scale of the process, e.g. Eb = 0(Q) 
for DIS [304]. The exponential suppression is less severe than in the electroweak case, 
because a s > « w . 



9.3 Instantons at HERA 



In recent years, it has been realized [304, 305 , |307 , |308| , 309 , 310] that quantitative 



calculations are possible for processes induced by QCD instantons in DIS due to 
the presence of a hard scale, Q 2 . In DIS, events due to QCD instantons I (and 
anti-instantons /) are predominantly produced in a photon-gluon fusion processes^ 
(Fig. El) 

1* + 9 ^^2(Qr + Qr) + n g9 1* + 9 + + n g g. (9.5) 

nf nf 

In each event, quarks and antiquarks of all rif active flavours are found, with n g 
gluons in addition. 

The kinematics is depicted in fig. |9.3| . The DIS variables Bjorken x and Q 2 can 
be measured from the scattered electron, q = e — e'. A measurement of the other 
variables is more challenging. A measurement of the invariant mass of the hadronic 
system, excluding the remnant, would determine s. If the outgoing "current jet" 
could be identified and measured, it's 4-momentum q" would determine q' = q — q", 
and thus the variables x' and Q' 2 which characterise the instanton subprocess. In 
practice, when not all of the five independent invariants (for example x, Q 2 , x', Q' 2 , s) 
can be measured, they are being integrated over. 

The instanton induced cross section is given by a convolution of the probability 
to find a gluon in the proton P g / P , the probability that the virtual photon splits into 

1 Quark initiated processes have not yet been considered. Due to the large gluon content of the 
proton in the HERA domain at small x, they are expected to be of minor importance. In addition, 
they are expected to be suppressed by 0{a 2 s ) with respect to the gluon initiated processes. 
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P 




DIS variables: 

Q 2 --=-q 2 

x:=Q 2 /(2P-q) 

W 2 := (q + P) 2 = Q 2 (l-x)/x 

§ ■= (q + g) 2 

i = x(l + s/Q 2 ) 

Variables of instanton subprocess: 

Q' 2 :=-q' 2 

x> := Q' 2 /(2 g ■ q') 

s > ■= ( q > + gf = Q* (l- x ')/x' 



Figure 9.3: Kinematics of instanton induced processes in DIS. The labels denote the 4- vectors 
of the particles. A virtual photon 7* (4-momentum q — e — e') emitted from the incoming electron 
fuses with a gluon (4-momentum (?) from the proton (4-momentum P). The gluon carries a 
fraction £ of the proton momentum. The virtual quark q* entering the instanton subprocess has 
4-momentum q' ', and the outgoing quark from the 7* — * qq splitting has 4-momentum q" . The 
invariant masses squared of the 7*3 and q*g systems are s and s' . W is the invariant mass of the 
total hadronic system (the -f*p system). 0<a;<a;/^<x'<l holds. For completeness, we note 
y := (Pq) / (Pe) = Q 2 /(sx), where s = (e + P) 2 is the ep invariant mass squared. 



a quark-antiquark pair in the instanton background P g ,L», and the cross section 



a q * g {x' , Q' 2 ) of the instanton subprocess [p05| , [3071| . Multi-gluon emission enhances 
the cross section ||302 | 



(I) 

a q*g;n g (x 



n 



9- 



exp(— 4:ir/a s 



(9.6) 



The cross section of the instanton induced subprocess is then |p05 



a 



(I) 
1*9 



x',Q 



f2\ 



E 



(-0 

q g;n g 



47T 



Efa;') / 



■2± 

2 



-47T 



cxp 



F(x') 



(9.7) 



It depends critically on the functions F(x') (called the "holy grail" function), which 
modifies the exponent in the suppression factor exp(— 4ir/a s ), and on S(x'), which 
depends on F(x'). There exists also a scale dependence due to the choice of the 
renormalization scale fi{Q')- 
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F(x') can be estimated reason- p( x ') ^ 



ably well (see fig. |9.4|) for x' not too 
small, 



ex- 
x' is 



x' > 0.2 pQ5l . The 
trapolation to lower values of 
unreliable due to inherent ambigu- 
ities. In addition, multi-instanton 
effects should be avoided by limit- 
ing the instanton size pi (the spa- 
tial region occupied during the in- 
teraction) to pi < 2 GeV -1 with a 
25 GeV 2 [BDSL |307f 




> 



cut-off Q' 2 
That requirement ensures also that 
a s (p(Q')) stays small enough to ap- 



Fig. 9.4 
small s' 
applied. 



The holy grail function F(x') [|305|. For 
(x' ps 1), instanton perturbation theory is 
The calculation with the valley method 



matches smoothly with the perturbative result. 



ply instanton perturbation theory. 

The resulting instanton induced 
subprocess cross section a q l g (x', Q' 2 ) 

(see fig. PT5| ) is peaked at Q' « 5 GeV and exponentially grows with decreasing x' . 
The integrated instanton induced ep DIS cross section (see fig. |9l]) is sizeable; for 
x > 0.001 and x' > 0.2 it is of O(10 pb). The cross section is approximately scaling 
(depends only on x, not on Q 2 for large Q 2 ) ||307|| . It grows towards small x, and 



increases dramatically when the lower x' cut-off is relaxed. Eventually higher order 
instanton effects have to dampen the growth of the cross section. 
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a) 

Figure 9.5: The instanton subprocess cross section [ |309| for q*g — > hadrons as a function of a) 
x' for Q' = 5 GeV and b) Q' for x' = 0.2. 



Two kinematic regions have to be distinguished. For x' > 0.2 the predictions 
are relatively safe, allowing the instanton theory to be tested. Either instantons are 
discovered at the predicted level - including the substantial theoretical uncertainties, 
which still need to be quantified -, or the theory has to be revised. For x' < 0.2 
the cross section presumably continues to grow, but the extrapolation is extremely 
uncertain. For a discovery, this is the favourable region due to the large cross 
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Figure 9.6: The instanton induced DIS cross section [ 309 1 ep — > e'X integrated over Bjorken 



x > x mm , y > 0.1, and over the regions Q' > 5 GeV and x' > x' min as indicated. 

section. A negative result however cannot be turned against the theory, it would 
rather restrict the unknown behaviour of F(x') at small x'. Most promising is the 
kinematic region of small Bjorken-x, because both the total DIS cross section and 
the predicted fraction of instanton induced events increase towards small x (see 
fig. Gjjb). 



9.4 Experimental Signatures 

In the theoretically safe region, x' > 0.2, the expected fraction of instanton events 
in all DIS events is of O(10~ 3 — 10~ 4 ) (compare fig.p.9|b), too small to be detected 
in inclusive cross section measurements (i.e. the structure function F 2 ). Instead, 
dedicated searches for the characteristic features of instanton events in the hadronic 
final state have to be performed. A Monte Carlo generator (QCDINS [|311|| ) to sim- 
ulate the hadronic final state of instanton events in DIS is available. In general, the 
event shape predictions are more stable than the rate predictions, because poorly 
known factors cancel. The instanton event properties can be contrasted with pre- 
dictions from event generators for normal DIS events (ARIADNE [K(J , LEPTO 



and HERWIG which give an overall satisfactory description of the DIS final 



state properties 123 



In the q*g rest frame 2nj — 1 quark and antiquarks and n g gluons are emitted 
isotropically from the instanton subprocess. n g is Poisson distributed with ||307| , |309 



2?r ,AF(x') 
( n ~ — x > i _ x > —A/. (9.1 

a* dx' 
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After hadronization, this leads to a spherical system with a high multiplicity of 



hadrons, depending mainly on the available centre of mass energy y s' = Q'yl/x' — 1 . 

For a typical situation (V = 0.2, Q' = 5 GeV =>• sfs' = 10 GeV), (n g ) = 0(2). 
About ra p = 10 partons and n = 20 hadrons are expected. The expected par- 
ton momentum spectrum is semi-hard ||304jl with transverse momentum (jpr) ~ 
(7r/4)(^/(n p ». 

Hadronic final state properties are conveniently being studied in the centre of 
mass system (CMS) of the incoming proton and the virtual boson, i.e. the CMS of 
the hadronic final state. Longitudinal and transverse quantities are calculated with 
respect to the virtual boson direction (defining the +z direction). The pseudora- 
pidity 7] is defined as rj = — lntan(0/2), where 9 is the angle with respect to the 
virtual photon direction. When boosted to the CMS, the hadrons emerging from 
the instanton subprocess occupy a band in pseudorapidity of half width Ar) ~ 1, 
which is homogeneously populated in azimuth ||304j . 

The characteristics of instanton events by which they can be distinguished from 
normal DIS events are therefore: high multiplicity with large transverse energy; 
spherical event configuration (apart from the current jet); and the presence of all 
flavours (twice!) that are kinematically allowed in each event. One would therefore 
look for events which in addition to the other characteristics are rich in K°s, charm 
decays, secondary vertices, muons etc.. In general, the strength of instanton signals 
in the hadronic final state increases somewhat towards low x' and large Q' 2 due to 
the increasing "instanton mass" a/s 7 = Q'Jl/x' — 1. 




Figure 9.7: Transverse energy Et in the hadronic CMS ]312[ . a) The Et flow vs. -q. The proton 
remnant direction is to the left. The standard QCD model (nDIS=ARIADNE) and different 
instanton scenarios are confronted with the HI data [114]. The excluded scenario J312[ | with an 
instanton fraction // > 11.8% for x 1 > 0.2 is indicated, b) The Et distribution, where the 
transverse energy is measured in the CMS rapidity bin < rj < 2, for two instanton scenarios, and 
the standard QCD model (nDIS= ARIADNE). The plots are normalized to the total number of 
events N that enter the distributions. 
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The "instanton band" shows up in the flow of hadronic transverse energy Et as 
a function of 77 (fig. |9~77|a). It's height and position depends on x' and Q' 2 (and also 
on x and Q 2 ). In normal DIS events on average an Et of 2 GeV per 77 unit is 
observed. In instanton induced events, that number may go up to 10 GeV for low 
x'. A possible search strategy could involve the Et distribution in a selected rapidity 
band (fig. |9.7|b), looking for high E T events in the tail of the distribution [ |313| , |312| |. 

Further discrimination can be obtained ||313|| from the fact that for instanton 



events the Et should be distributed isotropically, while normal DIS events are jet- 
like, in particular for large Et ■ One defines 



E n 



mm 



n\ 



E\ r 



n 



(9.9) 



The sum runs over all final state hadrons i with momentum fy. n is the unit vector 
perpendicular to the virtual photon axis which minimizes E ont and thus defines the 
event plane, n' lies in the event plane and is normal to both n and the virtual 
photon axis. It is easy to show that for an ideal isotropic "instanton decay" , E out = 
a/s'/2 [ 513 ]. The "instanton mass" y/s' can thus be reconstructed experimentally 
(fig. |9T8|a). Normal DIS events, either "1+1" or "2+1" jet events (the +1 refers to 
the unobserved proton remnant) are contained in the event plane, E out <C E in , in 
contrast to instanton events with E ont pa E\ n (see fig. 
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Figure 9.8: a) The correlation between 2-E' out and the "instanton mass", Wj = yP (top), and 
the resolution for y/~P that can be achieved (bottom) |313| . The primes indicate additional cuts in 
77 to minimize higher order QCD radiation which may wash out the relation between -E ut and sfp. 
b) E' out vs. E' in for normal (top, HERWIG) and instanton induced events (bottom, QCDINS) 
|313| . Both distributions are taken in the hadronic CMS for events with 0.001 < x < 0.01, 
0.1 < y < 0.6 and 20 GeV 2 < Q 2 < 70 GeV 2 . 

Instanton events are characterized by a large particle density localized in rapidity. 
In normal DIS events there are about 2 charged particles per unit of pseudorapidity 
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118||, rather uniformly distributed in 77. For a low x' cut-off, that number goes up to 



10 in the peak of the instanton band ||312|| . Very sensitive to instanton events is the 
charged particle multiplicity distribution ||312|| , see fig. |9~9j a. A significant fraction of 
the instanton events would lead to charged multiplicities which are very unlikely to 
be found in normal DIS events. Furthermore, particle-particle correlation functions 
should be influenced by instanton effects | |3 14| . 
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Figure 9.9: a) The probability distribution P{n) of the charged particle multiplicity n from 
the CMS pseudorapidity range 1 < 77 < 5 for events with 185 GeV < W < 220 GeV. Shown 



are the unfolded HI data [116 , the expectation from a standard DIS model (nDIS=ARIADNE)), 
and the predictions for instanton events with different cut-off scenarios [312 . b) The maximally 
allowed fraction fn m of instanton induced events in DIS for Q' 2 > 25 GeV 2 and x' > 0.2 from 



transverse energy flows and the multiplicity distribution as function of x [312|. Regions above 



the lines are excluded at 95% C.L.. The numbers give the average Q 2 values in GeV 2 for the 



bins. The theory prediction, calculated with QCDINS |TT), for 10 GeV 2 < Q 2 < 80 GeV 2 is 
superimposed (full line, label INSDIS). 



9.5 Searches for Instanton Processes 

The fact that instanton events look very different from the expectation for stan- 
dard QCD events can be exploited to search for instanton signals in the HERA 
data. One strategy is to compare the shape of hadronic final state distributions to 
the expectation from standard QCD events (nDIS) with an admixture of instanton 
events (INSDIS) of fraction // . In case the measured distribution agrees with the 
standard QCD expectation, a limit on the fraction of instanton induced events in 
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DIS fi < /ii m can be set. The caveat of this method is that one has to make an 
assumption on what standard QCD looks like. In particular at small x that issue 
is under debate ||278| , |246| , [315|| , see section [7l| and chapter |[ There exists a dan- 
ger that an instanton effect is tuned or explained away by stretching the standard 
QCD predictions by generator tuning, introducing BFKL effects etc.. A good un- 
derstanding of standard QCD will therefore be crucial for the positive identification 
of instanton effects. 

In the first search for instanton events [ |117| ] an anomalous K° yield has been 
looked for. For x > 10 -3 about 0.12 K° mesons (including K°) have been measured 
per event and unit pseudorapidity, with a relatively flat 77 distribution. For instanton 
events with x' > 0.2, a peaked distribution with about 0.55 K° per event and unit rj 



is expected. From the comparison with standard QCD event generators |56], gl], HI 
and the instanton generator ||311|| an upper limit of /j < /n m 
obtained 



6% at 95% C.L. is 



HZ I- 



The charged particle multiplicity distribution P(n) in high energy reactions can 
often by described by a negative binomial distribution (NBD). Also the DIS data 
are relatively well described by NBDs [TT£\. The multiplicity distribution from the 
CMS interval 1 < i] < 5 for events with W = 80 — 115 GeV (corresponding to 
x > 0.0007) can be parametrized with an NBD of mean (n) = 6.90 ± 0.33. Possible 
deviations from an NBD allow for an instanton fraction of at most // =2.7% at 95% 
C.L. 
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Other measured event shapes have been systematically analysed in terms of their 
sensitivity to instanton events ||312|| , and their dependence on the kinematic variables 
x, Q 2 , x f , Q' 2 . The most sensitive distributions were the transverse energy flows ||114 



the pseudorapidity distribution of charged particles and their p T spectra [|118|| . For 
example, the flow has been measured over a wide range of x and Q 2 , allowing 
to extend the search region down to x — 0.0001. From a shape analysis ||312|| (see 
fig. |9.7| ), instanton fractions // between 5 and 13 % can be excluded for x' > 0.2 (see 
figs. |9.9| b, |9.10p . For lower x' the signal is more prominent, and somewhat better 
limits are obtained. 

The fact that HI ||116|| did not observe any events above a certain multiplicity 
n max has been exploited |[312 j to place more stringent limits on instanton produc- 
tionF]. A significant fraction of instanton induced events would have multiplicities 
n > n max (compare fig. p.9| a). Instanton fractions // > 0.4 — 0.6% can therefore 
be excluded for x' > 0.2 (see figs. |9.9|b, |9.10| ), and somewhat lower // values for a 
lower cut-off x' > 0.1 ||312|| . This search method has the advantage that, in contrast 
to the previous shape comparisons, it does not rely on assumptions for standard 
QCD event topologies, since no background needs to be subtracted. Unavoidable of 
course is the dependence on the expected instanton event shape, which may be even 
more uncertain than the standard QCD event shapes. 

The available bounds on instanton production are summarised in tab. 0. The 

were derived from the shape of the 



2 The previous limits from the HI multiplicity analysis |11 
multiplicity distribution for n < n max • 
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Figure 9.10: 95% C.L. limits on instanton production with Q' 2 > 25 GeV 2 and x' > 0.2 
|312| . The cross-section limits (cr lim ) together with the maximally allowed instanton fraction /u m 
are shown in the x -Q 2 plane. They are obtained from the flow analysis (open fields), and 
from the multiplicity analysis (shaded fields) with their numbers at the right edge. The boundaries 
implied by the analysis cuts of the energy flow analysis in the angle and energy of the scattered 
electron, 9 e < 173° and E e > 12 GeV, and by the requirement W 2 > 4400 GeV 2 are indicated. 



most stringent limits for the theoretically "safe" scenario x' > 0.2 are still a factor 
20 higher than what is predicted from the instanton theory, see fig. |9.9| b. Limits for 
other scenarios can be found in |[312j| . For x' > 0.1 they are already below the naive 
extrapolation into the theoretically uncertain region, providing a constraint for the 
theory and the holy grail function F(x'). 



9.6 Conclusion 

Instanton transitions, a yet unexplored facette of non-abelian gauge field theories, 
have been discussed. While in the electroweak theory the B + L violating effects in- 
duced by instantons are expected only at energies > 10 TeV , their chirality violating 
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analysis 


DIS kinematics covered 


instanton scenario 


limit 




Q 2 (GeV a ) 


X 


W (GeV) 


Q' 2 (GeV") 


x' 


/lim 


K u 11171 




10-70 


0.001 -0.01 


95 - 230 


> 1 


£ 0.2 


6 % 


multipl. flllGf 


10 - 80 


0.0007 - 0.012 


80 - 115 


> 1 


> 0.2 


2.7% 


E T flows |312| 


5-50 


0.0001 - 0.01 


65 - 230 


> 25 


> 0.2 


5 - 13 % 


multipl. |3 


if 


10-80 


0.0001 - 0.01 


80 - 220 


> 25 


> 0.2 


0.4 - 0.6 % 



Table 9.1: Limits on QCD instantons in DIS. A fraction fj > f\ im of instanton induced events 
in DIS is excluded at 95% C.L.. 



pendant in QCD could lead to striking signatures already at present day colliders. 
In DIS at HERA, these are a high particle multiplicity with large transverse energy 
localized in rapidity, and s, c and possibly b quarks in the final state. 

The expected contribution to DIS events from instantons is of 0(1CT 3 — 1CT 4 ), 
with substantial theoretical uncertainties. First analysis of HERA data taken in the 
years < 1994, corresponding to an integrated luminosity of (9(1.3 pb _1 ), are still a 
factor Rj 20 above the prediction. With higher statistics data samples (O(40 pb _1 ) 
up to the end of 1997) and improved search strategies, a fundamental discovery at 
HERA appears to be in reach. 

As the study of the hadronic final state at HERA has shown, model predictions 
for normal DIS final states are not unique. It will probably be safe to take some 
uncertainty into account also for the prediction of the instanton final state. A good 
understanding of QCD in DIS will be vital for instanton searches. 

It might be possible to exploit also other reactions than DIS, such as photo- 
production, where the hard scale needed for reliable instanton calculations could be 
provided by the pt of a jet. Instantons have even be proposed ||316|1 as an explanation 
of the possible excess of high Q 2 events at HERA 
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Chapter 10 
Summary 



HERA has entered a new kinematic regime in deep inelastic scattering: at small 
x, large W and both small and large Q 2 . After five years of HERA operation, the 
hadronic final state has been measured in great detail. The large body of data has 
been confronted with QCD, the theory of the strong interaction. 
We have seen progress happening at different stages: 

• QCD predictions are tested experimentally, and parameters of the theory are 
being extracted. An example are the measurements of jet rates and the ex- 
traction of a s or of the gluon density in the proton. 

• Where the theory is less well established, there is an interplay between ex- 
perimental data and theory, in which new ideas are born and developed (or 
rejected). Examples are soft colour interactions as an explanation for rapidity 
gaps and large E T flows, or power corrections for hadronization effects. 

• The data are not yet understood. Either there exist no predictions, or there 
are conflicting explanations. An example are the measured forward jet rates, 
which could be explained by BFKL evolution, or by hadronziation effects, or 
by a resolved structure of the virtual photon. 

The availability of as much data as possible over a wide kinematic range is very 
important for the scientific process to work efficiently A new idea which explains 
one aspect of the data can be rejected early if it fails in another. The data have 
severely constrained QCD models for the hadronic final state which varied widely 
before the data had become available. 

General event properties in deep inelastic scattering at HERA 

In the hadronic CMS, the produced hadrons cover about 5 units of rapidity in either 
hemisphere, depending mainly on the available CM energy of the hadronic system, 
W . There are about 2.5 charged particles per unit rapidity, and the transverse 
energy is on average 2 GeV per unit rapidity. Multiplicity distributions exhibit the 
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familiar KNO scaling, and Bose-Einstein correlations are measured similar to other 
reactions. A large rapidity gap is found in about 10% of the events. Around 25% 
of the events contain charm. In the current hemisphere about 0.9 K° mesons and 
0.15 A particles are measured per event. In 5 to 10% of the events more than one 
jet is found (apart from the remnant). 

The global event features are largely determined by the available phase space 
given by the CM energy W of the hadronic system. The dependence on the photon 
virtuality Q 2 has also been studied. In the current region the transverse activ- 
ity (Et flow, px spectra, jets) increases with Q 2 . In the central rapidity region 
Q 2 dependencies are small. Here the data are similar to what has been measured in 
photoproduction (Q 2 =0) or hadron-hadron collisions. The particle density per unit 
rapidity increases slower than in e + e~ annihilation with the available CM energy 
W. When compared with DIS data at lower W, a faster than logarithmic rise of the 
average multiplicity is observed. 

Comparisons with QCD 

Transverse energy and momentum distributions provide evidence for hard QCD 
radiation beyond the LO matrix element. More Et has been measured than antici- 
pated before HERA data taking. It is not yet clear whether soft colour interactions 
are required to describe the transverse energy flow, or whether not yet well estab- 
lished perturbative QCD effects are responsible for the large E T measured. Soft 
colour interactions are of interest also because they provide a mechanism for pro- 
ducing rapidity gaps. Alternatively, colour dipole radiation or BFKL evolution can 
be employed to produce as much Et as seen in the data. 

The scaled charged particle momentum spectra measured either in the current 
hemispheres of the hadronic CMS with W as scale or the Breit frame with Q as 
scale are in general well described by QCD models in various approximations of the 
perturbative shower evolution in conjunction with non-perturbative hadronization 
models. Where applicable, perturbative QCD calculations in NLO employing frag- 
mentation functions measured elsewhere are in agreement with the data. Scaling 
violations of the spectra will provide a means to measure a s . 

The growth of multiplicity with energy is described by perturbative QCD calcu- 
lations in the MLLA approximation, assuming local parton hadron duality (LPHD). 
The shape of the momentum spectra in the "hump backed" form is approximately 
Gaussian. It's evolution with Q is consistent with the MLLA+LPHD expectation 
assuming soft colour coherence. Models without soft colour coherence describe the 
data less well. When allowance is made for the boson-gluon fusion process in DIS, 
the particle spectra are similar to what has been measured in e + e~ annihilation. 
Evidence is found for the prediction that the soft part of the Lorentz invariant spec- 
trum is independent of the energy of the fragmenting parton. The assumption of a 
constant coupling a s is in stark conflict with the data. 

Global event shape variables like thrust, jet mass and jet broadening in the 



148 



Breit frame current hemisphere show an increasing event collimation with energy 
Q. "Power corrections" supplementing the NLO calculations provide a good de- 
scription of the data. The power corrections oc 1/Q parametrize higher orders and 
hadronization effects and appear to be universal within 20% - the same normaliza- 
tion parameter can be used for different shape variables and ep as well as e + e~ re- 
actions. The concept of "power corrections" may turn out to be very fruitful for 
the understanding of hadronization effects. Assuming power corrections, the strong 
coupling a s has been determined with competitive precision. 

Of similar precision is the a s determination from the rate of 2+1 jet events. The 
strong coupling a s has been measured as a function of Q 2 from 2+1 jet events in a 
region where NLO QCD is able to describe the data, namely where Q 2 is not too 
small. The data are consistent with the expected Q 2 dependence of the running 
coupling. However, the optimal strategy for a s measurements from jets has not yet 
been found. A better understanding of the connection between measured hadron jets 
and perturbative QCD predictions for parton jets appears to be needed. Smaller 
theoretical uncertainties are anticipated with other jet algorithms and increased 
statistics at higher Q 2 or jet pr- 

At small x and Q 2 and in particular towards the proton remnant the measured 
jet rates are larger than expected from either NLO calculations or conventional QCD 
models. "Conventional" stands for deep inelastic scattering, where a virtual photon 
scatters on a parton in a proton, and where the partons in the proton obey DGLAP 
evolution with Q 2 as scale. A good description is obtained either by the colour 
dipole model, or by a model in which the partonic content of the virtual photon is 
resolved by the large jet pt ■ 

A strong growth of the gluon density towards small x has been inferred from 
DGLAP analyses of structure function measurements. This is confirmed by direct 
measurements from jet and charm production, by which photon-gluon fusion pro- 
cesses can be tagged. The errors are still large, however. It will be very important 
to extend such direct measurements towards smaller x values. 

Small x physics has received considerable interest. It has been speculated that 
the rise of F 2 with decreasing x is a sign for BFKL evolution, but DGLAP evolution 
was found to be able to account for the data as well. Traces of BFKL evolution have 
therefore been searched for in less inclusive observables in the hadronic final state, 
namely enhanced transverse activity in between the current and the target frag- 
mentation regions: transverse energy flow, high px particles, "forward jets" . The 
data could not be described by conventional DGLAP evolution. Though BFKL 
calculations can explain the measured enhanced activity, there exist competing ex- 
planations, again making use of a virtual photon structure. It remains to be seen 
though whether or not such an alternative model would be in conflict with other 
HERA data. In case both explanations work equally well, one might speculate that 
they have more in common than is apparent. 

Finally, the possibility exists to discover QCD instanton effects at HERA. Such 
a discovery would open a new field of non-perturbative QCD, and would imply 
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also baryon and lepton number violation in the electroweak sector already within 
the standard model. Observable effects are predicted for the hadronic final state, 
and first bounds on instanton production have been placed. The predicted cross 
sections are such that it appears likely that with increased sensitivity instantons 
can be discovered at HERA, or the theory in its present form can be rejected. 
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